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utilizing a library double precision numerical integration 
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to provide the user with a highly accurate time response of 
a system which has been designed on the OPTSYSX Frogram. 
This series of programs permits the user to rapidly design, 
analyze and test all types of Optimal Systems Control prob- 
lems. Examples of the various types of problems are worked 
through to illustrate all of the capabilities available. 
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I. IBTBODDCTION 

The purpose of this thesis is to describe and demon- 
strate the modification and additions to the existing 
FOETEAN program (0P1SYSX) which is used in the study, 
design, and application of Optimal Systems Control theory. 

The CPTimal systems control program (OPTSYS) was origi- 
nated in 1971 by Hall [Ref. 1] to support his research in 
rotary-wing aircraft control systems. The most recent 

program modifications were made by Walker [Ref. 2] and Liu 
[Ref. 3 ] of Stanford University and are designated OPTSYS 4 
and CEISYS 5, respectively. The OPTSYS modifications made 
by Hoden [Ref. 4] were primarily devoted to creating a 
user-friendly interactive version (OPTSYSX) of the OPTSYS 4 
Program . 

The goal of this thesis work was to develop a program 
set which will operate in an interactive mode and plot the 
time response of a State Variable Control System which has 
been developed and/cr analyzed using the OPTSYSX Program. 
!1inor modifications tc the OPTSYSX Program were necessary to 
allow the user to build a data file of the matrices required 
for the time response calculations plus additional matrices 
which cculd be used again upon reentry to the OPTSYSX 
Program without the laborious and time-consuming task of 
reentering each data element in all of the required 
matrices . 

It is assumed that the reader/user is' familiar with the 
basic concepts of Control Theory and Optimal Systems Design. 
The symb cl/naming conventions of Bryson [Ref. 5] are used in 
the program operation discussion and in the examples of the 
problems solved using this system. A glossary of the 
symbols and abbreviations used in this discussion is 
provided on page 8. 
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An explanation of the OPTSYSX capabilities and a program 
set overview are presented first. 

This work concludes with examples of varioas types of 
problems demonstrated in the interactive mode, including a 
copy of each terminal session with the final results. A set 
of complete program listings are included in Appendices A, 
E, C an d D . 
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II. THE OPTSYSX COMPOTER PROGRAM 



A. GENERAL BACKGROO HD 

OPTSYSX is an interactive, double-precision FORTRAN 
program employing modern control theory analysis techniques. 
Its extensive capabilities include the synthesis and anal- 
ysis of filters and regulators as well as eigensystem anal- 
ysis, modal distribution, transfer function analysis and 
power spectral density computations. The modifications to 
the OPTSYS Program introduced by this thesis work have not 
affected any of the program’s original capabilities. 

B. PROGRAM OVERVIEW 

OPTSYSX is an extremely large and complex program with 
over 2800 lines of code. In order to use this program in 
its small (set up for 32X32 matrices) version, the user is 
required to extend the IBM 370/3033 virtual machine (VM) 
memory capacity beyond 720 kilobytes which is ths default VM 
memory size. A significant increase in the size of the 
OPTSYSX program would make the program too large to operate 
on a cne megabyte VM, the largest virtual memory available 
or. a user's virtual machine. The high resolution plotting 
software is limited to single precision variables. 
Therefore double precision library routines cannot be called 
from the plotting program. For these reasons the task of 
obtaining the time response of a system was divided between 
three programs, OPTSYSX, 0PTCA1C and OPTPLOT. An Executive 
program (OPTSYS EXEC) was written to make the interfacing of 
the three programs transparent to the user. 

Minor modifications were made to the OPTSK SX program 
including the addition of three subroutines to handle the 
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input and output of matrix data to and from a data file on 
the user’s disk. The OPTCALC program performs the double 
precision numerical integration of the system Df equations 
over time and creates another disk data file of the state 
variable variation with time. OPTPLOT takes this time 
response data and presents it in a graphical format on the 
TEK 618 graphical display or as a VERSATEC pen plot. 

C. OPTSYS EXEC 

The OPTSYS EXEC is written in the EXEC 2 language. This 
language allows the EXEC to issue almost any command that 
can be entered in the direct mode at the terminal. Therefore 
an EXEC is the ideal controller for the "black box" type of 
system where the user is not aware of what is actually 
taking place within the program (s). The OPTSYS EXEC was 
written to complete all of the required interfacing between 
the three programs (OPTSYSX, OPTCALC and OPTPLOT) , without 
the direct guidance or control of the user. By answering 
guestions presented on the terminal screen, the user deter- 
mines the logic flow through the EXEC while the EXEC estab- 
lishes the appropriate FILEDEFs and loads the programs 
required by the user’s desires. 

D. OPTSYSX MODIFICATIONS 

Three subroutines (RDMATF, RDMAT, and FRTMATi were added 
to the CPTSYSX Program for data file read error check, 
matrix input from a data file and matrix output to a data 
file, respect ively . „ These three subroutines provide the 

user with the opportunity to save the [ F ], [G], [ H], 

[GAMMA], [A] and [B] matrices for use in a subsequent run of 
the OPTSYSX Program. The WRTMAT subroutine also saves the 
[C] and [K] matrices for system time response calculation 
and plotting by the OPTCALC and OPTPLOT programs. 
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1 . RDM ATF S ubro utine 

The P.DKATF subroutine is used to check for the exis- 
tence of a previously generated file containing matrix data. 
Seven flags may be set by this subroutine. Six of these 
flags correlate with the six matrices that the user may save 
for reuse later in the OPTSYSX Program. The remaining flag 
(IRDMAT) must be set to enable the RDMAT subroutine to read 
matrix information from the data file. A READ statement of 
the form 

READ (9,11 1 , ERR=2 22, END=3 33 ) A, 3 

(where "111", "222" and "333" designate line numo ers for the 
FORMAT statement and branch on ERROR or branch on 

END-CF- FILE routines, respectively) is used for the data 

file check. The ncnexistence of the file or premature 
END- OF- FILE are detected by the ERR and END o hecks which 
cause a branch to a routine that sets the IRDMAT flag to "0" 
and returns to the calling program. 

If no error is detected during the initial read 
attempt, the variable B is checked for the sentinel "1". 
This second check is to help ensure that the file is actu- 
ally a file which contains valid matrix elements. The user 
is then presented a message which asks if he/she wants to 
use the matrices which are available. The user may respond 
with cn cf three answers: 

1 . Use all of the matrices. 

2. Use selected matrices. 

3. Use none of the matrices. 

If the answer is "1" or "2", the subroutine reads the matrix 
dimensions (Ns, Nc, No and Ng) from the data file and 
changes the IRDMAT flag to 1 to key the RDMAT subroutine to 
read the matrix elements from the data file. If the answer 
is "1" all of the matrix-save flags are set to "1". 
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matrices. When the user has completed the analysis/design of 
the system of interest, the WRTMAT subroutine asks the user 
if he/she wants to calculate the time response of the system 
which the user just designed. If the user answers YES, the 
WRTMAT subroutine generates a data file of appropriate 
matrix information and halts execution of the OPTSYSX 
program. Control then reverts to the OPTSYS EXEC. If the 
user answers NO, the WRTMAT subroutine returns control to 
the main program and normal OPTSYSX program operation 
continues. 

The information v T ritten to the data file consists of 
2 "1"s (which are used as a sentinel or flag by the RDMATF 
subroutine (as previously explained) and in a similar manner 
by the OPTSYS EXEC), followed by the matrix dimensions (Ns, 
Nc, Nc and Ng) and then by the [ F ], [G], [H], [GAMMA], [C], 
[K], [A] an< ^ [B] matrices. These matrix elements are written 
to the OPTMAT DATA file using a 4D20.13 format is a compro- 
mise between the maximum feasible accuracy of data exchange 
between the double- precisi on programs and the use of a 
moderate amount of the user’s disk space. 

E. OPTCALC PROGRAM 

1 • System I nt eg rat ion 

The OPTCALC program is a FORTRAN interactive 

double-precision system integration routine. This program 
uses the International Mathematical & Statistical Library 
(IMSL) subroutine DG EAR to perform the numerical integration 
of the system under analysis. The stiff system mode of DGEAR 
is used in order tc provide the capability to do time 

response calculations of the X-29A longitudinal axis back-up 
mode system which is an 98 X 98 stiff system. 
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2 . System Equati on Rep res e ntatio n 
The OPTCALC program uses the state variable format such as 

xdot = [F]*x + [3]*uc (2.1) 

to define the system. In this system the [F] matrix is the 
open-loop dynamics matrix (system or plant) and the [G] 
matrix is the control matrix. The variable assignments are 
x as the state vector and uc as the control input vector. It 
follows that xdot is the time derivative of x. 

Various forms of equation 2.1 are used for all the 
time response calculations. The [ F] matrix is modified to 
[F+G^C] 1 for closed-loop (regulator only) system calcula- 
tions as in equations 2.2 and 2.3. 

Xdot = [F + G*C]*x + [G]*uc (2.2) 



u = [ C ]*x + uc (2.3) 

For this closed loop system, the £C] matrix is the control 
gain or regulator gain matrix and u is the total input 
vector. 

The combined filter and regulator systems can be 
represented by equations 2.2, 2.4, 2.5, 2.6 and 2.7. The 

[H] matrix is the measurement scaling (observer output) 
matrix and the [K] matrix is the estimator or Kalman filter 
observer gain matrix. The variables xe, xedot and x are the 
state estimate vector, the derivative of the state estimate 
vector and the state reconstruction error, respectively. 



1 The OPTSYSX sign convention for the C matrix is the 
negative of the standard normally used in controls. 
Therefore [F+G*C] has the correct sign for OPTSYSX matrices. 
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xdot = [F + G*C]*x + [G]*uc 



( 2 . 2 ) 



2 = [ H ]*X 



(2.4) 



xedot = [ F J*xe + [G]*u + [K]*(z - [ H]*xe) (2.5) 



u = [ C ]*xe 



( 2 . 6 ) 



x = x - xe 



(2.7) 



Equations 2.2, 2.4, 2.5, and 2.6 can be combined into the 

augmented matrix form of equation 2.8. 



Equation 2.8 is an augmented equation in which the Ns dimen- 
sion has been doubled and the state and state estimate 
vectors have been combined into one vector of 2*Ns length. 
The [G] matrix is also augmented by repeating the first Ns 
rows again beginning at row Ns+1, making the new [G] matrix 
dimensions (2*Ns,Nc). 



sented by the same equations with the [C] matrix set to 0. 
These equations are: 



xdot 



F 




( 2 . 8 ) 



xedot 



K*H 



In a similar manner, a filter only system can repre- 



xdot = [F]*x + [G]*uc 



( 2 . 1 ) 



2 = [ H ] *X 



(2.4) 
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xedot = [F]*xe + [G]*u + [K]*(z - [ H ]*xa ) 



(2.5) 



x = x - xe (2.7) 

The combination of eguations 2.1, 2.4, 2.5 and 2.7 into an 

augmented [F] matrix is similar to equation 2.8 but with 
fewer terms and the upper right quadrant equal to zero. 
This filter-only augmented system equation is shown below. 



xdot 




xedot 





F 

K*H 




(2.9) 



a. System Selection 

When the CPTCALC program is run, the |F|, |G| 

and 1 C | matrices (and |H| and |K| matrices, if available) 
are presented on the terminal as a check and as a reminder 
of the characte rist ics of the system that has been passed 
from OPTSYSX to OPTCAIC. The user is requested to select the 
type cf system response to be calculated. 

1. CPE N 1C0P TIME RESPONS E. 

2. C10SED LOOP TIME RESPONSE. 

3. CIOSED LOOP FILTER ONLY TIME RESPONSE. 

4. CLOSED LOOP FILTER + REGULATOR TIME RESPONSE. 

Selection of 2, 3 or 4 forms the appropriate system matrix 

equations 2.2, 2.9 or 2. 8 , respectively and doubles the 

length of the |G| matrix, if required. 

b. Defining Calculation Limits and Inputs 

After the user determines the type of system 
under study, the OPTCALC program prompts for the integration 
start and stop times and the number of data points desired. 
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The user has some control over the tradeoff between curve 
fidelity and computer time used by varying the number of 
data points calculated. Computer time use is normally a 
factor only on very large systems. If less than 200 points 
are calculated, the CPTPLOT program uses a curve smoothing 
function which may cause minor inaccuracies in the plotted 
curve but avoids the sharp peaks and irregular appearance 
generated by plotting straight lines between an insufficient 
number of data points. When 200 or more points are calcu- 
lated no smoothing is done. The points are connected by very 
short straight lines which has the appearance of a smooth 
curve. 

Step or ramp functions are available as control 
inputs. Only one type of function can be used for each 
control, but the magnitude and start and stop times can be 
set as desired. 

The state and state estimate initial conditions 
can be set to any value by the user. The control inputs and 
non-zero initial conditions can be used simultaneously, if 
desired. Eefore the time response calculations begin, the 
user is given the opportunity to make changes in any area of 
the system integration initial conditions that have been 
previously selected. 

3 • S ystem T im e Re spon se 

Equation 2.1 is evaluated directly in the open loop 
system response calculations. The FCN subroutine was 
written to evaluate the system of equations for the DGEAP. 
IMS1 subroutine. Each time the FCN subroutine is called by 
EGEAF., it updates the control inputs (ac) and then evaluates 
each state derivative by summing all the terms across that 
row of the [F]*x and [G]*uc matrices. The same FCN subrou- 
tine is used for all system integrations. As explained in 
the previous section, the [F] matrix is replaced by the 
[F+G*C] matrix for closed loop system response problems. 
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a. Systems With Filters 

Augmented equations 2.9 and 2.8 are ased for the 
time response evaluation of systems with filters only and 
systems with filters plus regulators, respectively. The 
augmented matrix is developed as a dummy matrix and is then 
inserted as the [F] matrix with the dimensions doubled 
(2*Ns). The [G] matrix is also augmented by repeating the 
first Ns rows again beginning at row Ns + 1, macing the new 

[G] matrix dimensions (2*Ns,Nc). 

The augmented system can be evaluated by simply 
doubling the old system row and column dimension (Ns) and 
calling the DG EAR integration subroutine. Using this 
method, the existing FCN subroutine requires no changes to 
evaluate the augmented system. 

b. OPTCALC Output 

The OPTCALC program uses FILEDEF 8 for the data 
file output as well as FILEDEF 5 to read and write to the 
terminal. The output data file contains the following 
discrete information: the matrix dimensions Ns and Nc, the 

augmented matrix dimensions, the number of data points 
calculated and a flag to indicate that an augmented matrix 
was calculated. The [C] matrix is passed to permit the 
calculation of u the total control input to the system. The 
final portion of the data file is individual data points of 
time, external control input (uc) and each state (x) and 
state estimate (xe) . This data file provides all of the 
data required by the OPTPLOT program to make a smooth graph- 
ical response curve. 



F. OPTPLOT PROGRAM 

The OPTPLOT program is a FORTRAN interactive plotting 
program using the Display Integrated Software System and 
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Plotting Language (DISSPLA) by Integrated Software Systems 
Corporation. This program provides the user a high resolu- 
tion graphical display of the system’s time response ani if 
desired will provide a VERSATEC pen-plot of the same graph. 

1 . Gen era l Oper ation 

Plotting data is received from the program OPTCALC 
via a data file on FILEDEF 8. The types of data provided in 
the file are discussed in the previous section. The program 
presents the user with a series of questions to determine: 



1. The number of curves to plot. 

2. Select the type cf variable for each curve. 

3. Select the variable subscript for each curve. 

4. Select the number of headings and contents of each. 



The program then plots the selected variables and provides 
the user a graphical display on the Tektronics 618 (TEK 618) 
display. 

The following Main Menu is then presented which 
provides the user with the major decision points of the 
program. 

1. BEGIN NEW GRAPH CF OTHER CONTROLS , STATES, OR ESTIMATES. 

2. REPLOT PREVIOUSLY SAVED GRAPH DATA. 

3. EDIT THE CURRENT GRAPH. 

4. PLOT REVISED GRAPH ON THE TEK618. 

5. QUIT AND/OR MAKE METAFILE OF THE CURVES 

PREVIOUSLY SAVED. 

The purpose of each of the selections is self- 
explanatory, however the methods of their use may not be. If 
number 3 is selected, the user is then presented the 
following Edit Menu of items to make additions, deletions or 
corrections to the curves that are plotted on the TEK 618 
screen: 



22 



1. CHANGE VARIABLES OR ADD A CURVE ON THE CURRENT PLOT. 

2. DELETE CURVE FROM CURRENT PLOT. 

3. EDIT CURVE TITLE (S) . 

4. EDIT PAGE HEADING (S). 

5. CHANGE THE Y-AXIS SCALE. 

6. CHANGE THE TIME AXIS SCALE. 

7. CHANGE THE PLOT SIZE. 

8. CHANGE THE LETTERING HEIGHT. 

9. CHANGE POSITION OF THE LEGEND. 

10. EDITING COMPLETE. 



This extensive list of modification capabilities provides 
the user with the tocl to make almost any imaginable altera- 
tion to an existing plot. 

Since the OPTFLOT program receives the time response 
data from the OPTCALC program, item 6 of the Edit Menu 
cannot be used to expand the time scale beyond the time span 
previously calculated. Therefore the time axis change 
feature can be used only to select a subset of the original 
data. 



2. VERSATEC Pen -plo ts 

The VERSATEC pen- plots are provided through the 
DISSPOP portion of DISSPLA. In order to use the DISSPOP 
feature, a device independent plot file called a metafile 
must be generated. 

To view the graphical time response plot on the TEK 
618 terminal, the TEK618 option of DISSPLA must be called 
within the OPTPLOT plotting program. The graphical image 
data is then sent to the TEK 618 display screen. 

A metafile is created when the COMPR3 option of 
DISSPLA is called by the OPTPLOT plotting program. When the 
plotting program is executed with the COMPRS option, the 
graphical image data is sent to a metafile on the user’s 
disk. The TEK618 option and COMPRS option ace mutually 
exclusive (only one can be active at a time) , therefore 
graphical data cannot go to both the terminal screen and the 
metafile, concurrently. As a further complication, the 
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TEK61S and COMPRS options ca nn ot be used in alternating 
pattern, first to originate and edit each graph and then to 
add this graph to a metafile possibly containing several 
ether graphs. Therefore if more than one pen-plot per 
terminal session is desired, some type of capability must be 
provided to save the information required to reproduce a 
given graph. 

When the user attempts to leave the current plot 
(ie. selecting items 1, 2 or 5 of the Main Menu) the program 
asks the user to save the current graph data fur later use 
in generating a metafile. This feature provides the capa- 
bility to save any desired graph data in order to later make 
a metafile and obtain a pen-plot. When the COMPRS option is 
used (by selecting item 5 of the Main Menu) , any number of 
graphs may be added to the metafile up to the limit of 
available user disk space (provided graph data has been 
previously saved). After exiting the OPTPLOT program, the 
CPTSYS EXEC asks the user if he/she wants a hard copy of the 
metafile that had been generated during the session. If the 
user answers YES, the OPTSYS EXEC calls the DISSPO? EXEC 
with the VRSTEC option. 

When the user exits the DISSPOP EXEC, the CPTSYS 
EXEC gives the user the options to: 

1 . RUN OPISYS AGAIN. 

2. RUN OPTCALC AGAIN. 

3. Q0II. 

The option to run OPTSYSX again allows the user to use all 
or part of the matrices that had been saved in tae data file 
without manually reentering each element. The OPTCALC 
option could be exercised if the user wants to use the same 
system matrices again, but change the contrul input or 
initial conditions or change the type of system (open. 
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was 



closed, filter only cr filter plus regulator) that 
uated on the previous run. 



eval 
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III. SYS TEM USE AND EXAMPLES 



This chapter contains several basic examples of the four 
types of problems which may be solved using OPTSifSX, OPTCALC 
and OPTPIOT under control of the OPTSYS EXEC. Included with 
these examples are copies of each recorded terminal session. 

A. OPEN-LOOP SYSTEM TIME RESPONSE 

The following open-loop system example was taken from 
[Ref. 6, pp 5.3 - 5.7]. 

The full terminal session is recorded below, with user 
input at the left margin in lower case l etter s or numbers 
below each "?" . 



record on 

BEGIN RECORDING OF TERMINAL SESSION 
E; 1 = 0.01/0.02 1 9: 58 :26 
opt sy s 



THE OPTSYS EXEC CONTROLS A TRIO OF PROGRAMS: 



1. OPTS YS X FORTRAN 

2. OPTCALC FORTRAN 

3. OPTPLOT FORTRAN 



SYSTEM ANALYSIS) 

CALCULATE TIME RESPONSE) 
DISSPLA PLOTTING ROUTINE) 



EACH PROGRAM PASSES INFORMATION TO THE NEXT PROGRAM 
THROUGH A DATA FILE WRITTEN TO THE USERS DISK. IN THIS 
CASE, THESE FILES ARE "OPTMAT DATA" AND "OPTPLOT DATA". 
THE SIZE OF THESE EILES VARY WITH THE SYSTEM ORDER. AND 
CAN USE ABOUT 20£ OF THE USERS DISK SPACE. THEREFORE 
ENSURE THAT SUFFICIENT DISK SPACE IS AVAILABLE. 



- TYPE "E" TO EXIT, ANY OTHER ENTRY TO CONTINUE - 



******:££:**£ * * *** * rfr*******#***************;}:** *** *** * £ 

YOU HAVE A DATA FILE NAMED 'OPTMAT DATA' ON YOUR A 
DISK THAT WAS PREVIOUSLY GENERATED BY THE OPTSYS 
PROGRAM AND CONTAINS THE F, G, H, GAMMA, A AND B 
MATRICES FROM THAT RUN. 

IF YOU WOULD LIKE TO USE THESE SAME MATRICES FOR 
THIS RUN, THE OPTSYS PROGRAM WILL READ IN THE 
DESIRED DATA AT THE APPROPRIATE TIME, 

IF YOU TYPE (Y) ES. 

ANY OTHER INPUT WILL RESULT IN THAT FILE BEING ERASED! 
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************************************************************** 

y 



******************* ******** *********************************** 

DO YOU WANT THE NUMERICAL OUTPUT FROM OPTSYSX TO GO 
TO YOUR TERMINAL S(CREEN) OR TO A D (ISK) FILE? 

(S OR D) 

************** ************************************************ 
S 



************************************************************** 

OUTPUT WILL COME TO YOUR TERMINAL SCREEN. 
************** ************************************************ 

LOADING CPTSYS. . .. 

EXECUTION BEGINS... 

OPTSYSX IS A COMPLETELY INTERACTIVE OPTIMAL SY STEMS CONTROL 
PROGRAM. IT WILL SOLVE NUMEROUS CONTROL PROBLEMS ON THE 
FOLLOWING TYPES OF SYSTEMS CONTROL EQUATIONS; 

X DOT = {F} *X + {G}*U + {GAM J * (W + W 0) 

MEASUREMENT EQUATION — 

Z = {H}*X + {D} *W + V 

REGULATOR PERFORMANCE INDEX— 

J = 1/2 * INTEGRAL (Y *{A}*Y + U * {B} *U) DT 

STATE FEEDEACK GAIN DEFINITION — 

U = - {C} *X 

DO YOU WISH TO CONTINUE? TYPE "YES" CR ''NO". 

y 

— BATA ENTRY — 

ALTHOUGH OPTSYSX IS SPECIFICALLY DESIGNED TO READ 
ALL MATRIX DATA INTERACTIVELY. SEVERAL ALTERNATE 
METHODS ARE AVAILABLE TO USERS: 

METHOD 1— THE "F»,"G", AND "GAMMA" MATRICES 
MAY BE READ FROM SEPARATE DATA FILES. 

METHOD 2— THE "F»,»G", AND "GAMMA" MATRICES MAY BE 
EXPLICITLY DEFINED WITHIN SUBROUTINE "SETUP". 

{NOTE : IN EITHER CASE, THE USER SHOULD OBTAIN A COPY 
OF THE FEOGRAM LISTING AND EXAMINE 
THE EXAMPLES CONTAINED IN S/R "SETUP".} 

DO YOU WISH TO CONTINUE? TYPE "YES" OP. "NO". 

y 



DO YOU WISH TO INPUT THE "F", "G", AND "GAMMA" 
MATRICES FFCM SUBROUTINE "SETUP" IAW THE 
METHOD DESCRIBED ON THE PREVIOUS SCREEN? 

TYPE "YES" OR "NO". 

I) 
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GENERAL OPTSYSX OPTIONS: 

OPTION 1 — SYSTEM ANALYSIS WITHOUT 

OPEN-LOOP EIGENSYS TEM CALCULATIONS. 

OPTION 2 -- SYSTEM ANALYSIS WITH OPEN-LOOP 
EIGEN SYSTEM CALCULATIONS. 

OPTION 3 — OPEN- LOOP EIGENSYSTEM FOUND 
AND PROGRAM TERMINATES. 

{"F"- MATRIX ENTRY FOLLOWS IMMEDIATELY.} 

OPTION 4 — MODAL DISTRIBUTION MATRICES COMPUTED 
WITHOUT FILTER OR REGULATOR SYNTHESIS 
OR STEADY-STATE ANALYSIS. 

SELECT AN OPTION: 1 , 2 , 3 , OR 4. 



DO YOU DESIRE RMS VALUES OF STATE AND CONTROL? 

TYPE "YES" GR "NO", 
n 



r> 

i 



OPEN- ICO? 
OPTION 
OPTION 
CP T ION 
OPTION 
SELECT 



TRANSFER FUNCTION OPTIONS: 

1 — _ NO OPEN-LOOP TRANSFER FUNCTIONS COMPUTED. 

2 — PCLES, RESIDUES, AND ZEROS COMPUTED. 

3 — ONLY POLES AND ZEROS COMPUTED. 

4 — ONLY POLES AND RESIDUES COMPUTED. 

AN OPTION : 1 , 2 , 3, OR 4. 



i 



NCI SE TRANSFER 
OPTION 1 — 
OPTION 2 — 
OPTION 3 — 
OPTION 4 -- 
SELECT AN 



FUNCTION OPTIONS: 

NO NOISE TRANSFER FUNCTIONS COMPUTED. 
PCLES, RESIDUES, AND ZEROS COMPUTED. 
ONLY POLES AND ZEROS COMPUTED. 

CNLY POLES AND RESIDUES COMPUTED. 
OFTION : 1, 2, 3, OR 4. 



COMPENSATOR 
OPTION 1 
OPTION 2 
CPTION 3 
CPTION 4 
{NO 



TRANSFER FUNCTION OPTIONS: 

-- NO COMP. TRANSFER FUNCTIONS COMPUTED. 

— PCLES, RESIDUES, AND ZEROS COMPUTED. 

— CNLY POLES AND ZEROS COMPUTED. 

-- CNLY POLES AND RESIDUES COMPUTED. 

TE: A COMPENSATOR TRANSFER FUNCTION CAN BE 

COMPUTED ONLY IF BOTH A REGULATOR 
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U )-0 0*0 



-p 

i 



AND FILTER ARE SYNTHESIZED 
AND/OR INPUT.} 

SELECT AN OPTION: 1 , 2, 3, OR 4. 



WILL A FEED-FORWARD DISTRIBUTION MATRIX 
("D" - MATRIX} BE INPUT ? 

TYPE "YES" OR "NO". 



THIS OPTION DETERMINES THE CRITERIA FOR DECIDING WHEN A 
MARKCV PARAMETER IS ZEEO-THE MARKOV PARAMETER INDICATES 
THE CRD EE OF THE NUMERATOR POLYNOMIAL OF EACH TRANSFER 
FUNCTION. 



ALL "N" ZEROS CF THIS POLYNOMIAL ARE PRINTED OUT AND 
THIS TEST TELLS HOW MANY EXTRA ROOTS EXIST AT Z = 0. 
LESS THAN 10 . 0 ** {-IE} IS CONSIDERED ZERO. 



THE DEFAULT VALUE OF THIS PARAMETER 
IN OTHER WORDS/ IE = 1.0E-6. 



(IE} IS 6. 



IF YOU DESIRE A DIFFERENT MARKOV CRITERIA/ 

TYPF THE INTEGER VALUE. 

IF YCU DESIRE THE DEFAULT VALUE, TYPE "0" {ZERO} 



POWER SPECTRAL DENSITY {PSD} OPTION 1 : 

OPTION 1 — COMPUTE THE PSD OF THE OUTPUTS AND/DR THE 

CONTROLS OF THE CONTROLLED SYSTEM WHEN FORCED BY 
PROCESS AND MEASUREMENT NOISE. {NOTE: BOTH A 
REGULATOR AND A FILTER MUST BE RESIDENT IN THE 
PROGRAM TO USE THIS OPTION.} 

OPTION 2 — SAME AS OPTION 1 ABOVE BUT ONLY PRINT THE 
RESIDUES OF EACH TRANSFER FUNCTION 
USED IN THE PSD COMPUTATION. 

OPTION 3 — NOT DESIRED. 



SELECT AN OPTION: 1, 2, OR 3. 



THE "F", "G", "H" , "GAM", "A" AND "B" MATRICES 

FROM YOUR PREVIOUS OPTSYS RUN WERE SAVED. 

THE FOLLOWING OPTIONS ARE AVAILABLE: 

1. USE ALL OF THE SAME MATRICES AGAIN. 

2. USE SELECTED MATRICES AGAIN. 

3. INPUT ALL NEW MATRICES. 

ENTER 1, 2, OR 3. 

NOTE: EACH SAVED MATRIX WILL BE REDISPLAYED AT 

THE PROPER INPUT SEQUENCE INTERVAL 
AND YOU WILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 
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FLAG/PARAMETER SETTINGS FOR THIS RUN ARE AS FOLLOWS: 



IOL IQ IR ISS IM 




ITF1 IT 


F2 ITF 3 


IFDFW 


IE 


IDEBUG 


3 0 0 0 


1 


0 


0 0 


0 


0 


0 


ISET IDSTAE IPSD 


IYU INORM 


IREG NS 


NC 


N03 


NG 


0 0 0 


0 


0 


0 4 


0 


0 


0 


ORDER OF SYSTEM = 


4 












NUMBER CF CONTROLS 


= 


0 











NUMEER CF OBSERVATIONS = 0 

NUMBER CF PROCESS NOISE SOURCES = 0 

THE SYSTEM MATRIX {"F"-MATRIX} . . . 



0. 0 


1 .00000 


0.0 


0.0 


0. 0 


-0.41500 


-0.01 11 0 


0.0 


9.80000 


-1 .43000 


-0 . 01980 


0.0 


0. 0 


0.0 


1 . 00000 


0.0 



DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" OR "NO". 



OPEN LOO? DYNAMICS MATRIX 



0 . 0 0.100 0 E +0 1 

0.0 -0.4 15 0D+00 

0. 98C0D + 01 -0.143 0D+0 1 

0.0 0.0 



OPEN LOOP EIGENVALUES. 
0. 0 :-6. 80767D-0 1: 



0.0 0.0 

-0. Ill 0D —01 0.0 

-0.1980D-01 0.0 

0.1000D+01 0.0 



DET (SI-F) . . 

1.22984D-01, 3. 80349D- 0 1: 



OPEN LOOP RIGHT EIGENVECTOR MATRIX T 

0.0 -3 .44 94S3D-02 - 1 . 375658 D- 02 9.725766D-03 

0.0 2.34 830 ID-02 -5. 391 019D-03 -4. 036 1 9 3D-03 

0. 0 5.6 22534D-0 1 1.229836D-01 3. 803490D-01 

1. 000000D+00 -8 .259 115D-01 1.000000D+00 0.0 



OPEN LOOP 

3. 7 38739D+ 01 
-5. 8586C5D+00 
— 4 . 2 2 26 0 8D+ 01 

2. 23 1407D+01 



LEFT EIGENVEC 

9.009 C09D + 0 1 
2 .4 23 391D+0 1 
-7.007 502D+0 1 
-1 .316561D+01 



TOR MATRIX 

-4 . 260481 D- 15 
4. 06 97 40 D- 01 
3.361245D-01 
1.918868D+00 



T-INV. . 

1. 000 0 0 0D + 00 
0 . 0 
0.0 
0.0 



MODAL MEASUREMENT SCALING MATRIX H (BAR) *T. . 

0.0 0.0 0.0 0.0 
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DO YOU WISH TO OBTAIN A TIMS RESPONSE 
OF THE SYSTEM YOU ARE EVALUATING? 

(Y OR N) 

NOTE: YOU MUST BE LOGGED ON AT A DUAL SCREEN 
(TEK 618) TERMINAL TO UTILIZE THIS MODE. 

THE F (SYSTEM) , G (CONTROL). H (OBSERVABLES), GAM (NOISE), 
A (OUTPUT COST) AND B (CONTROL COST) MATRICES WILL BE 
SAVED FOR REENTRY TO THE MAIN OPTSYS PROGRAM. 



*** * * ********* ** ** ****** ********************* **** ******** ** * * * 
IF YOU ARE DISSATISFIED WITH THE RESULTS 
THUS FAR AND WOULD LIKE TO EXIT TO CMS, 

-TYPE * Y * TO EXIT- 

(ANY OTHER INPUT TO CONTINUE) 

******** ****** ************* ****************** **** ******** *** * * 

LOADING CPTCALC. ... 

EXECUTION BEGINS... 



DURING THIS SECTION OF THE PROGRAM YOU WILL: 

- SELECT THE TYPE OF SYSTEM RESPONSE TO PLOT 
(OPEN LOOP, CLOSED LOOP, OR FILTER/REGULATOR) 

- PROVIDE START AND STOP TIME FOR PLOTTING CALCULATIONS 

- SELECT THE TYPE CE DRIVING FUNCTION IS) (STEP OR RAMP) 

- PROVIDE START AND STOP TIMES FOR THE DRIVING FUNCTION (S) 

- PROVIDE DRIVING FUNCTION MAGNITUDE (S) . 

CLEAR THE SCREEN TO CONTINUE 



THE F MATRIX 



0.0 


1.00000 


0.0 


0.0 


0 . 0 


-0. 41500 


-0.01 11 0 


0.0 


S. 80000 


-1 . 43000 


-0.01980 


0.0 


0.0 


0.0 


1.00000 


0.0 



THE G MATRIX 

0.0 

0.0 

0.0 

0.0 



THE C MATRIX 

0 . 0 0.0 0.0 0.0 



THE FOLLOWING PLOTTING OPTIONS ARE AVAILABLE IF THE 
REQUIRED MATRICES WERE CALCULATED IN OPTSYSX: 

1. OPEN LOOP TIME RESPONSE 

XDOT = [F}*X + {G}*UC 

2. CLOSED LOOP TIME RESPONSE 

XDOT = (F-G*C)*X + { G} *UC , U = { C} *X 

OPTIMIZED FILTER CLOSED LOOP SYSTEM RESPONSE. 
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-? 

i 



XDOT =_ (F|*X + 



XHDOT 
FILTER + 



J%?5 



C, 



z 

{K} *(Z 



(F) *XH -f {0} U + 

REGULATOR CLOSED LOOP 
XDOT = {F+G*C} *X + f G} *UC , 
XHDOT = (F}*XH + {G}U + {K}*{Z 



: i^* x 



:*XH} 

SYSTEM RESPONSE 
- H *XH} , 



Z = {H} *X 

U = (C) *X H 



SELECT 1, 2, 3 OR 4. 



7 

6 



AT WHAT TIME DO YOU WANT TO START 
THE TIME RESPONSE CALCULATIONS? 

INPUT START TIME IN SECONDS. (NORMALLY 0.0) 
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AT WHAT TIME DO YOU WANT TO STOP 
THE TIME RESPONSE CALCULATIONS? 

INPUT STOP TIME IN SECONDS. 



500 



THIS PROGRAM DIVIDES THE TIME INTERVAL YOU SAVE 
JUST SPECIFIED INTO UP TO 500 SMALL INTERVALS FOR 
THE INTEGRATION AND PLOTTING ROUTINES. IN ORDER 
TO SAVE COMPUTER TIME, THE NUMBER OF POINTS CAN BE 
CAN BE REDUCED WITH SOME LOSS IN CURVE FIDELITY. 

HOW MANY POINTS DO YOU WANT TO CALCULATE? 



n 



DOES THE SYSTEM UTILIZE A DRIVING FUNCTION (CONTROL INPUT)? 

(Y) ES OR (N) O 



DOES THE SYSTEM START WITH ALL INITIAL CONDITIONS = 0.0 ? 

(Y) ES OR (N) O? 



n 



*? 

6.02 



“7 

6 

■7 

6 

6 



WHAT 


IS 


THE INITIAL 


CONDITION 


FOR 


X ( 


11 


? 


WHAT 


IS 


THE INITIAL 


CONDITION 


FOR 


X ( 


2) 


7 


WHAT 


IS 


THE INITIAL 


CONDITION 


FOR 


X ( 


3) 


7 


WHAT 


IS 


THE INITIAL 


CONDITION 


FOR 


X ( 


4) 


7 


THIS 

MAKE 

1 . 


IS YOUR LAST OPPORTUNITY 
CHANGES IN THE FOLLOWING 

SELECT ANOTHER TYPE OF S 


TO 

AREAS. 
YSTEM TO 


PLOT 
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(OPEN/ CLOSED/ FILTER OR FILTER/3 EGULATCR) 

2. START AND STOP TIMES 

3. DRIVING FUNCTIONS 

4. INITIAL CONDITIONS 

5. CONTINUE 

SELECT A NUMEER BETWEEN 1 AND 5. 



THE FOLLOWING INFORMATION IS PROVIDED ONLY 
FOR AN INDICATION OF PROPER PROGRAM OPERATION. 

AIL CONTROLS, STATES AND STATE ESTIMATES CAN BE PLOTTED. 



TIME 


0(1) 


0.0 


0 . 0 


0.50 


0 . 0 


1. 0 0 


0 . 0 


1.50 


0 . 0 


2. 00 


0 . 0 


2. 50 


0 . 0 


3. 00 


0 . 0 


3.50 


0 . 0 


4.00 


0 . 0 


4. 50 


0 . 0 


5. 00 


0 . 0 


5.50 


0 . 0 


6.00 


0 . 0 


6.5 0 


0.0 


7. 00 


0 . 0 


7.50 


0 . 0 


8. 00 


0 . 0 


8.50 


0 . 0 


9.00 


0 . 0 


9.50 


0 . 0 


10.00 


0 . 0 


10. 50 


0 . 0 


11. 00 


0 . 0 


11.50 


0 . 0 


12. 00 


0 . 0 


12. 50 


0 . 0 


13.00 


0 . 0 


13.50 


0 . 0 


14. 00 


0 . 0 


14.50 


0 . 0 


15.00 


0 . 0 


15. 50 


0 . 0 


16. 00 


0 . 0 


16.50 


0 . 0 


17. 00 


0 . 0 


17.50 


0 . 0 


18. 00 


0 . 0 


18.50 


0 . 0 


19. 00 


0 . 0 


1 9.50 


0 . 0 


20.00 


0 . 0 


20.50 


0 . 0 


21.00 


0 . 0 


21. 50 


0 . 0 


22. 00 


0.0 


22. 50 


0 . 0 


23. 00 


0 . 0 


23.50 


0 . 0 



X(1) X(2) 

0.20000 00D-01 0.0 0. 

0.19957 04D-0 1-0. 2532513 D- 03 0. 
0.1967379D-01-0.9447647D-03 0. 
0. 1895459D-01-0. 1 982623D-02 0. 
0. 17647 46D-01-0. 328281 9 D- 02 0. 
0. 15641 62D-0 1-0. 4763866D-02 0. 
0. 12867 15D-01-0. 634347 ID-02 3. 
0.9295838D-02-0.7936863D-02 0. 
0. 4942574D-02-0.9456467D-02 0. 
-0. 13339 96D-03-0. 1081 268D-01 0. 
■0.58 27 877D-02-0.1 191 55 9D- 01 0. 
•0.11 99 21 2D -01-0. 1267744 D- 01 0. 
•0. 1843467D-01-0.1301571D-01 0. 
- 0.24 924 87D-0 1-0. 1 285666 D- 01 0. 
■0.3119815D-01-0. 121391 0D -01 0. 
•0.36963 32D-01-0. 1081832D-01 0. 
•0.41911 73D-0 1-0. 886 990 2D- 02-3 . 
•0.45728 36D-0 1-0 . 6 2 9323 8D- 02- 0 . 
-0.4810447D-01-0.3 1 1 4 54 6D- 02-0 . 
•0.48751 63D-0 1 0 . 6 10 803 9D- 03- 0 . 
•0.47416 70D-0 1 0 . 4 7 96 89 0D-02-0 . 

• 0. 43 897 25D-0 1 0. 9 326 9 1 7D- 02-0 . 
•0.38056 76D-01 0. 1 4054 16D-0 1-3 . 
•0.29839 14D-01 0. 1 880428D-01-0 . 
•0. 19281 59D-01 0 . 2 3 379 1 3D- 0 1-0 . 

• 0. 65 252 73D-0 2 0 . 2 756 20 ID- 01 -3 . 
0.8177 0 35D-0 2 0 . 3 1 1 2 44 5D- 01- 0 . 
0. 24457 06D-0 1 0 . 3 3 8343 8D- 01-0 . 
0.41831 51D-01 0.3546543D-01-0 . 
0.5970 7 49D-0 1 0.3580735D-01-0. 
0.77393 26D-0 1 0 . 3 46769 3D- 01 -3 . 
0.94114 63D-0 1 0.3 192943D-01-0 . 
0. 10903 70D + 00 0 .2 746 96 1 D-0 1 -3 . 
0.12129 26D + 00 0 . 2 1 26 23 4D-0 1 0. 
0. 13001 34D + 00 0. 1 334194D-01 0. 
0. 1343671D + 00 0 . 3 8 1 976 4D- 02 0. 
0.13359 8 0D + 00-0. 7 1 10370 D- 02 0. 
0. 1270681D+00-0. 1916965D-01 0. 
0. 1143002D + 00-0.3 199492D-01 3. 
0.9501896D-01-0.4514306D-01 0. 
0.6918881D-01-0. 580993 0D-01 0. 
0.37047 08D-0 1-0.702895 8D- 01 0, 

•0.8704588D-03-0. 8109723D- 01 0. 
•0.4371323D-01-0.8988342D-01 0. 
•0.90310 68D-01-0. 9601 136D-01 3. 
■0. 1391800D+00-0.9887356D-01 0. 
•0. 1885480D+00-0.9792139D-01 0. 
•0.2363885D+00-0.9269602D-01 0. 



X (3) 



9752469D-0 1 
1937236D+00 
28722 00D+00 
3757431 D+0 0 
4563456 D+00 
52561 08D+00 
57 98568D+0 0 
6 153434 D+00 
62 84 849D+ 0 0 
6 1 606 55D+0 0 
5754551D+00 
50481 90D+00 
40331 60D+00 
2712 751 D+00 

I 103423D+00 
7641071 D- 01 
2844264D+00 
5076571 D+00 
7386238D+00 
96854 50D+00 

I I 87536D+0 1 
13848 85 D+ 01 
15493 96D+0 1 
1669798D+0 1 
1 7352 16D+0 1 
1735684 D+01 
16 62 679D+ 0 1 
1509 680D+0 1 
12 72 699D+0 1 
9507871 D+00 
54 646 64D+0 0 
6608227 D-0 1 
4799671 D+00 
1077141 D+01 
1 706 829D+0 1 
23466 10 D+0 1 
29706 98D+01 
35 505 99D+0 1 
4055 S63D+0 1 
4455624 D+0 1 
4718817 D+01 
4816533 D+01 
4722976 D+01 
4417079 D+0 1 
38840 14D+01 
31 16637 D+01 
21 16795 D+01 
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2a. oo o.o 

24.50 0.0 

25. 00 0.0 



-0.2804747D+00-0.8285960D-01 0. 8964261 D+ 00 
-0.31 844 75D + 00-0. 6 82254 5D-01- 0 . 5216259D+00 
- 0.34 789 8 1D + 0 0-0. 4 87 857 1D-0 1-0 . 2 1 0 3 1 53 D+0 1 



*********** ********* ************************* **** ************* 



IF YOU ARE DISSATISFIED WITH THE RESULTS 
THUS FAR AND WOULD LIKE TO EXIT TO CHS, 

-TYPE ’ Y ' TO EXIT- 



(ANY OTHER INPUT TO CONTINUE) 

************************************* ************************* 



E 

C 

•p 



12 0 ) 
121 ) 
122 ) 



R/O 

R/O 

R/O 



... Your Fortran program is now being loaded ... 
... execution will soon follow ... 
EXECUTION BEGINS... 



THIS PCRTICN OF THE PROGRAM PLOTS: 

_ fup FTATF^ 

- EXTERNAL "c 6 NT ROL INPUTS, 

- FEEDBACK CONTROL INPUTS, 

- STATE ESTIMATES AND 

- RECONSTRUCTION ERRORS 

FROM THE DATA THAT YOU JUST CALCULATED. 

THE CAPABILITY IS ALSO AVAILABLE TO REVIEW ANY 
GRAPHS THAT YOU HAD PREVIOUSLY SAVED AS DATA 
FILES ON YOUR DISK. 

CLEAR THE SCREEN TO CONTINUE. 



THE FOLLOWING OPTIONS ARE AVAILABLE: 

1. PL CT THE DATA YOU JUST CALCULATED. 

2. PLOT A CURVE THAT YOU PREVIOUSLY SAVED. 

ENTER 1 OR 2 

•? 

i 



YOU MAY PLOT UP TO 4 SYSTEM VARIABLES VS TIME. 
HOW MANY VARIABLES DC YOU WISH TO PLOT? 



WHICH TYPE OF V ARIAEIE DO YOU WISH TO PLOT AS CURVE NUMBER 1? 

1. STATE VARIABLE (IE., XI. X2, ETC) 

2. EEEDEACK CONTROL (IE., U = -C*X) 

3. CCNTROL INPUT (IE., U1, U2, ETC.) 

4. STATE ESTIMATE (OBSERVER) (IE., XHAT1, XHAT2, ETC.) 

5. STATE RECONSTRUCTION ERROR (IE., X1-XHAT1, 

X2-XHAT2, ETC) 

ENTER 1,2, 3,4 OR 5 
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1 



WHAT IS THE SUBSCRIPT OF THE STATE VARIABLE THAT 
YCU WANT TO PLOT AS THE NUMBER 1 CURVE VS TIME? 

•7 

i 



WHAT IS THE CURVE LAEEL FOR THIS VARIABLE? 



NOTE: 



state yl 



1. 40 CHARACTERS MAX LENGTH 

2. GREEK SYMBOLS WILL BE PRINTED FOE ANY LETTERS 
ENCLOSED IN PARENTHESES. 



IE. 



A) => 
E| => 
F) => 
Q => 



ALPHA 

BETA 

PHI 

THETA 



WHICH TYPE OF VARIAEIE DO YOU WISH TO PLOT AS CURVE NUMBER 2? 



-? 

i 



1 , 

2 , 
■3 

i 

4, 

5, 



STATE VARIABLE (IE., XI, X2, ETC) 
FEEDEACK CONTFCL (IE.^ U = -C*X) 



CONTROL INPUT (IE., Ul, U2, ETC.) 

STATE ESTIMATE (OBSERVER) (IE., X HAT 1 , XHAT2, ETC.) 
STATE RECONSTRUCTION ERROR (IE., Xl-XHATl, 



X2-XHAT2, ETC) 

ENTER 1,2, 3,4 OR 5 



o 

2 



WHAT IS THE SUBSCRIPT OF THE STATE VARIABLE THAI 
YOU WANT TO PLOT AS THE NUMBER 2 CURVE VS TIME? 



WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 



NOTE: 



state y2 



1. 40 CHARACTERS MAX LENGTH 

2. GREEK SYMBOLS WILL BE PRINTED 
ENCLOSED IN PARENTHESES. 



FOR ANY LETTERS 



IE. 



A) => ALPHA 

E) => BETA 

F) => PHI 

Q) => THETA 



WHICH TYPE OF VARIABLE DO YOU WISH TO PLOT AS CURVE NUMBER 3? 

1. STATE VARIABLE (IE., XI, X2, ETC) 

2. FEEDEACK CONTFOL (IE.. 0 = -C*X) 

3. CONTROL INPUT (IE., Ul, U2, ETC.) 

4. STATE ESTIMATE (OBSERVER) (IE., XHAT1, XHAT2, ETC.) 

5. STATE RECONSTRUCTION ERROR (IE., XI -XH ATI , 

X2-XHAT2, ETC) 

ENTER 1,2, 3,4 OR 5 
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•5 

3 



W HAT IS THE SUBSC 
YOU WANT TC PLOT 



RIPT or THE STATE VARIABLE THAT 
AS THE NUMBER 3 CURVE VS TIME? 



WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 



NOTE; 



state y3 



1. 40 CHARACTERS MAX LENGTH 

2. GREEK SYMBOLS WILL BE PRINTED 
ENCLOSED IN PARENTHESES. 



FOR ANY LETTERS 



I E. 




ALPHA 

BETA 

PHI 

THETA 



WHICH TYPE OF VARIABLE DO YOU WISH TO PLOT AS CJ RVE NUMBER 4? 



1 . 
2 . 

3. 

4. 

5. 



-> 

i 



STATE VARIABLE (IE.. XI, X2, ETC) 

EEEDEACK CONTROL (IE., 0 = -C*X) 

CONTROL INPUT (IE., Ul, 02, ETC.) 

STATE ESTIMATE (OBSERVER) (IE., XHAT1, XH&T2, 
STATE RECONSTRUCTION ERROR (IE., X1-XHAT1, 
X2-XHAT2, ETC) 



ETC. ) 



ENTER 1, 2, 3,4 OR 5 



4 



WHAT IS THE SUBSC 
YOU WANT TC PLOT 



RIPT OF THE STATE VARIABLE THAT 
AS THE NUMBER 4 CURVE VS TIME? 



WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 



NCTE: 



state y4 



1. 40 CHARACTERS MAX LENGTH 

2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 



IE. 




ALPHA 
BETA 
PHI 
T HETA 



3 



YCU MAY USE UP TO 3 HEADINGS. 

HOW MANY HEADINGS DO YOU DESIRE ON THIS GRAPH? 
0, 1, 2 OR 3 



WHAT IS THE DESIRED HEADING NUMBER 1? 



NOTE: 1. 

2 . 



open loop system 



40 CHARACTERS MAX LENGTH 

GREEK SYMBOLS WILL BE PRINTED FOR 

LETTERS ENCLOSED IN PARENTHESES. 



ANY 



IE. 



(A) 

3 

F ; 

(Q 



=> 

=> 

=> 

=> 



ALPHA 

BETA 

PHI 

THETA 
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WHAT IS THE DESIRED HEADING NUMBER 2? 



example 2 



NOTE: 1. 40 CHARACTERS MAX LENGTH 

2. GREEK SYMBOLS WILL BE PRINTED FOR 
LETTERS ENCLOSED IN PARENTHESES. 
IE. (A) => ALPHA 

’B) => BETA 

F) => PHI 
Q) => THETA 



ANY 



WHAT IS THE DESIRED HEADING NUMBER 3? 



NOTE: 1. 40 CHARACTERS MAX LENGTH 

2. GREEK SYMBOLS WILL BE PRINTED FOR 
LETTERS ENCLOSED IN PARENTHESES. 
IE. (A) => ALPHA 

(B) => BETA 

(F) => PHI 

(Q) => THETA 

modern control theory 

» USING A PRE-ALLOCATED DATASET FOR UNIT FT17F001. 



ANY 



TEE FOLLOWING OPTIONS ARE AVAILABLE. 

1. BEGIN NEW GRAPH OF OTHER CONTROLS. STATES, 

2. REPLCT PREVIOUSLY SAVED GRAPH DATA. 

3. EDIT THE CURRENT GRAPH. 

4. PLOT REVISED GRAPH ON THE TEK618. 

5. QUIT AND/OR MAKE METAFILE OF THE CURVES. 
PREVIOUSLY SAVED. 



OR ESTIMATES. 



SELECT 



A NUMBER EETWEEN 1 AND 5. 



THE GRAPH EDIT MENU 

1. CHANGE VARIABLES OR ADD A CURVE ON THE CURRENT PLOT. 

2. DELETE CURVE FROM CURRENT PLOT. 

3. EDIT CURVE TITIE(S). 

4. EDIT PAGE HEADING (S) . 

5. CHANGE THE Y-AXIS SCALE. 

6. CHANGE THE TIME AXIS SCALE. 

7. CHANGE PLOT SIZE. {DEFAULT IS 8 . 5 X 6.0) 

3. CHANGE THE LETTERING HEIGHT. 

9. CHANGE POSITION OF THE LEGEND. 

10. EDITING COMPLETE. 

SELECT A NUMBER EETWEEN 1 AND 10. 



-4 



HOW MANY 
(LEFT OR 
MOVE THE 

NOTE: 1. 

2 . 

3. 



INCHES IN THE X DIRECTION 

RIGHT), DO YOU WANT TO MOVE 

LEGEND BOX FROM ITS PRESENT POSITION 

DEFAULT PLOT SIZE IS 3.5 X 6.0 
LEFT IS NEGATIVE 
RIGHT IS POSITIVE 



HOW MANY INCHES IN THE Y DIRECTION 

(UP OR DOWN). DO YOU WANT TO MOVE 

MOVE THE LEGEND BOX FROM ITS PRESENT POSITION 
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NOTE: 1. DEFAULT PAGE SIZE IS 8.5 X 6.0 

2. DCWN IS NEGATIVE 

3. OP IS PO SITIVE 



THE GRAPH EDIT MEND 



1. CHANGE VARIABLES OR ADD A CURVE 

2. DELETE CURVE FROM CURRENT PLOT. 

3. EDIT CURVE TITIE(S). 

4. EDIT PAGE HEADING (5). 

5. CHANGE THE Y-AXIS SCALE. 

6. CHANGE THE TIME AXIS SCALE. 

7. CHANGE PLOT SIZE. (DEFAULT IS 8.5 X 6.0) 

8. CHANGE THE LETTERING HEIGHT. 

9. CHANGE POSITION OF THE LEGEND. 

10. EDITING COMPLETE. 



ON THE CURRENT PLOT. 



SELECT A NUMBER BETWEEN 1 AND 10. 

io 



THE FOLLOWING OPTIONS ARE AVAILABLE. 

1. BEGIN NEW GRAPH CF OTHER CONTROLS, STATES, OR ESTIMATES. 

2. REPLCT PREVIOUSLY SAVED GRAPH DATA. 

3. EDIT THE CURRENT GRAPH. 

4. PLOT REVISED GRAPH ON THE TEK618. 

5. QUIT AND/OR MAKE METAFILE OF THE CURVES. 

PREVIOUSLY SAVED. 

SELECT A NUMBER BETWEEN 1 AND 5. 



THE FOLLOWING OPTIONS ARE AVAILABLE. 

1. BEGIN NEW GRAPH CE OTHER CONTROLS, STATES, OR ESTIMATES. 

2. REPLCT PREVIOUSLY SAVED GRAPH DATA. 

3. EDIT THE CURRENT GRAPH. 

4. PLOT REVISED GRAPH ON THE TEK618. 

5. QUIT AND/OR MAKE METAFILE OF THE CURVES. 

PREVIOUSLY SAVED. 

SELECT A NUMBER BETWEEN 1 AND 5. 



DO YOU WANT TO SAVE THE CURRENT GRAPH DATA TO 
EE USED LATER TO GENERATE A METAFILE? 

Y OR N 

NOTE: A METAFILE IS REQUIRED FOR SMOOTH VERS ATEC PLOTS. 
THERE WILL BE AN OPPORTUNITY TO GENERATE A METAFILE 
JUST BEFORE EXITING THIS PROGRAM. 



WHAT FILE NAME 



openlcop 



DC YOU WANT THE CURVE DATA STORED UNDER? 
(8 CHARACTERS MAX) 



THE CURVE DATA IS BEING FILED UNDER OPE NLOOP DATA 
END CF DISSPLA 9.0 -- 26506 VECTORS GENERATED IN 2 PLOT FRAMES 
PROPRIETARY SOFTWARE PRODUCT OF ISSCO, SAN DIEGO, CA. 

3442 VIRTUAL STORAGE REFERENCES; 6 READS; 0 WRITES. 
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THE FOLLOWING OPTIONS ARE AVAILABLE: 



1. MAKE METAFILE OF PREVIOUSLY SAVED DURVE. 

2. QUIT. 

ENTER 1 OR 2 

-7 

i 



WHAT FILE NAME IS THE DATA STORED UNDER? 

openloop 

THE CURVE DATA IS BEING LOADED FROM FILE OPENLOOP DATA 
» USING A PRE-ALLOCATED DATASET FOR UNIT FT18F001. 



THE FOLLOWING OPTIONS ARE AVAILABLE: 

1. MAKE METAFILE OF PREVIOUSLY SAVED GURVE. 

2. QUIT. 

ENTER 1 OR 2 

7 

2 

END OF DISSPLA 9.0 -- 13197 VECTORS GENERATED IN 1 PLOT FRAMES 
PROPRIETARY SOFTWARE PRODUCT OF ISSCO, SAN DIEGO, CA. 

1817 VIRTUAL STORAGE REFERENCES; 5 READS; 0 WRITES. 

DASD 121 DETACHED 
DASD 122 DETACHED 
DASD 120 DETACHED 



*** * * ** * * *** * ** * **££****:£*******:{£* *** ********* * ** * £ 

DO YOU WANT A VRSTEC PLOTTER SMOOTH COPY OF THE 
TEE DISSPLA METAFILE THAT YOU JUST CREATED? 

(Y OR N) 

I*********** *** * * ** 4 I******************************* **** *** ** * * * 



E (120) E/O 

DASD 001 LINKED R/O ; R/W BY MVS 
Z (00 1) E/C - OS 
DASD 001 DETACHED 
CREATING NEW FILE: 

CREATING NEK FILE: 

PUN FILE 6680 TO MVS COPY 001 NOHOLD 

DASD 120 DETACHED 



**************** ** * * ****** ***************** ** *** ******* ** *** * * 
YO'JE GRAPH (S) CAN BE PICKED UP AT THE COMPUTER CENTER. 

TEE GRAPH (S) WILL BE ADDRESSED TO "POP (USER ID) ". 
**************** *** ******** ********************** ******** *** * * 



*********** ***** ** ************************************* ** *** * * 
DO YOU WANT TO 

1. RUN OPTSYSX AGAIN 

2. RUN THE PLOT PROGRAM USING THE SAME MATRICES? 
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(TC PLOT ANOTHER TYPE OF SYSTEM (OPEN/CLOSED)) 

3. QUIT 

ENTER 1 , 2 OR 3 

********************* ****** ********************** ******** *** * * 
3 

*********** *** ************************************************ 

HAVE A GOOD DAY! ! 

*************** *** ********* ****************** ** ** ************* 
R; T=11. 05/17.26 20:C8:15 
record off 

END RECORDING OF TERMINAL SESSION 

The graphical output generated by this example follows 
as figure 3.1. 
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OPEN LOOP SYSTEM 
EXAMPLE 2 

MODERN CONTROL THEORY 



© 




TA 31V1S 



I 1 I I I l I 

050*0 520 0 000*0 520*0- 050*0- 5Z.00- 001*0- 

ZX 31V1S 

f 1 1 1 T 1 1 1 1 



0*5 


o'* oc o z ot oo n- o'z- o - e- 

CA 31V AS 


0*02 


! ] i 1 1 1 

O'ST 00t 09 0 0 0 9- 001- 



U 

W 

00 

I 

u 

2 



f A 31V IS 
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Figure 3.1 Open-loop Time Response 



E. CICSED-LOOP SYSTEM TIME HESPONSE 

The following closed-loop system example was taken from 
[Ref. 6, pp 5.8 - 5.19]. 

The full terminal session is recorded below, with user 
input at the left margin in lower case l etter s or n u mbers 
below each . 



record on 

EEGIN RECORDING OF TERMINAL SESSION 

R; T = C. 01/0.02 20: 1 9 :44 

optsys 



IKE OPTSYS EXEC CONTROLS A TRIO OF PROGRAMS: 



1. OPTS YS X FORTRAN 

2. OPTCALC FORTRAN 

3. OPTPLOT FORTRAN 



SYSTEM ANALYSIS) 

CALCULATE TIME RESPONSE) 
DISSPLA PLOTTING ROUTINE) 



EACH PROGRAM PASSES INFORMATION TO THE NEXT PROGRAM 
THROUGH A DATA FILE WRITTEN TO THE USERS DISK. IN THIS 
CASE, THESE FILES ARE "OPTMAT DATA" AND "OPTPLOT DATA". 
THE SIZE OF THESE FILES VARY WITH THE SYSTEM ORDER, AND 
CAN USE ABOUT 20? OF THE USERS DISK SPACE. THEREFORE 
ENSURE THAT SUFFICIENT DISK SPACE IS AVAILABLE. 



- TYPE "E" TO EXIT, ANY OTHER ENTRY TO CONTINUE - 



******** *************%********************* ** *** * * £ * * * * ** *** * * 

YOU HAVE A DATA FILE NAMED ’OPTMAT DATA’ ON YOUR A 
DISK THAT WAS PREVIOUSLY GENERATED BY THE OPTSYS 
PROGRAM AND CONTAINS THE F, G, H, GAMMA, A AND B 
MATRICES FROM THAT RUN. 

IF YOU WOULD LIKE TO USE THESE SAME MATRICES FOR 
THIS RUN, THE OPTSYS PROGRAM WILL READ IN THE 
DESIRED DATA AT TBE APPROPRIATE TIME, 

IF YOU TYPE (Y) ES. 

ANY OTHER INPUT WILL RESULT IN THAT FILE BEING ERASED! 

** * * * ** ** ** *** ***** ******** ****************************** *** * * 



*** * * ****** *** ** ** * ******** ****************** **** ************* 

DC YCU WANT THE NUMERICAL OUTPUT FROM OPTSYSX TO GO 
TO YOUR TERMINAL S(CREEN) OR TO A D (ISK) FILE? 

{S OR D) 

************************************************************** 



************************************************************** 
OUTPUT WILL COME TO YOUR TERMINAL SCREEN. 
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************************************************************** 



LOADING CPTSYS. . .. 

EXECUTION BEGINS... 

OPTSYSX IS A COMPLETELY INTERACTIVE OPTIMAL SYSTEMS CONTROL 
PROGRAM. IT FILL SOLVE NUMEROUS CONTROL PROBLEMS ON THE 
FOLLOWING TYPES OF SYSTEMS CONTROL EQUATIONS: 

X DOT = {F}*X + (G}*U + £GAM}*(W+W0) 

MEASUREMENT EQUATION— 

Z = £H}*X + {D} *U + V 

REGULATOR PERFORMANCE INDEX — 

J = 1/2 * INTEGRAL (Y *{A}*Y + U * {B} *0) DT 

STATE FEEDBACK GAIN DEFINITION — 

U = - {C} *X 

DO YOU WISH TO CONTINUE? TYPE "YES" OR "NO". 
— DATA ENTRY — 

AITHOUGH OPTSYSX IS SPECIFICALLY DESIGNED TD BEAD 
ALL MATRIX DATA INTERACTIVELY. SEVERAL ALTERNATE 
METHODS ARE AVAILABLE TO USERS: 

METHOD 1— THE "F" «G» AND "GAMMA" MATRICES 
MAY BE READ FROM SEPARATE DATA FILES. 

METHOD 2 — THE "F","G», AND "GAMMA" MATRICES MAY BE 
EXPLICITLY DEFINED WITHIN SUBROUTINE "SETUP". 

{NOTE: IN EITHER CASE, THE USER SHOULD 03TAIN A COPY 
OF THE PROGRAM LISTING AND EXAMINE 
THE EXAMPLES CONTAINED IN S/R "SETUP".} 

DO YOU WISH TO CONTINUE? TYPE "YES" OR "NO". 

y 



DO YOU WISH TO INPUT THE "F" , "G", AND "GAMMA" 
MATRICES FROM SUBROUTINE "SETUP" IAW THE 
METHOD DESCRIBED ON THE PREVIOUS SCREEN? 

TYPE "YES" OR "NO". 

GENERAL OPTSYSX OPTIONS: 



OPTION 

OPTION 

OPTION 



1 

2 

3 



CPTION 4 — 



SYSTEM ANALYSIS WITHOUT 

CPEN-L OOP EIGENS YSTEM CALCULATIONS. 

SYSTEM ANALYSIS WITH OPEN-LOOP 
EIGENS YSTEM CALCULATIONS. 

CPEN-L OOP EIGENS YSTEM FOUND 
AND PROGRAM TERMINATES. 

{"F"-MATRIX ENTRY FOLLOWS IMMEDIATELY.} 

MODAL DISTRIBUTION MATRICES COMPUTED 
WITHOUT FILTER OR REGULATOR SYNTHESIS 
OR STEADY-STATE ANALYSIS. 



SELECT AN OPTION: 1,2,3, OR 4. 
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1 



n 



i 



DO YOU DESIRE RMS VALUES OF STATE AND CONTROL? 
TYPE "YES" OR "NO". 



CPTSYSX LQR/CLA 

OPTION 1 — OPTIMA 1 FILTER 
SYNTHESIS WITH 
MATRIX INPUT. 

OPTION 2 — OPTIMAL FILTER 
SYNTHESIS WITH 
MATRIX INPUT. 

OPTION 3 — OPTIMAL FILTER 
SYNTHESIS WITH 
MATRIX INPUT. 

OPTION 4 — OPTIMAL FILTER 
SYNTHESIS WITH 
MATRIX INPUT. 

SELECT AN OPTION: 1, 2, 3 



SSICAL OPTIONS : 

AND/OR REGULATOR 
NO EXTERNAL " C» OR "K" 



AND/OR REGULATOR 
EXTERNAL "C" 



AND/OR REGULATOR 
EXTERNAL "K" 



AND/OR REGULATOR 
EXTERNAL "C» AND "K" 



, OR 4. 



DO YOU WISH TO DETERMINE THE STEADY-STATE RESPONSE 
FOR A CONSTANT DISTURBANCE? 

TYPE "YES" CR "NO", 
n 

DO YOU WISH TO DETERMINE THE MODAL DISTRIBUTION 
AND GAIN MATRICES? 

TYPE "YES" CR "NO", 
n 

CPE N-LCOP TRANSFER FUNCTION OPTIONS: 

OPTION 1 — NO OPEN-LOOP TRANSFER FUNCTIONS COMPUTED. 
CPTION 2 — POLES, RESIDUES, AND ZEROS COMPUTED. 
OPTION 3 ▼— ONLY POLES AND ZEROS COMPUTED. 

CPTION 4 -- ONLY POLES AND RESIDUES COMPUTED. 

SELECT AN OPTION: 1, 2, 3, OR 4. 

1 

NOISE TRANSFER FUNCTION OPTIONS: 

CPTION 1 — NO NOISE TRANSFER FUNCTIONS COMPUTED. 
CPTION 2 — POLES, RESIDUES, AND ZEROS COMPUTED. 
CPTION 3 — ONLY POLES AND ZEROS COMPUTED. 

CPTION 4 -- CNLY POLES AND RESIDUES COMPUTED. 

SEIECT AN OPTION: 1, 2, 3, OR 4. 
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0*0 



COMPENSATOR TRANSFER FUNCTION OPTIONS: 



i 



OPTION 1 — NO COMP. TRANSFER FUNCTIONS COMPUTED. 

CETION 2 — POLES, RESIDUES, AND ZEROS OOMPUTED. 

CPTION 3 — CNLY POLES AND ZEROS COMPUTED. 

OPTION 4 — CNLY POLES AND RESIDUES COMPUTED. 

{NOTE: A COMPENSATOR TRANSFER FUNCTION CAR BE 
COMPUTED ONLY IF BOTH A REGULATOR 
AND FILTER ARE SYNTHESIZED 
AND/OR INPUT.} 

SELECT AN OPTION: 1, 2, 3, OR 4. 

i 



WILL A FEED-FORWARD DISTRIBUTION MATRIX 
{"D" - MATRIX} EE INPUT ? 

TYPE "YES" OR "NO”, 

n 

THIS OPTION DETERMINES THE CRITERIA FOR DECIDING WHEN A 
MARKOV PARAMETER IS ZERO-THS MARKOV PARAMETER INDICATES 
THE CRDER OF THE NUMERATOR POLYNOMIAL OF EACH TRANSFER 
FUNCTION. 

ALL "N" ZEROS CF THIS POLYNOMIAL ARE PRINTED OUT AND 
THIS TEST TELLS HOW MANY EXTRA ROOTS EXIST AT Z = 0 . 
LESS THAN 10.0** {”12} IS CONSIDERED ZERO. 

THE DEFAULT VALUE OF THIS PARAMETER {IE} IS 6. 

IN CTREE WORDS, IE = 1.0E-6. 

IF YOU DESIRE A DIFFERENT MARKOV CRITERIA, 

TYPE THE INTEGER VALUE. 

IF YOU DESIRE THE DEFAULT VALUE, TYPE "0" {ZERO} 



DO YOU DESIRE TO SYNTHESIZE A STABLE FILTER {OR REGULATOR} BY 
DESTABILIZING TEE ORIGINAL SYSTEM? 

{NOTE: WORKS FOR FILTER OR REGULATOR BUT NOT FOE EOT H 
IN THE SAME RUN.} 

TYPE "YES" CR "NO", 

n 

DC YOU DESIRE TO PRINT THE EULER -LAGR ANG E El GENSYSTEM 
PRICE TO DECCMPOSIT ION {FOR CHECKING THE PROGRAM}? 

TYPE "YES" OR "NO", 

n 



POWER SPECTRAL DENSITY {PSD} OPTION 1 : 

OPTION 1 -- COMPUTE THE PSD OF THE OUTPUTS AND/OR THE 

CONTROLS OF THE CONTROLLED SYSTEM WHEN FORCED BY 
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(xJ*0 



OPTION 2 



PROCESS AND MEASUREMENT NOISE. {NOTE: BO 
REGULATOR AND A FILTER MUST EE RESIDENT I 
PROGRAM TO USE THIS OPTION.} 

-- SAME AS OPTION 1 ABOVE BUT ONLY PRINT THE 
RESIDUES OF EACH TRANSFER FUNCTION 
USED IN THE PSD COMPUTATION. 

OPTION 3 — NOT DESIRED. 



SELECT AN OPTION: 1, 2, OR 3. 



DC YOU DESIRE REGULATOR SYNTHESIS ONLY? 
TYPE ''YES" OR "NO". 



y 



THE "F", '' G", "H" , "GAM", "A" AND " E" MATRICES 
FROM YOUR PREVIOUS OPTSYS RUN WERE SAVED. 

THE FOLLOWING OPTIONS ARE AVAILABLE: 

1. USE ALL OF THE SAME MATRICES AGAIN. 

2. USE SELECTED MATRICES AGAIN. 

3. INPUT ALL NEW MATRICES. 

ENTER 1, 2, OR 3. 

NOTE: EACH SAVED MATRIX WILL BE REDISPLAYED AT 

THE PROPER INPUT SEQUENCE INTERVAL 
AND YOU WILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 

2 



DO YOU WISH TO SAVE THE "F''-MATRIX FROM THE LAST 
RUN TO BE USED IN THIS RUN? 



NOTE: THE MATRIX WILL BE REDISPLAYED AT 
THE PRCPER INPUT SEQUENCE INTERVAL 
AND YOU WILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 

TYPE "YES" OR "NO". 



DO YOU WISH TO SAVE THE "A''-MATRIX FROM THE LAST 
RUN TO BE USED IN THIS RUN? 

NOTE: THE MATRIX WILL BE REDISPLAYED AT 
THE PRCPER INPUT SEQUENCE INTERVAL 
AND YOU WILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 

TYPE "YES" OR "NO". 



DO YOU WISH TO SAVE THE "E"-MATRIX FROM THE LAST 
RUN TO BE USED IN THIS RUN? 
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NCTE; SHE MATRIX WILL BE REDISPLAYED AT 
THE PROPER INPUT SEQUENCE INTERVAL 
AND YOU KILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 



n 



TYPE "YES" OR "NO". 



ENTER THE # OF CONTROLS fNC} OF THE CONTROL SYSTEM MODEL 
{" G"-MATRIX} . 

i 



ENTER THE # OF MEASUREMENTS OR OBSERVATIONS {NO} OF THE 
{"H "-M ATR IX} . 



ENTER THE # OF PROCESS NOISE SOURCES {NG} OF THE 
{"G AM M A"-MATRIX} . 

6 

FIAG/PAEAME TER SETTINGS FOR THIS RUN ARE AS FOLLOWS: 
IOL IQ IR ISS IM ITF1 ITF2 ITF3 IFDFW IE IDEBUG 
000000 0 0 0 00 
ISET IDSTAB IPSD IYU INORM IREG NS NC NOB NG 
0 0 00 0 14140 



ORDER OF SYSTEM = 4 

NUMBER OF CONTROLS = 1 

NUMBER OF OBSERVATIONS = 4 

NUMEER OF PROCESS NOISE SOURCES = 0 



THE SYSTEM MATRIX {"F"-MATRIX} . . . 



0 . 0 


1 . 00000 


0.0 


0.0 


0 . 0 


-0.41500 


-0.01110 


0.0 


9.80000 


-1. 43000 


-0 . 01 980 


0.0 


0.0 


0.0 


1.00000 


0.0 



DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" OR "NO". 



OPEN LOO? DYNAMICS MATRIX F 



0.0 


0.100 0D+0 1 


0 . 0 


0.0 


0.0 


-0.4 15 CD+00 


-0. 11 10D-01 


0.0 


0. 98C0D+01 


-0.143 OE+0 1 


-0. 1980D-0 1 


0.0 


0 . 0 


0.0 


0. 1000D + 01 


0.0 



ENTER THE MEASUREMENT SCALING MATRIX {" H"-M ATRIX} . 

DIMENSION = # OBSERVATIONS {NO} X # STATES {NS} 
THE ELEMENT H( 1, 1)= 
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THE 


EL EMENT 


H ( 


1, 


2) = 


THE 


ELEMENT 


H( 


1, 


3) = 


THE 


ELEMENT 


H ( 


1/ 


4) = 


TEE 


ELEMENT 


H( 


2, 


1) = 


THE 


ELEMENT 


H ( 


2/ 


2) = 


THE 


ELEMENT 


H ( 


2/ 


3) = 


THE 


ELEMENT 


H( 


2/ 


4) = 


THE 


ELEMENT 


H{ 


“3 

9 


1) = 


THE 


ELEMENT 


H( 


3/ 


2) = 


THE 


ELEMENT 


H( 


/ 


3) = 


THE 


ELEMENT 


H( 


3, 


4) = 


THE 


ELEMENT 


H( 


4, 


1) = 


THE 


ELEMENT 


H( 


4, 


2) = 


THE 


ELEMENT 


H ( 


4, 


3) = 


THE 


ELEMENT 


H( 


4, 


4) = 



THE MEASUREMENT SCALING MATRIX (" H" -MAT P IX} 



1 . 00000 
0 . 0 
0 . 0 
0. 0 



0.0 

1 . 00000 
0.0 
0.0 



0.0 

0.0 

1 . 00000 
0.0 



0.0 

0.0 

0.0 

1.00000 



DO YCU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT 
TYPE "YES” OR "NO". 



MEASUREMENT SCALING MATRIX H. 



0. 1 0 C0E + 0-1 
0. 0 
0. 0 
0.0 



0. 0 
0. 0 

0. 1000D + 0 1 

0. 0 



0.0 

0 .1000D+01 

0.0 



0.0 

0.0 

0. 1000D+01 



0*0 0*0 0*0 0*0 0*0 0*0 0*0 0*0 0*0 0*0 0*0 0*0 



ENTER THE OUTPUT MEASUREMENT COST MATRIX ("A "-MATRIX} 
.. 0NS x # OBSERVATIONS {NO} 



o 

6 

o 

6 

o 

6 

o 

6.25 



THE 


ELEMENT 


A( 


1, 


1) = 


THE 


EL EMENT 


A ( 


1/ 


2) = 


THE 


ELEMENT 


A { 


1, 


3) = 


THE 


ELEMENT 


A ( 


1/ 


4) = 


THE 


ELEMENT 


A ( 


2, 


1) = 


THE 


ELEMENT 


A ( 


2/ 


2) = 


THE 


ELEMENT 


A ( 


2/ 


3) = 


THE 


EL EMENT 


A ( 


2 / 


4) = 


THE 


ELEMENT 


A ( 


3/ 


1) = 


THE 


ELEMENT 


A ( 


- t 


2) = 


THE 


ELEMENT 


A ( 


3, 


3) = 


THE 


ELEMENT 


A ( 


3, 


4) = 


THE 


ELEMENT 


A ( 


4, 


1) = 


THE 


ELEMENT 


A ( 


4, 


2) = 


THE 


ELEMENT 


A ( 


4, 


3) = 


THE 


ELEMENT 


A ( 


4 / 


4) = 



THE OUTPUT 


MEASUREMENT 


COST MATRIX 


£" A"-MAT3 IX} .. . 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0 . 0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.25000 


DO YOU KISH 


TO CHANGE 


THE VALUE OF 


ANY MATRIX ELEMENT? 



TYPE "YES" CP "NO". 



OUTPUT COST MATRIX. 

0.0 0.0 
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0.0 



0.0 



0.0 

0.0 

0.2500D+00 



0.0 0.0 

0.0 0.0 

0 . C 0.0 



0 . 

0 . 

0. 



0 

0 

0 



6 

-> 

6.27 

S.8 



■j 

6 



INTER THE CONTRCI DISTRIBUTION MATRIX {»G"-MATRIX} . 



THE 


DIMENSION 
ELEMENT G ( 


i. 


# ^ STATES {NS} X # CONTROLS {NC} 


THE 


EL EMENT 


G( 


2, 


1) = 


THE 


ELEMENT 


G ( 


3, 


1) = 


THE 


ELEMENT 


G ( 


4, 


1) = 



THE CONTROI DISTRIBUTION MATRIX {"G"-MA TR IX} . . . 



0 . 0 

6.27000 
S. 8 COO 0 
0.0 



DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE n YES " OR "NO". 



ENTER THE CCNTRCI COST WEIGHTING MATRIX P'S" -MATRIX} 
DIMENSION = # CONTROLS {NC} X # CONTROLS £NC} 
THE ELEMENT B( 1, 1) = 



TEE CONTROL COST MATRIX B 

131.30000 



n 



DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" CR "NO". 



THE CONTROL DISTRIBUTION MATRIX G. . 

0.0 

0. 6270C+01 
0. 98C0D+01 
0 . 0 

TEE CONTROL COST MATRIX B.. 

0. 1313B+03 

EIGENSYSTEM OF OPTIMAL REGULATOR 

C-LOCP OPTIMAL REG. E- V ALUES. . . DE T ( SI- F+G *C) . . 
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-1. 23385D + 00 / 5. 54 5 4 6D-0 1 : - 4. 1 983 5D- 0 1 , 1 . 1 35 32D + 0 0 : 



C-1CCP EIGHT EIGENVECTOR MATRIX. 



M. 



-1. 0 19344D- 01 
1. 1 2969 ID-01 
1. 000030D+00 
-6. 742684D-01 



2.3087 17D-0 2 -8. 
-8.501 340D-02 -9. 
0 .0 1 . 
-3 .030 447D-0 1 -2. 



155484D-02 1. 

31 7336D-02 -1. 
000000D+ 00 0 



1222 64D-01 
3973 74D-01 
0 



865351 D- 01 -7. 748 4S9D-0 1 



CONTROL EIGENVECTOR MATRIX C*H. . 

-5.46431 4D- 03 2.109409D-02 2. 7139250-02 -1 . 6 7 63 34D-0 2 

C-LCCP OPT. REG. LEFT E-VECTOR MATRIX. .H-INV. . 

-3. 764753D+ 00 2.578703D+00 -3. 562309D-0 1 -1 . 01 022 CD+00 

-3. 4 2160 5D+ 01 - 9 .4 86 653D + 00 -4 . 604269D+00 -3 . 2452 6 1 D + 0 0 
3. 764753D+00 -2 .5 78703D+00 1.356231D+00 1.010220D+00 

1. 526591D+01 2 .4 1 9 663D+00 1.609196D+00 4.8415 48D-01 

THE CFTIMAL FEEDBACK GAIN CONTROL MATRIX. . .C=BINV*GT*S.. . 
-8. 5492E-0 1 -3.2475D-01 -8.53450-02 -4.3635D-02 



TEE CLOSED LOOP DYNAMICS MATRIX F-G*C. 



0.0 

-5. 3 6 033 7D+ 00 
1. 421803D+00 
0 . 0 



1 .0 00 0 00D + 00 
-2.45 1 197D+00 
-4 .6 12 572D + 00 
0.0 



0.0 0 . 
-5.4621 16D-01 -2. 
-8. 5617 86 D-0 1 -4. 
1.000000D+00 0. 



0 

73593 ID-01 
276 2 56D-0 1 
0 



DO YOU WISH TO OBTAIN A TIME RESPONSE 
OF THE SYSTEM YOU ARE EVALUATING? 

(Y OR N) 

NOTE: YOU MUST EE LOGGED ON AT A DUAL SCREEN 

(TEK 6 18) TERMINAL TO UTILIZE THIS MODS. 

THE F (SYSTEM) , G (CONTROL). H (OBSERVABLES), GAM (NOISE), 
A (OUTPUT COST) AND B (CONTROL COST) MATRICES FILL BE 
SAVED FOR REENTRY TC THE MAIN OPTSYS PROGRAM. 

y 



****** ******** * * ** * * ******* ****************** **** ****** ** *** * * 
IF YOU ARE DISSATISFIED WITH THE RESULTS 
THUS FAR AND WOULD LIKE TO EXIT TO CMS, 

-TYPE ’Y’ TO EXIT- 

(ANY OTHER INPUT TO CONTINUE) 

IfLif. IfLif. #£=£### jfcjjc SjCjJf # if. if. 

LOADING OPTCALC. ... 

EXECUTION BEGINS .. . 



DURING THIS SECTION OF THE PROGRAM YOU WILL: 

- SELECT THE TYPE CF SYSTEM RESPONSE TO PLOT 
(OPEN LOOP, CLOSED LOOP, OR FI LTER/REGU1ATOR ) 

- PROVIDE START AND STOP TIME FOR PLOTTING CALCULATIONS 
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SELECT THE TYPE CF DRIVING 
PROVIDE START AND STOP TINES FOR”TH£ 
EROVIDE DRIVING FUNCTION MAGNITUDE (S) 



FUNCTION (S) (STEP 
DRIVING 



OR RAH?) 
FUNCTION (S) 



CLEAR THE SCREEN TO CONTINUE 



THE F MATRIX 



0.0 


1.00000 


0.0 


0.0 


0.0 


-0.4 1500 


-0.01110 


0.0 


9.80000 


-1 . 43000 


-0.01980 


0.0 


0. 0 


0.0 


1.00000 


0.0 



THE G MATRIX 

0 . 0 

6.27000 
S. 80000 
0.0 



-0. 85452 



1. 00000 
0. 0 
0 . 0 
0. 0 



THE C MATRIX 
-0.32475 -0.08534 

THE H MATRIX 



0.0 

1 . 00000 
0.0 
0.0 



0.0 

0.0 

1.00000 

0.0 



-0.04364 



0.0 

0.0 

0.0 

1.00000 



THE K MATRIX 



0 . 0 


0 . 0 


0.0 


0.0 


0 . 0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 



THE FOLLOWING PLOTTING OPTIONS ARE AVAILABLE IF THE 
REQUIRED MATRICES WERE CALCULATED IN CPTSYSX: 

1. OPEN LOOP TIME RESPONSE 

XDOT = {F}*X + {G}*UC 

2. CLOSED LOOP TIME RESPONSE 

X COT = {F-G*C}*X + {G}*UC, U = {C} *X 

3. OPTIMIZED FILTER CLOSED LOOP SYSTEM RESPONSE. 

XDOT = {F}*X + {G}*UC, Z = { H} *X 

XHDOT = {F}*XH + {G}U + {K}*{Z - H*XH) 

4. FILTER + REGULATOR CLOSED LOOP SYSTEM RESPONSE. 

XEOT = {F+G*C}*X + f G} * UC. Z = {H} *X 

XHDOT = {F}*XH +{G}U + {K}*{Z - H*XH} , U = {C}*XH 

SELECT 1, 2, 3 OR 4. 



0 . 0 

-5. 36034 



THE AUGMENTED F MATRIX (F + G*C) 

1.00000 0.0 0.0 

-2.45120 -0.54621 -0.27359 
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1.42180 -4.61257 -0.85618 -0.42763 

0.0 0.0 1.00000 0.0 



AT WHAT TIKI DO YOU WANT TO START 
TEE TIME RESPONSE CALCULATIONS? 

INPUT START TIME IN SECONDS. (NORMALLY 0.0) 

6 



AT WHAT TIME DO YOU WANT TO STOP 
THE TIME RESPONSE CALCULATIONS? 

INPUT STOP TIME IN SECONDS. 

25 



THIS PROGRAM DIVIDES THE TIME INTERVAL YOU HAVE 
JUST SPECIFIED INTO UP TO 500 SMALL INTERVALS FOR 
THE INTEGRATION AND PLOTTING ROUTINES. IN ORDER 
TO SAVE COMPUTER TIME, THE NUMBER OF POINTS CAN BE 
CAN BE REDUCED WITH SOME LOSS IN CURVE FIDELITY. 

HOW MANY POINTS DO YOU WANT TO CALCULATE? 

500 

DOES THE SYSTEM UTIIIZE A DRIVING FUNCTION (CONI ROL INPUT)? 

( Y) ES OR ( N) O 
n 



DOES THE SYSTEM START WITH ALL INITIAL CONDITIONS =0.0? 



(Y) ES OR ( N) O? 
n 



*7 


WHAT 


IS 


THE 


INITIAL 


CONDITION 


FOR 


X ( 


11 


? 


6.02 






















WHAT 


IS 


THE 


INITIAL 


CONDITION 


FOR 


X ( 


2 ) 


• 


6 






















WHAT 


IS 


THE 


INITIAL 


CONDITION 


FOR 


X ( 


3) 


7 


6 




















6 


WHAT 


IS 


THE 


INITIAL 


CONDITION 


FOR 


X ( 


4) 


7 


THIS 


IS 


YOUR 


LAST OPPORTUNITY 


TO 










MAKE 


CHANGES 


IN I HE 


FOLLOWING 


AREAS. 








1 . 


SELECT 


ANOTHER 


TYPE OF 


SYSTEM TO 


PL 



(OPEN, CLOSED, FILTER OR FILTER /R EGULATOR) 

2. START AND STOP TIMES 

3. DRIVING FUNCTIONS 

4. INITIAL CONDITIONS 
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uvo 



5. CONTINUE 



SELECT A NUMBER BETWEEN 1 AND 5. 



TEE FOLLOWING INFORMATION IS PROVIDED 0 NLY 
FOR AN INDICATION OF PROPER PROGRAM OPERATION. 

AIL CONTROLS, STATES AND STATE ESTIMATES CAN BE PLOTTED. 



TIME 


0(1) 


0 . 0 


0.0 


0. 50 


0 . 0 


1.00 


0 . 0 


1. 50 


0 . 0 


2.00 


0 . 0 


2. 50 


0 . 0 


3. 00 


0 . 0 


3.50 


0 . 0 


4. 00 


0 . 0 


4.50 


0 . 0 


5. 00 


0 . 0 


5.50 


0 . 0 


6. 00 


0 . 0 


6. 50 


0 . 0 


7. 00 


0 . 0 


7.50 


0 . 0 


8. 00 


0 . 0 


8.50 


0 . 0 


9. 3 0 


0 . 0 


9.50 


0 . 0 


10. 00 


0 . 0 


1 0.50 


0 . 0 


1 1.00 


0 . 0 


1 1. 50 


0 . 0 


12. 00 


0 . 0 


12. 50 


0 . 0 


13. 00 


0 . 0 


13. 50 


0 . 0 


14. 00 


0 . 0 


14.50 


0 . 0 


15. 00 


0 . 0 


15.50 


0 . 0 


16. 00 


0 . 0 


16.50 


0 . 0 


17.00 


0 . 0 


17.50 


0 . 0 


18. 00 


0 . 0 


1 8. 50 


0 . 0 


19.00 


0 . 0 


19.50 


0 . 0 


20.00 


0 . 0 


20.50 


0 . 0 


21. 00 


0 . 0 


21.50 


0 . 0 


22. 00 


0 . 0 


22. 50 


0 . 0 


23.00 


0 . 0 


23. 50 


0 . 0 


24. 00 


0 . 0 


24.50 


0 . 0 


25.00 


0 . 0 



X(1) 

0.20000 00D-01 
0. 11031 08D^*0 1- 
■0.31 240 39D-02- 
■0. 12 142 95D-0 1- 
■0.13720 04D-01 
•0.97362 74D-02 
■0.34552 04D-02 
0.22423 63D-02 
0.56808 82D-02 
0.64850 95D-02- 
0.5231 7 22D-02- 
0.29374 93D-02- 
0. 59766 76D-03- 
■0. 1106274D-02- 
■0.19035 54D-02- 
■0.18757 04D-02 
■0. 13141 53D-02 
•0.56469 48D-03 
0.91051 25D-04 
0.4951 8 28D-0 3 
0.61687 27D-03- 
0.51774 99D-03- 
0. 30233 67D-03- 
0.73033 14D-04- 
■ 0. 98853 95D-0 4- 
■0. 1831320D-03- 
■0. 1853 8 26D- 03 
■0. 1330746D-03 
• C. 60089 64D-04 
0.53076 24D-05 
0. 46744 27D-04 
0.60412 64D-04 
0.51867 17D-04- 
0.31203 03D-0 4- 
0.85704 07D-05- 
■0.87893 15D-05- 
■0.17647 26D-04- 
■0.1 834 8 2 1D-0 4 
■0. 1348538D-04 
•0.63771 76D-05 
0.1 4426 43D-06 
0.43880 1 1D-05 
0.59036 99D-05 
0 . 51 86 8 79D-05- 
0.32109 31D-05- 
0.98082 65D-0 6- 
■0.76927 17D-06- 
■0. 16961 69D-05- 
0. 18131 70D-05 
0. 13640 03D-05 
■0.67 30920D-06 



X (2) 

0.0 0 
0 . 2 835 004D- 01 0 
0.249901 4D- 0 1 0 
•0. 1 044492D- 01 0 
0.341597 4D- 02 0 
0.1 1 3728 5D- 01-0 
0.1 27 6 48 2D- 01-0 
0.9465182D-02-0 
0.4 17680 ID-02-0 
0. 7495987D-03-0 
0.3909040D-02-0 
0.4 926 279D- 02 0 
•0.4 20285 ID-02 0 
0.2 52 4 967D- 02 0 
■0 .6 96430 ID-03 0 
0. 7052577D-03 0 
0. 1 421 627D-02 0 
0.1 48037 0D- 02-0 
0.1 08 99 6 4D-02-0 
0.5 17557 3D- 03-0 
0.8 50624 ID- 05-0 
0.3515887D-03-0 
0.4 74953 2D-03 0 
0.4 18178 ID-03 0 
0.2595 192D-03 0 
0.799512 8D- 04 0 
0.6 1 2523 0D- 04 0 
0. 1362740D-03 0 
0. 1460367 D- 03-0 
0. 1 10081 4D-03-0 
0.545142 ID-04-0 
0.2181 162D-05-0 
•0. 3284278 D- 04-0 
0. 463302 ID-04 0 
0.417515 6D- 04 0 
0.266291 0D- 04 0 
0.8 965253D-05 0 
0 . 5243 773D- 05 0 
0. 1305990D-04 0 
0.1 44071 8D-04-0 
0.1 11 1254D-04-0 
0.572330 6D- 05-0 
0.5212 149D-06-0 
0.304858 5D- 05-0 
0.451 0007D-05 0 
0.4 161 88 8D- 05 0 
0.272561 8D- 05 0 
0.990971 2D- 06 0 
0.4 36 470 8D- 06 0 
0.1 247888D-05 0 
0.1 419068D-05-0 



X (3) 



. 4433230D-0 1 
.79231 24D-01 
. 6 8503 92D-0 1 
. 2513011 E-0 1 
. 24639 73 D- 01 
. 5991416 D- 01 
. 71 72835D-01 
.62097 10D-0 1 
. 3954511D-01 
. 1411885 D- 01 
. 6308558D-02 
. 1 776 2 40 D-0 1 
.20100 80E-01 
. 1580986 D-01 
. 3417729 D- 02 
. 1131661 D- 02 
.39745 C7D- 02 
. 6 1 703 45 E-0 2 
. 5819287D-02 
. 38981 43D-02 
. 1 5 03 636D-02 
. 5066495 D-03 
. 16 80906D-02 
. 1 9 64 9 39 D-0 2 
. 1581498 D- 02 
. 8706565D-0 3 
. 1 5090 12D-03 
.3658826 D-03 
. 59934 13D-03 
. 57887 39 D-0 3 
. 39 74740D-03 
. 16 29461 D-0 3 
.38451 10D-04 
. 1 5964 41 D-03 
. 1929 683D-03 
. 15888 05 D-0 3 
. 90382 70D-04 
. 1 9 14584D-04 
. 33222 73D- 04 
.5799797 D-04 
. 57451 93D-04 
. 4042431 D-04 
. 1 750 6 G8D-0 4 
. 26 33 878C-0 5 
. 1510481 D-04 
. 1891 776D-04 
. 1 5 93 4 80 D- 04 
. 9 351 3 35 D- 0 5 
.2311 875D-05 
. 2984 895D-05 



X*** ******* ************************************************** 1 
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IF YOD ARE DISSATI SFI ED WITH THE RESULTS 
THUS FAR AND WOULD LIKE TO EXIT TO CMS, 

-TYPE ’ Y ’ TO EXIT- 

CANY OTHER INPUT TO CONTINUE) 

**************** ** ********* ****************** ** ** * ***** ** *** * * 



E 

C 

E 



12 0 ) 
12T 
122 ) 



P/O 

R/C 

R/O 



... Your Fortran program is now being loaded ... 
... execution will soon follow ... 
EXECUTION BEGINS.. . 



THIS PCRTICN OF THE PROGRAM PLOTS: 

- the states 

- EXTERNAL c6nTROL INPUTS, 

- F E EEEACK CONTROL INPUTS, 

- STATE ESTIMATES AND 

- RECONSTRUCTION ERRORS 

FROM THE DATA THAT YOU JUST CALCULATED. 

THE CAPABILITY IS ALSO AVAILABLE TO REVIEW ANY 
GRAPHS THAT YOU HAD PREVIOUSLY SAVED AS DATA 
FILES ON YOUR DISK. 

CLEAR THE SCREEN TO CONTINUE. 



THE FOLLOWING OPTIONS ARE AVAILABLE: 

1. PLOT THE DATA YOU JUST CALCULATED. 

2. PLOT A CURVE THAT YOU PREVIOUSLY SAVED. 

ENTER 1 OR 2 

i 



YOU MAY PLOT UP TO 4 SYSTEM VARIABLES VS TIME. 
HOW MANY VARIABLES DC YOU WISH TO PLOT? 

4 



WHICH TYPE OF VARIAEIE DO 



i 



YOU 

XI 



WISH TO PLOT AS CURVE NUMBER 1? 



1. STATE VARIABLE (IE., XI. X2, ETC) 

2. REEDBACK CONTROL (IE. , U = -C*X) 

3. CONTROL INPUT (IE., U1, U2, ETC.) 

4. STATE ESTIMATE (OBSERVES) (IE., XHAT1, aha 

5. STATE RECONSTRUCTION ERROR (IE., X1-XHAT1, 
X2-XHAT2, ETC.) 

ENTER 1,2, 3,4 OR 5 



XHAT2, ETC.) 



-> 

i 



WHAT IS THE SUBSCRIPT OF THE STATE VARIABLE THA 
YOU WANT TO PLOT AS THE NUMBER 1 CURVE VS TIME 
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•Oh-1 



IjJ*0 



WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 



NOTE: 



state yl 



1. 40 CHARACTERS MAX LENGTH 

2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 



IE. 



=> ALPHA 
=> BETA 
= > PHI 
=> THETA 



0 

1 



WHICH TYPE OF VARIAEIE DO YOU WISH TO PLOT AS CURVE NUMBER 2? 

1. STATE VARIABLE (IE., XI, X2, ETC) 

2. FEEDBACK CCNTRC1 (IE., U = -C*X) 

3. CONTROL INPUT (IE.. U1, U2 , ETC.) 

4. STATE ESTIMATE (OESERVER) (IE., XHAT1, XHAT2, ETC.) 

5. STATE RECONSTRUCTION ERROR (IE., X1-XHAT1, 

X2-XHAT2, ETC) 

ENTER 1,2, 3,4 OR 5 



•5 

2 



WHAT IS THE SUBSCRIPT OF THE STATE VARIABLE THAT 
YOU WANT TC PLOT AS THE NUMBER 2 CURVE VS TIME? 



WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 



KCTE 



state y2 



1. 40 CHARACTERS MAX LENGTH 

2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 



IE. 



=> ALPHA 
=> BETA 
= > PHI 
= > THETA 



WHICH TYPE OF VARIAEIE DO YOU WISH TO PLOT AS CURVE NUMBER 3? 



i 



i 

2 , 
3 , 

4 , 

5 



STATE VARIABLE (IE., XI, X2, ETC) 

FEEDBACK CONTROL (IE. , 0 = -C*X) 

CCNTRCL INPUT (IE., U1, U2, ETC.) 

STATE ESTIMATE (OBSERVER) (IE., XH AT 1 , XHAT 2 , ETC.) 
STATE RECONSTRUCTION ERROR (IE., X 1 - X HA T 1 , 

X2-XHAT2, ETC) 

ENTER 1,2, 3,4 OR 5 



WHAT IS THE SUBSCRIPT OF THE STATE VARIABLE THAT 
YOU WANT TC PLOT AS THE NUMBER 3 CURVE VS TIME? 



WHAT IS THE CURVE LABEL FOB THIS VARIABLE? 

NOTE: 1. 40 CHARACTERS MAX LENGTH 

2. GREEK SYMBOLS WILL BE PRINTED FOE ANY LETTERS 
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ENCLOSED IN PARENTHESES. 



IE. 



state y 3 



= > 
= > 
= > 
=> 



ALPHA 

BETA 

PHI 

THETA 



WHICH TYPE OF VARIABLE DO YOU FISH TO PLOT AS CJRVE NUMBER 4? 



0 

1 



!c 



1. STATE VARIABLE (IE., XI. X2, ETC) 

2. FEEDBACK CCNTRCI (IE. , U = -C*X) 

3. CONTROL INPUT (IE., U 1 , U2 , ETC. 

4. STATE ESTIMATE (OBSERVER) (IE., 

5. STATE RECONSTRUCTION ERFOR (IE., 
X2-XHAT2, ETC) 

ENTER 1,2, 3,4 OR 5 



HAT 1 , XHAI 2 , ETC.) 
X1-XHAT1, 



~> 

4 



W HAT IS THE SUBSCRIPT OF THE STATE VARIABLE THAT 
YOU WANT TO PLOT AS THE NUMBER 4 CURVE VS TIME? 



WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 



NOTE: 



state y4 



1. 40 CHARACTERS MAX LENGTH 

2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 



IE. 



= > 
= > 
= > 
= > 



ALPHA 

BETA 

PHI 

THETA 



3 



YOU MAY USE UP TO 3 HEADINGS. 

HOW MANY HEADINGS DO YOU DESIRE ON THIS GRAPH? 
0, 1, 2 OR 3 



WHAT IS THE DESIRED HEADING NUMBER 1? 



NOTE: 



1. 

2 . 



40 CHARACTERS MAX LENGTH 
GREEK SYMBOLS WILL BE PRINTED 
ENCLOSED IN PARENTHESES. 



FOR ANY LETTERS 



IE. 



closed loop system 



A 

E 

F 

Q) 



= > 
=> 
= > 
=> 



ALPHA 

BETA 

PHI 

THETA 



WHAT IS THE DESIRED HEADING NUMBER 2? 



NOTE: 



1. 40 CHARACTERS MAX LENGTH 

2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 



IE. 



example 3 



A) => ALPHA 

E) => BETA 

F) => PHI 

Q) => THETA 
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WHAT IS THE DESIRED HEADING NUMBER 3? 



NOTE: 



1 . 

2 . 



40 CHARACTERS MAX LENGTH 
GREEK SYMBOLS WILL BE PRINTED 
ENCLOSED IN PARENTHESES. 



FOR ANY LETTERS 



IE. 




A LPH A 
B ETA 
PHI 
THETA 



modern control theory' 

» USING A PRE-ALLOCATED DATASET FOR UNIT FT17F001 



THE FOLLOWING OPTIONS ARE AVAILABLE. 

1. BEGIN NEW GRAPH OF OTHER CONTROLS. STATES. 

2. REELOT PREVIOUSLY SAVED GRAPH DATA. 

3. EDIT THE CURRENT GRAPH. 

4. PLCT REVISED GRAPH ON THE TEK618. 

5. QUIT AND/OR MAKE METAFILE OF THE CUEVES. 
PREVIOUSLY SAVED. 



OR ESTIMATES. 



SELECT 

7 



A NUMBER EETWEEN 1 AND 5. 



DO YOU FA NT TO SAVE THE CURRENT GRAPH DATA TO 
EE USED LATER TO GENERATE A METAFILE? 

Y OR N 

NOTE: A METAFILE IS REQUIRED FOR SMOOTH VERSATEC PLOTS. 
THERE WILL BE AN OPPORTUNITY TO GENERATE A METAFILE 
JUST BEFORE EXITING THIS PROGRAM. 



WHAT FILE NAME DO 
closedlp 



YOU WANT THE CURVE 
(8 CHARACTERS MAX) 



DATA STORED UNDER? 



THE CURVE DATA IS BEING FILED UNDER CLOSEDLP DATA 
END OF DISSPLA 9.0 -- 16300 VECTORS GENERATED IN 1 PLOT FRAMES 
PROPRIETARY SOFTWARE PRODUCT OF ISSCO, SAN DIEGO, CA. 

1883 VIRTUAL STORAGE REFERENCES; 6 READS; 0 WRITES. 



THE FOLLOWING OPTIONS ARE AVAILABLE: 

1. MAKE METAFILE OF PREVIOUSLY SAVED EURVE. 

2. QUIT. 

ENTER 1 OR 2 



WHAT FILE NAME IS THE DATA STORED UNDER? 

closedlp 

THE CURVE DATA IS BEING LOADED FROM FILE CLOSEDLP DATA 
» USING A PRE-ALLOCATED DATASET FOR UNIT FT18F001. 



TEE FOLLOWING OPTIONS ARE AVAILABLE: 

1. MAKE METAFILE OF PREVIOUSLY SAVED E URVE. 

2. QUIT. 

ENTER 1 OR 2 
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✓ 

2 

END CF DISSPLA 9.0 — 16260 VECTORS GENERATED IN 1 PLOT FRAMES 
PROPRIETARY SOFTWARE PRODUCT OF ISSCO, SAN DIEGO , CA. 

1886 VIRTUAL STORAGE REFERENCES; 5 READS; 0 WRITES. 

EASD 121 DETACHED 
DASD 122 DETACHED 
DAS D 120 DETACHED 

************************************************************** 

DO YOU HA NT A VRSTEC PLOTTER SMOOTH COPY OF THE 
THE DISSPLA METAFILE THAT YOU JUST CREATED? 

(Y OR N) 

************** ********** ** *********************** ******** **** * 
1 (1201 E/O 

DASD 001 LINKED R/O; R/W BY MVS; R/O BY 0085P 
Z (0011 B/C - OS 
DASD 001 DETACHED 
CREATING NEW FILE: 

CREATING NEW FILE: 

PUN FILE 6749 TO MVS COPY 001 NOHOLD 

EASD 120 DETACHED 



$$$$$$$$$$$$$$ + £ $ i$c 

YOUR GRAPH (S) CAN BE PICKED UP AT THE COMPUTER CENTER. 

THE GRAPH ( S) WILL EE ADDRESSED TO "POP (USER ID)". 
************** ******************************************* ***** 



********* ** ***** ** ********* ****************** ** ** ******** *** * * 
DO YOU WANT TO 

1 . RUN OPTSYSX AGAIN 

2. RUN THE PL CT PROGRAM USING THE SAME MATRICES? 

(TO PLOT ANOTHER TYPE OF SYSTEM (OPEN/CLOSED)) 

3. QUIT 

-ENTER 1 , 2 OR 3 
3 



## 3fc# # #### $$$$ 3}e 2$C 3(c £ J*C $ ## ## # # # 



HAVE A GOOD DAY! ! 



************************************************************** 
R; T=9. 59/15.78 20:35:04 
record off 

END RECORDING OF TERMINAL SESSION 

The graphical output generated by this example follows 
as figure 3.2. 
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CLOSED LOOP SYSTEM 
EXAMPLE 3 

MODERN CONTROL THEORY 



o 
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940*0 0500 920*0 
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000*0 
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920*0- 090*0- 
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Figure 3.2 Closed-loop Time Response 



c 



FILTER CLOSED- LOOP SIMULATION 



The following filter s 
pp. 332 - 334]. 

In its present ccnfigur 
requires the input of a [ C 
regulator (if a [G] matri 
initiating the optimal esti 
] K j matrix evaluation. I 
program sequencing conven 
difficulties which may not 
system model, optimal filt 
by entering the identity m 
sequences requiring the ent 
matrix. 

The full terminal sessi 
input at the left margin i 
kelov each "?" . 



record on 

EEGIN RECORDING OF TERMINAL 
R; T= 0.0 1/0. 02 2 0:55:40 
opt sys 



THE OPTSYS EXEC CONTR 

1. OPTS YS X FORTRAN 

2. OPTCALC FORTRAN 

3. OPTPLOT FCRTRAN 

EACH PROGRAM PASSES INF 
THROUGH A DATA FILE WRI 
CASE, THESE FILES ARE " 
THE SIZE OF THESE FILES 
CAN USE ABOUT 20?? OF TH 
ENSURE THAT SUFFICIENT 

- TYPE "E" TO EXIT, 






YOU HAVE A DATA FILE N 
DISK THAT WAS PREVIOU 
PROGRAM AND CONTAINS 
MATRICES FROM THAT RUN 

IF YOU WOULD LIKE TO 
THIS RUN, THE CPTSYS 
DESIRED DATA AT THE AP 



imulation was taken from [Ref. 7 

ation, OPTSYSX program sequencing 
] matrix or design of an optimal 
x has been provided) , pricr to 
mator synthesis or user provided 
n order to comply with built-in 
tions, and circumvent program 
be specified in the particular 
er synthesis may be accomplished 
atrix [I] in those program input 
ry of an output cost (weighting) 

on is recorded below, with user 
n lower case letters or nu mbers 



SESSION 



OLS A TRIO OF PROGRAMS : 

(SYSTEM ANALYSIS) 

(CALCULATE TIME RESPONSE) 
(DISSPLA PLOTTING ROUTINE) 

ORMATION TO THE NEXT PROGRAM 
TTEN TO THE USEPS DISK. IN THIS 
OPTMAT DATA" AND "OPTPLOT DATA". 

VARY WITH THE SYSTEM ORDER, AND 
E USERS DISK SPACE. THEREFORE 
DISK SPACE IS AVAILABLE. 

ANY OTHER ENTRY TO CONTINUE - 



*******£*$****;}:****$:£ ****** ******$* 

AMED ’OPTMAT DATA’ ON YOUR A 
SLY GENERATED 3Y THE OPTSYS 
THE F, G, H, GAMMA, 1 AND B 



USE THESE SAME MATRICES FOR 
PROGRAM WILL READ IN THE 
PROPRIAIE TIME, 
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IF YOU TYPE (Y) ES . 

ANY OTHER INPUT KILL RESULT IN THAT FILE 3EING ERASED! 
*********** *************************************************** 



*********** *** ************************************************ 

DC YOU KANT THE NUMERICAL OUTPUT FROM OPTSYSX TO GO 
TO YOUR TERMINAL S(CREEN1 OR TO A D (ISK) FILE? 

(S OR D) 

************************************************************** 



********************************************************* *** * * 

OUTPUT RILL COME TO YOUR TERMINAL SCREEN. 
************************************************************** 

LOADING CPTSYS 

EXECUTION BEGINS... 

OPTSYSX IS A COMPLETELY INTERACTIVE OPTIMAL SYSTEMS CONTROL 
PROGRAM. IT WILL SOLVE NUMEROUS CONTROL PROBLEMS ON THE 
FOLLOWING TYPES OF SYSTEMS CONTROL EQUATIONS: 

X DOT = {F}*X + [G}*U + [GAM}* (W+WO) 

MEASUREMENT EQUATION — 

Z = [H}*X + {D} *U + V 

REGULATOR PERFORMANCE INDEX— 

J = 1/2 * INTEGRAL (Y *[A}*Y + U * {B} *U) DT 

STATE FEEDEACK GAIN DEFINITION — 

U = - {C} *X 

DO YOU WISH TO CONTINUE? TYPE "YES" CR "NO". 

Y — DATA ENTRY — 

ALTHOUGH OPTSYSX IS SPECIFICALLY DESIGNED TO READ 
ALL MATRIX DATA INTERACTIVELY, SEVERAL ALTERNATE 
METHODS ARE AVAILABLE TO USERS: 

METHOD 1 — THE "F","G", AND "GAMMA" MATRICES 
MAY BE READ FROM SEPARATE DATA FILES. 

METHOD 2 — THE "F","G", AND "GAMMA" MATRICES MAY BE 
EXPLICITLY DEFINED WITHIN SUBROUTINE "SETUP". 

[NOTE: IN EITHER CASE, THE USER SHOULD OBTAIN A COPY 
OF THE PROGRAM LISTING AND EXAMINE 
THE EXAMPLES CONTAINED IN S/E "SETUP".} 

DO YOU KISH TO CONTINUE? TYPE "YES" OR "NO". 

y 

DO YOU WISH TO INPUT THE "F", "G", AND "GAMMA" 

MATRICES FFCM SUBROUTINE "SETUP" IAW THE 
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METHOD DESCRIBED ON THE PREVIOUS SCREEN? 



n 



i 



TYPE "YES" OR "NO". 

GENERAL OPTSYSX OPTIONS: 

OPTION 1 -- SYSTEM ANALYSIS WITHOUT 

CPEN-L OOP EIGENS YSTEM CALCULATIONS. 

OPTION 2 — SYSTEM ANALYSIS WITH OPEN-LOOP 
EIGENS YSTEM CALCULATIONS. 

OPTION 3 — CPEN-L OOP EIGENSYSTEM FOUND 
AND PROGRAM TERMINATES. 

{"F"-MATRIX ENTRY FOLLOWS IMMEDIATELY.} 

OPTION 4 — MODAL DISTRIBUTION MATRICES COMPUTED 
WITHOUT FILTER OR REGULATOR SYNTHESIS 
CR STEADY-STATE ANALYSIS. 

SELECT AN OPTION: 1,2,3, OR 4. 



n 



DO YOU DESIRE RMS VALUES OF STATE AND CONTROL? 
TYPE "YES" CR "NO". 



OPTSYSX L QR/CLASS 



OPTION 


1 


— OPTIMAL FILTER 
SYNTHESIS WITH 
MATRIX INPUT. 


OPTION 


2 


— OPTIMAL FILTER 
SYNTHESIS WITH 
MATRIX INPUT. 


GPTION 


3 


— OPTIMA I FILTER 
SYNTHESIS WITH 
MATRIX INPUT.' 


OPTION 


4 


— OPTIMAL FILTER 
SYNTHESIS WITH 
MATRIX INPUT. 


SELECT 


AN 


OPTION: 1, 2, 3, 



ICAL OPTIONS: 

AND/OR P.EGULATOF 
NO EXTERNAL "C" OR "K" 



AND/OR REGULATOR 
EXTERNAL "C" 



AND/OR REGULATOR 
EXTERNAL "K" 



AND/OR REGULATOR 
EXTERNAL "C" AND "K" 



OR 4. 



DO YOU WISH TO DETERMINE THE STEADY-STATE RESPONSE 
FOR A CONSTANT DISTURBANCE? 

TYPE "YES" CR "NO". 

DO YOU WISH TO DETERMINE THE MODAL DISTRIBUTION 
AND GAIN MATRICES? 

TYPE "YES" CR "NO". 



OPE N-IOCP TRANSFER FUNCTION OPTIONS: 

OPTION 1 — NO OPEN-LOOP TRANSFER FUNCTIONS COMPUTED. 
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^Ui-9 



OPTION 2 -- POLES, RESIDUES, AND ZEROS COMPUTED. 
OPTION 3 — ONLY POLES AND ZEROS COMPUTED. 

OPTION 4 — ONLY POLES AND RESIDUES COMPUTED. 

SELECT AN OPTION: 1, 2, 3, OR 4. 

i 

NOISE TRANSFER FUNCTION OPTIONS: 

OPTION 1 — NO NOISE TRANSFER FUNCTIONS COMPUTED. 
OPTION 2 — POLES, RESIDUES, AND ZEROS COMPUTED. 
OPTION 3 — ONLY POLES AND ZEROS COMPUTED. 

OPTION 4 -- ONLY POLES AND RESIDUES COMPUTED. 

SELECT AN OPTION: 1, 2, 3, OR 4. 

i 

COMPENSATOR TRANSFER FUNCTION OPTIONS: 

OPTION 1 — NO COMP. TRANSFER FUNCTIONS COMPUTED. 

OPTION 2 — POLES, RESIDUES, AND ZEROS COMPUTED. 

OPTION 3 — ONLY POLES AND ZEROS COMPUTED. 

OPTION 4 — ONLY POLES AND RESIDUES COMPUTED. 

{NOTE: A COMPENSATOR TRANSFER FUNCTION CAN BE 
COMPUTED ONLY IF BOTH A REGULATOR 
AND FILTER ARE SYNTHESIZED 
ANE/OR INPUT.} 

SELECT AN OPTION: 1, 2, 3, OR 4. 

i 

WILL A FEED-FORWARD DISTRIBUTION MATRIX 
{"D" - MATRIX} EE INPUT ? 

TYPE "YES" OR "NO". 

HIS OPTION DETERMINES THE CRITERIA FOR DECIDING WHEN A 
ARKOV PARAMETER IS ZERO-THE MARKOV PARAMETER INDICATES 
HE ORDER OF THE NUMERATOR POLYNOMIAL OF EACH TRANSFER 
FUNCTION. 

ALL "N" ZEROS OF THIS POLYNOMIAL ARE PRINTED OUT AND 
THIS TEST TELLS HOW MANY EXTRA ROOTS EXIST AT Z = 0. 
LESS THAN 10. 0**{-IE} IS CONSIDERED ZERO. 

THE DEFAULT VALUE OF THIS PARAMETER {IE} IS 6. 

IN OTHER WORDS, IE = 1.0E-6. 

IF YOU DESIRE A DIFFERENT MARKOV CRITERIA, 

TYPE THE INTEGER VALUE. 

IF YOU DESIRE THE DEFAULT VALUE, TYPE "0" {ZERO} 
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K )*0 



0 



DO YCU DESIRE TO SYNTHESIZE A STABLE FILTER [OR REGULATOR} BY 
DESTABILIZING THE ORIGINAL SYSTEM? 

[NOTE: WORKS FOR FILTER OR REGULATOR BUT NOT FOR BOTH 
IN THE SAME RUN.} 

TYPE "YES" CR "NO". 

DC YOU DESIRE TO PRINT THE EULER-LAGRANGE El GEN SYSTEM 
PRIOR TO DECOMPOSITION [FOR CHECKING THE PROGRAM}? 

TYPE "YES" CR "NO". 



POWER SPECTRAL DENSITY (PSD} OPTION 1 : 

OPTION 1 — COMPUTE THE PSD OF THE OUTPUTS AND/3 R THE 

CONTROLS OF THE CONTROLLED SYSTEM WHEN FORCED BY 
PROCESS AND MEASUREMENT NOISE. [NOTE: BOTH A 
REGULATOR AND A FILTER MUST BE RESIDENT IN THE 
PROGRAM TO USE THIS OPTION.} 

OPTION 2 — SAME AS OPTION 1 ABOVE BUT ONLY PRINT THE 
RESIDUES OF EACH TRANSFER FUNCTION 
USED IN THE PSD COMPUTATION. 

OPTION 3 -- NOT DESIRED. 



SELECT AN OPTION: 1, 2, OR 3. 



DO YOU DESIRE REGULATOR SYNTHESIS ONLY? 
TYPE "YES" CR "NO". 



ENTER THE # OF STATES (NS} OF THE SYSTEM MAI RIX 
{"F"-M ATRIX} . 



ENTER THE # OF CONTROLS [NC} OF THE CONTROL SYSTEM MODEL 
{"G"-MATRIX} . 

i 

ENTER THE # OF MEASUREMENTS OR OBSERVATIONS {NO} OF THE 
{"H"-MATRIX} . 

i 



•? 

i 



IOL 

0 



ENTER THE # OF PROCESS NOISE SOURCES [NG} OF THE 
{" GAM MA"-M ATRIX} . 

FLAG/ PA FAME TER SETTINGS FOR THIS RUN ARE AS FOLLOWS: 
IQ IE ISS IM ITF1 ITF2 ITF3 IF'DFW IE IDEBUG 
02000 0 0 0 00 



0 
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ISET IDSTAE IPSD IYU INORM IREG NS NC NOB NG 
0 0 0 0 0 0 2 1 1 1 

ORDER OE SYSTEM = 2 

NOMEER CF CONTROLS = 1 

NUMBER CF OBSERVATIONS = 1 

NUMBER CF PROCESS NCISE SOURCES = 1 

ENTER THE SYSTEM MATRIX { , 'F"-MATRIX} 

DIMENSION = # STATES {NS} X # STATES {NS} 

THE ELEMENT F( 1, 1)= 

6 

THE ELEMENT F ( 1, 2) = 

i 

THE ELEMENT F ( 2, 1) = 

6 

THE ELEMENT F{ 2, 2) = 

0 

-4. 6 

THE SYSTEM MATRIX {» F" -MATRIX} .. . 

0.0 1.00000 

0. 0 -4.60000 

DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 

TYPE "YES" CR "NO", 
n 

CPE N LOOP DYNAMICS MATRIX F. . 

0.0 0.10 00E+0 1 

0.0 - 0.46 00D+0 1 

ENTER THE MEASUREMENT SCALING MATRIX {" H"-M ATR IX} . 

DIMENSION = # OBSERVATIONS {NO} X # STATES {NS} 
THE ELEMENT H ( 1, 1) = 

1 

THE ELEMENT H( 1, 2 ) - 

6 

THE MEASUREMENT SCALING MATRIX £ "H"-MAr RIX} . . . 

1 . 00000 0.0 

DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 

TYPE "YES" CR "NO", 
n 
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0*0 



MEASUREMENT SCALING MATRIX 



H 



O.IOCOE+OI 0.0 

ENTER THE OUTPUT MEASUREMENT COST MATRIX {"A "-MATRIX} . 
DIMENSION = # OBSERVATIONS {NO} X # OBSERVATIONS (NO) 
THE ELEMENT A ( 1, 1) = 

i 

THE OUTPUT MEASUREMENT COST MATRIX {" A"-MATR IX} . . . 

1. 00000 

DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 

TYPE "YES" OR "NO". 

C 

OUTPUT COST MATRIX A.. 

0. 10C0D+01 

ENTER THE CONTROL DISTRIBUTION MATRIX {"G»- M A TP. IX} . 

DIMENSION = # STATES {NS} X # CONTROLS {NC} 

THE ELEMENT G ( 1 , 1) = 

THE ELEME'NT G( 2, 1) = 

6.787 

THE CONTROL DISTRIBUTION MATRIX ("G "-MA T RIX} . . . 

0.0 

0. 78700 

DC YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 

TYPE "YES" OR "NO", 

n 

ENTER THE CONTROL COST WEIGHTING MATRIX f"3" -MATRIX} 
DIMENSION = # CONTROLS {NC} X # CONTROLS {NC} 

THE ELEMENT B ( 1 , 1) = 

i 

THE CONTROL COST MATRIX B. . . 

1. 00000 

DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" OR "NO". 
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THE CONTROL DISTRIBUTION MATRIX 



G 



0.0 

0.787 OD+OO 

THE CONTROL COST MATRIX B 

0. 10C0D+01 



EIGENSYST2K OF OPTIMAL REGULATOR 

C-LOCP CPTIMAL REG. E-VALUES. .. DET(SI- F+G*C) . . 
-1.7 12C6E-01:-4. 596 8 1D+00: 



C-LOCP RIGHT EIGENVECTOR MATRIX 



9. 85658 8D-01 -2 . 1 257C3D-0 1 
-1.687503D-01 9.771458D-01 



M 



CONTROL EIGENVECTOR MATRIX C*M.. 

-9.4963 1SD-0 1 3.S57155D-03 



C-LCCP OPT. REG. LEFT E-VECTOR MATRIX. . M- INV. . 

1 . 0 5379 8D+ 00 2.292453D-01 

1. 8 1 987 9D- 01 1 .062 S79D+00 



THE CPTIMAL FEEDBACK GAIN CONTROL MATRIX. . . C=BIN V*GT*S . . . 
-1.00CCD + 00 -2.134 9E-0 1 



TEE CLOSED LOOP DYNAMICS MATRIX F-G*C. . 

0.0 1 .0 00 OCOD+OO 

-7. 8 700 00D- 01 -4.768018D+00 



ENIER THE PROCESS NOISE DISTRIBUTION 
MATRIX {"GAMMA"-MATRIX} . 

DIMENSION = # STATES {NS} X # PROCESS NOISE SOURCES {NG} 
THE ELEMENT GAM ( 1, 1) = 



THE ELEMENT GAM ( 2, 1) = 

6. i 



THE PROCESS NOISE DISTRIBUTION MATRIX 
{"GAM M A M - M ATR IX} . . . 

0.0 

C. 10000 



DO YOU NISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" OR "NO". 
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ENTER THE PROCESS NOISE PSD WEIGHTING MATRIX 
f"Q"MATRIX} . 

DIMENSION = # PROCESS NOISE SOURCES {NG} X 
#PROCESS NOISE SOURCES {NG} 

THE ELEMENT Q ( 1, 1) = 

io 

THE PROCESS NOISE WEIGHTING MATRIX Q.. 

1 C. 00000 

DO YOU WISH TO CHANGE THE VALUE OF ANY MATRE X ELEMENT? 

TYPE "YES" CR "NO", 

n 

PROCESS NOISE DISTRIBUTION MATRIX GAMMA.. 

0.0 

0. 10C0D+00 

POWER SPECTRAL DENSITY - PROCESS NOISE.. ..Q.. 

0. 1 0 00D + 02 

ENTER TEE MEASUREMENT NOISE DISTRIBUTION MATRIX {"R"MATRIX} . 

DIMENSION = # OBSERVATIONS {NO} X # OBSERVATIONS {NO} 

THE ELEMENT R( 1, 1) = 

6. ooooooi 

THE MEASUREMENT NOISE DISTRIBUTION MATRIX R. .. 

0. 00000 

DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 

TYPE "YES" CR "NO", 

n 

POWER SPECTRAL DENSITY- MEASUREMENT NOISE.. R.. 

0. 1 0C0D-06 

ENTER THE FEEDBACK GAIN ESTIMATOR MATPIX {"K"- MATRIX} . 

DIMENSION = # STATES {NS} X # OBSERVATIONS {NO} . 

THE ELEMENT K( 1, 1)= 

95. 4 

THE ELEMENT K( 2, 1)= 

456 1 

THE FEEDBACK GAIN ESTIMATOR MATRIX {"K" -MATRIX} 
95.40000 
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4561. COOOO 



DO YOU KISH TO CHANGE THE VALOE OF ANY MATRIX ELEMENT? 
TYPE "YES" OR "NO". 

FILTER STEADY STATE GAINS K. . . 

9.5400C0D+01 
4.5610 00D+ 03 



THE CLOSED LOOP FILTER DYNAMICS MATRIX IS.... 

-9. 54 00 00D+ 01 1 .0 00 OCOD+OO 

-4. 56 1000D + 03 -4 .6 00 000D+00 



EIGENSYSTEM OF OETIMAL ESTIMATOR 

C-IOC? SUBOPT. EST. E-VALUES. .. DE T (SI-F+K*H) . . 
-5. OOOOOD+OI, 4.99984D+01: 



C-LCCP RIGHT EIGENVECTOR MATRIX 



9. 9 53957D-03 - 1 .0 9 62 16D-02 
1 . 0000 0CD+ 00 0 .0 



M 



MEASUREMENT EIGENVECTOR MATRIX H (BAR) *M. . 

9.953957D-03 - 1 . 09 62 16D -02 

C-IOCP SUBOPT. FILTER LEFT E-VECTOR MATRIX. . M- IN V . . 

0.0 1 .000 OCOD+OO 

-9. 1 22292D+C1 9.080291D-01 



THE COVARIANCE OF THE ESTIMATION ERROR P.. 

7. 1 5C503D-06 2.271000D-04 

2. 27 1000D- C4 1.181151D-02 



RMS VALUES OF THE ESTIMATION ERROR 

2. 6 7 404 2D- 03 1 .0868C8D-01 

DO YOU WISH TO OBTAIN A TIME RESPONSE 
OF THE SYSTEM YOU ARE EVALUATING? 

(Y OR N) 



NOTE: YOU MUST BE LOGGED ON AT A DUAL SCREEN 
(TEX 618) TERMINAL TO UTILIZE THIS MODE. 

THE F (SYSTEM) , G (CONTROL), H (OBSERVABLES), GAM (NCISE), 
A (CUTPUT COST) AND B (CONTROL COST) MATRICES WILL BE 
SAVED FOR REENTRY TO THE MAIN OPTSYS PROGRAM. 



y 
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*********** *** ************************************************ 
IF YOU AftE DISSATISFIED WITH THE RESULTS 
THUS FAR AN E WOULD LIKE TO EXIT TO C MS, 

-TYPE »Y» TO EXIT- 

(ANY OTHER INPUT TO CONTINUE) 

************************************************************** 

LOADING OPTCALC 

EXECUTION BEGINS .. . 

DURING THIS SECTION OF THE PROGRAM YOU WILL: 

- SELECT THE TYPE OF SYSTEM RESPONSE TO PLOT 
(OPEN LOOP, CLOSED LOOP, OR FILTER/REG ULATORl 

- PROVIDE START AND STOP TIME FOR PLOTTING CALCULATIONS 

- SELECT THE TYPE CF DRIVING FUNCTION (S) (STEP OR RAMP) 

- PROVIDE START AND STOP TIMES FOR THE DRIVING FUNCTION (S) 

- PROVIDE DRIVING FUNCTION MAGNITUDE (S) . 

CLEAR THE SCREEN TO CONTINUE 



THE F MATRIX 

0.0 1.00000 

0.0 -4.60000 



THE G MATRIX 

0.0 

0. 78700 



THE C MATRIX 
-1.00000 -0.21349 



THE H MATRIX 

1.00000 0.0 



THE K MATRIX 

95.40000 
456 1. 00000 



THE FOLLOWING PLOTTING OPTIONS ARE AVAILABLE IF THE 
REQUIRED MATRICES WERE CALCULATED IN OPTSYSX: 

1. OPEN LOOP TIME RESPONSE 

XDCT = {F} *X + £G] *UC 

2. CLOSED LOOP TIME RESPONSE 

XDCT = (F-G*C}*X + {G}*UC, U = {C} *X 

3. OPTIMIZED FILTER CLOSED LOOP SYSTEM RESPONSE. 

XDOT = {F}*X + (G}*UC, Z = {HI *X 

XHDOT = {F} *XH + {G}U + {K} * {Z - H*XH} 

4. FILTER + REGULATOR CLOSED LOO? SYSTEM RESPONSE. 

XDCT = fF + G*C}*X + {G}*UC, Z = {H} *X 

XHDOT = {F} *XH +{G}U + {K}*{Z - H*XH}, U = {C} *XH 
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SELECT 1 , 2, 3 OR 4 



7 

3 



95.40000 
4561. 00000 



THE (K*H) MATRIX 

0.0 

0.0 



0.0 

0.0 

95.40000 
4561. 00000 



THE C CMBINED SYSTEM F 

1.00000 0.0 

-4.60000 0.0 

0.0 -95.40000 

0.0 -4561.00000 



MATRIX (2*NS X 2*NS) 

0.0 

0.0 

1 . 00 000 
-4.60000 



THE AUGMENTED G MATRIX (2*NS X NC| 

0.0 

0. 78700 

0.0 

0. 78700 



AT WHAT TIME DO YOU WANT TO START 
THE TIME RESPONSE CALCULATIONS? 

INPUT START TIME IN SECONDS. (NORMALLY 0.0) 



AT WHAT TIME DO YOU WANT TO STOP 
TEE TIME RESPONSE CALCULATIONS? 

INPUT STCP TIME IN SECONDS. 

6.3 



500 

DOES 

y 



THIS PROGRAM DIVIDES THE TIME INTERVAL YOU HAVE 
JUST SPECIFIED INTO UP TO 500 SMALL INTERVALS FOR 
THE INTEGRATION AND PLOTTING ROUTINES. IN ORDER 
TC SAVE COMPUTER TIME, THE NUMBER OF POINTS CAN BE 
CAN BE REDUCED WITH SOME LOSS IN CURVE FIDELITY. 

HOW MANY POINTS DO YOU WANT TO CALCULATE? 



THE SYSTEM UTILIZE A DRIVING FUNCTION (CONTROL INPUT)? 

(Y) ES OR (N) 0 



TWO TYPES CF FUNCTIONS CAN BE USED AS DRIVERS. 

1. STEP INPUT 

2. RAMP INPUT 

ENTER YOUR SELECTION, 1 OR 2. FOR DRIVING FUNCTION NUMBER 1 

i 



AT WHAT TIME DO YOU DESIRE INPUT NUMBER 1 TO START? 
INPUT THE START TIME IN SECONDS. 
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Ul -0 0*0 0*0 



-> 

6 



6.4 



r> 

-10 



AT WHAT TIME DO YOU DESIRE INPUT NUMBER 1 TO STOP? 
INPUT THE STOP TIME IN SECONDS. 



WHAT IS THE MAXIMUM VAIUE OF 
DRIVING FUNCTION NUMBER 1 ? 



DOES THE SYSTEM START WITH ALL INITIAL CONDITIONS = 0.0 ? 

(Y ) ES OR (N) 0? 



6.i 

6 



.5 



WHAT IS THE INITIAL CONDITION FOR X ( 1 ) ? 

WHAT IS THE INITIAL CONDITION FOR XHAT ( 1) ? 

WHAT IS THE INITIAL CONDITION FOR X ( 2) ? 

WHAT IS THE INITIAL CONDITION FOR XHAT ( 2) ? 



THIS IS YOUR LAST OPPORTUNITY TO 
MAKE CHANGES IN THE FOLLOWING AREAS. 

1. SELECT ANOTHER TYPE OF SYSTEM TO PLOT 
(OPEN, CLOSED, FILTER OR FILTER/REGULATOR) 

2. START AND STOP TIMES 

3. DRIVING FUNCTIONS 

4. INITIAL CONDITIONS 

5. CONTINUE 

SELECT A NUMEER BETWEEN 1 AND 5. 



THE FOLLOWING INFORMATION IS PROVIDED ONLY 
FOR AN INDICATION OF PROPER PROGRAM OPERATION. 

AIL CCNTROLS, STATES AND STATE ESTIMATES CAN BE PLOTTED. 

'IME U(1) X(1) X (2) X ( 3) 



0.0 0.0 

0.01 - 0 . 10 00 000 D + 02 
0.01 -0. 10 000 00D+02 
0.02 -0. 1000000D+02 
0.02 -0. 1000000D+02 
0.03 -0. 1 OOOOOOD+02 
0. 04 -0. 1 0000 00D+ 02 
0. 04 -0. 1000000D+02 



0. 1000000D+00 0. 
0. 1028 186D+00 0. 
0. 1052810D+00 0. 
0. 1073 970D+00 0. 
0. 1091 759D+00 0. 
0. 1106 270D+00 0. 
0. 1 117 590D+00 0. 
0. 1 125809D+00 0. 



5000000D+00 
43981 44D+00 
3812672D+00 
3243 1 38D+00 
26891 08D+00 
21501 60D+00 
1 6258 84D +00 
11 15880D+00 



0, 

0, 

0. 

0 

0. 

0 , 

0 . 



497 3 46 3D- 0 1 
8502367D-01 
1 078565D+00 
1209 43 8D+ 0 0 
1 27031 8D+ 00 
128521 2D+ 0 0 



0.1 27303SD+00 
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0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0. 

0 . 

0 . 

0 . 

0. 

0. 

0. 

0. 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 



05 

05 

06 
07 

07 

08 
08 

09 
1 0 

10 
1 1 
1 1 
12 

13 
1 3 

14 

14 

15 
1 6 

16 
17 
17 
1 8 

19 

1 9 

20 
20 

2 1 
22 
22 
23 

23 

24 

25 

25 

26 
26 

27 

28 
28 
29 

29 

30 



-0. 1000000D+02 
-0. 1000000D+02 
-0.1 000000D+02 
-0.1 0000 0 0D+ 02 
-0.1 0000 0 0 D+0 2 
-0.1 000000D+02 
-0.1 0000 00 D+02 
-0. 1 0000 00D+ 02 
-0.10 000 00 D+0 2 
-0.1 000000D+02 
-0. 1000000D+02 
-0. 1000000D+02 
-0.1 0000 00 D+0 2 
-0.1 OOOOOOD+02 
-0.1 000000D+02 
-0. 1000000D+02 
-0. 1000000D+02 
-0. 1000000D+02 
-0.1 000000D+02 
-0. 1 0000 00D+ 02 
-0 .1000000D+02 
-0 .1000000D+02 
-0.1 000000D+02 
-0.10 000 00D+ 02 
-0.1 0000 00 D+0 2 
-0.10 COO 00D+02 
-0.10 000 00 D+02 
-0. 10C00 00D+02 
-0.1 0C0000D+02 
-0. 1000000D+02 
-0.1 0C0000D+02 
-0. 10C0000D+02 
-0. 1 000000D+02 
-0. 1 0C0000D+02 
- 0.10 000 00 D +02 
-0.10 COO 00D+ 02 
-0. 10 COO 00 D+0 2 
-0. 10000 00D+02 
-0.1 000000D+02 
-0. 1000000D+02 
-0.10 000 00 D+ 0 2 
-0.1 OOOOOOD+02 
-0.1 0 000 00 D+0 2 



C. 1 131 009D+00 
0. 1 133 273D+00 
0. 1 132 68 1D+00- 
0. 1 129310D+00- 
0. 1 123 237D+00- 
0. 1 114535D+00- 
0. 1103 275D+00- 
0. 1089 528D+00- 
0. 1 0 73 36 OD+O 0- 
0. 1054 837D+00- 
0. 1034 025D+0 0- 
0. 1010984D+00- 
0. 9857 769D-0 1- 
0. 9584 6 1 0D- 01- 
C. 9 290 943D-0 1- 
0.8977 326D-01- 
0.8644301D-01- 
0. 829 2 3 98D- 0 1- 
0.7922130D-01- 
0.7533997D-01- 
0. 7128485 D- 01- 
0.6706 068D-01- 
0.6267 205D-01- 
C. 5812 345D-01- 
0. 534 1 923 D-0 1- 
0.4856 362D-01- 
0. 4356 075D-01- 
0.3841 463 D-0 1- 
0. 331 2 915D-0 1- 
0.2770 81 1D-0 1- 
0. 2215 519D-0 1- 
0. 1647 400 D-0 1- 
0. 1066 802D-0 1- 
0.4740642D-02- 
-0. 1304821D-02- 
-0.7465 155D-02- 
-0. 1 373 724D- 0 1- 
-0.2011 802D-01- 
-0.2660 455D-01- 
-0.3319394D-01- 
-0. 3988340 D- 01- 
-0.4667 020D-01- 
-0.5355 170D-01- 



0. 61 97598D-01 
0. 137145 1D-0 1 
0. 33 233 1 5D-01 
0. 78902 78D-0 1 
0. 1233292D + 00 
0. 1665461D+00 
0. 2085866D+00 
0. 2494 827D+00 
0. 2892654D+00 
0. 32796 52D + 00 
0. 36561 14D + 00 
0. 4022329D+00 
0. 437 8574D + 00 
0. 4725121D+00 
0. 5062234D+00 
0.539017 OD + OO 
0. 57091 78D+00 
0. 601 95 03D + 00 
0. 6321380D + 00 
0 . 66 1 50 38D+00 
0.6900703D+00 
0. 71 78591D+00 
0. 74 4 89 1 4D+00 
0. 771 1879D+00 
0. 7967685D+00 
0. 821 6527D+00 
0. 8458595D+00 
0.8694073D+00 
0. 8923141D+00 
0. 9145973D+00 
0. 936273 9D+00 
0. 95 736 04D +00 
0. 97787 2 9D+00 
0. 9978270D+00 
0. 10 1 7238D+0 1 
0. 10361 20D+0 1 
0. 10544 89D + 01 
0. 1072357D+01 
0. 1 0 8973 9D + 01 
0. 1 106648D+01 
0. 11230 96D + 01 
0. 1 139097D+01 
0. 11546 62D + 0 1 



0. 


, 124741 6D + 


00 


0. 


, 121724 2D+ 


00 


0. 


, 1 187650D+ 


00 


0. 


, 1 161 034D + 


00 


0. 


, 1 137974D+ 


00 


0. 


1 1 17983D + 


00 


0. 


. 1 100069D+ 


00 


0. 


, 1083096D + 


00 


0. 


,106601 6D+ 


00 


0. 


, 1047973D+ 


00 


0. 


, 1C28 33 8D+ 


00 


0. 


, 1 006703D+ 


00 


0. 


, 9 828370D- 


01 


0. 


, 9566496D- 


0 1 


0. 


.9281 438D- 


01 


0. 


, 8973812D- 


01 


0. 


, 864454 5D- 


01 


0. 


,829467 1D- 


01 


0. 


,792521 3D- 


0 1 


0. 


,7537 11 4D- 


01 


0. 


.713120 5D- 


01 


0. 


.6708207D- 


0 1 


0. 


. 6 268 74 1D- 


01 


0. 


, 5813344D- 


01 


0. 


, 5342495D- 


01 


0. 


. 4 856 623D- 


01 


0. 


.4 35613 1D- 


01 


0. 


, 3841 398D- 


01 


0. 


. 3312793D- 


01 


0. 


, 2770 67 3D- 


01 


0. 


. 22 1 5 39 2D- 


01 


0. 


. 1647295D- 


01 


0. 


. 1066723D- 


01 


0. 


,4 740 10 6D- 


02 


-0. 


. 130514 9D- 


02 


-0. 


. 746532 6D- 


02 


-0. 


. 1373 73 CD- 


01 


-0. 


.201180 1D- 


01 


-0. 


.266045 0D- 


01 


-0 , 


. 33 1 938 8D- 


01 


- 0. 


. 3988334D- 


01 


-0. 


.466701 5D- 


0 1 


- 0. 


.5355 166D- 


01 



*** * * ** ******* ** ** ******** ***************** ** ************ *** ** 



IF YOU ARE DISSATISFIED WITH THE RESULTS 
THUS FAR AND WOULD LIKE TO EXIT TO CMS, 

-TYPE ’ Y ’ TO EXIT- 

(ANY OTHER INPUT TO CONTINUE) 

************************************************************** 



E 

C 

E 



12 0) 
121* 
122 



E/O 

E/O 

E/O 



... Your Fortran program is now being 
... execution will soon follow 
EXEC U TIC N BEGINS... 



loaded . . . 



THIS PORTION OF THE PROGRAM PLOTS: 
- THE STATES , 
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- EXTERNAL CONTROL INPUTS, 

- FEEEEACK CONTROL INPUTS, 

- STATE ESTIMATES AND 

- RECONSTRUCTION ERRORS 

FROM THE DATA THAT YOU JUST CALCULATED. 

THE CAPABILITY IS ALSO AVAILABLE TO REVIEW ANY 
GRAPHS THAT YOU HAD PREVIOUSLY SAVED AS DATA 
FILES ON Y CUR DISK. 

CLEAR THE SCREEN TO CONTINUE. 



THE FOLLOWING OPTIONS ARE AVAILABLE: 

1. PLOT THE DATA YOU JUST CALCULATED. 

2. PLCT A CURVE THAT YOU PREVIOUSLY SAVED. 

ENTER 1 OR 2 

•? 

i 



YOU MAY PLCT UP TO 4 SYSTEM VARIABLES VS TIME. 
HOW MANY VARIABLES DO YOU WISH TO PLOT? 

£ 

2 



WHICH TYPE OF VARIAEIE DO YOU WISH TO PLOT AS CURVE NUMBER 1? 



•? 

i 






1. STATE VARIABLE (IE., XI, X2, ETC) 

2. FEEDBACK CONTROL (IE. . U = -C*X) 

3. CONTROL INPUT (IE., U1, U2 , ETC. 

4. STATE ESTIMATE (OBSERVER) (IE., 

5. STATE RECONSTRUCTION ERROR (IE., 
X2-XHAT2, ETC) 

ENTER 1,2, 3,4 OR 5 



H AT 1 , XHAI2, ETC.) 
X 1-XHAT1 , 



WHAT IS THE SUBSCRIPT OF THE STATE VARIABLE THAT 
YOU WANT TC PLOT AS THE NUMBER 1 CURVE VS TIME? 

i 



WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 

NOTE: 1. 40 CHARACTERS MAX LENGTH 

2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 

IE. (A) => ALPHA 
E) => BETA 
i Fi => PHI 
'£) => THETA 

angular position - (x) 1 



WHICH TYPE OF VARIABLE DO YOU WISH TO PLOT AS CURVE NUMBER 2? 

1. STATE VARIABLE (IE., XI. X2, ETC) 

2. FEEDBACK CCNTRCL (IE. , U = -C*X) 

3. CCNTRCL INPUT (IE., U 1, U2, ETC.) 

4. STATE ESTIMATE (OBSERVER) (IE., XHAT1, XHAI 2 , ETC.) 

5. STATE RECONSTRUCTION ERROR (IE., X1-XHAT1, 

X2-XHAT2, ETC) 
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■7 

4 



ENTER 2, 3/4 OR 5 



”7 

i 



WHAT IS THE SUBSCRIPT OF THE STATE ESTIMATE THAT 
YOU WANT TO PLOT AS THE NUMBER 2 CURVE VS TIME? 



WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 



NOTE: 1. 40 CHARACTERS MAX LENGTH 

2. GREEK SYMBOLS WILL BE PRINTED 
ENCLOSED IN PARENTHESES. 

A LPH A 
BETA 
PHI 
THETA 
- (x)el 



FOR ANY LETTERS 



angular position 




-7 

3 



YOU MAY USE UP TO 3 HEADINGS. 

HOW MANY HEADINGS DO YOU DESIRE ON THIS GRAPH? 
0/ 1, 2 OR 3 



WHAT IS THE DESIRED HEADING NUMBER 1? 



NOTE: 



1. 

2 . 



40 CHARACTERS MAX LENGTH 
GREEK SYMBOLS WILL BE PRINTED 
ENCLOSED IN PARENTHESES. 



FOR ANY LETTERS 



IE. 



filter only closed 



;a) => 
e; => 

F) => 

JQ) => 

loop 



ALPHA 

BETA 

PHI 

THETA 



WHAT IS THE DESIRED HEADING NUMBER 2? 



NOTE: 



1. 

2 . 



40 CHARACTERS MAX LENGTH 
GREEK SYMBOLS WILL BE PRINTED 
ENCLOSED IN PARENTHESES. 



FOR ANY LETTERS 



IE. 



example 4. 1 




ALPHA 
BETA 
PHI 
I HETA 



WHAT IS THE DESIRED HEADING NUMBER 3? 

NOTE: 1. 40 CHARACTERS MAX LENGTH 

2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 

=> ALPHA 
= > BETA 
= > PHI 
=> THETA 
linear optimal control systems 

» USING A PRE- ALLOCATED DATASET FOR UNIT FT17F031. 

THE FCLLOWI NG OPTIONS ARE AVAILABLE. 

1. EEGIN NEW GRAPH OF OTHER CONTROLS, STATES, OR ESTIMATES, 

2. REPICT PREVIOUSLY SAVED GRAPH DATA. 
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3. EDIT THE CURRENT GRAPH. 

4. E1CI REVISED GRAPH ON THE TEK618. 

5. QUIT AND/OR MAKE METAFILE OF THE CURVES. 

PREVIOUSLY SAVED. 

SELECT A NUMBER BETWEEN 1 AND 5. 



THE GRAPH EDIT MENU 

1. CHANGE VARIABLES OR ADD A CURVE ON THE CURRENT PLOT. 

2. DELETE CURVE FRCM CURRENT PLOT. 

3. EDIT CURVE TITLE(S). 

4. EDIT PAGE HEADING (5). 

5. CHANGE THE Y-AXIS SCALE. 

6. CHANGE THE TIME AXIS SCALE. 

7. CHANGE PLOT SIZE. (DEFAULT IS 8.5 X 6.0) 

8. CHANGE THE LETTERING HEIGHT. 

9. CHANGE POSITION OF THE LEGEND. 

10. EDITING COMPLETE. 

SELECT A NUMBER EETWEEN 1 AND 10. 



CN WHICH CURVE EC YOU WANT TO CHANGE THE Y-5CALE? 
ENTER CURVE NUMBER- 1, 2, 3, OR 4 



WHAT IS THE NEW Y -MI N VALUE AT THE ORIGIN? 

•> 

-.075 

WHAT IS THE NEW Y-MAX VALUE? 

*? 

6. 1 5 

THE GRAPH EDIT MENU 

1. CHANGE VARIABLES OR ADD A CURVE ON THE CURRENT PLOT. 

2. DELETE CURVE FRCM CURRENT PLOT. 

3. EDIT CURVE TITLE (S) . 

4. EDIT PAGE HEAD ING (5). 

5. CHANGE THE Y-AXIS SCALE. 

6. CHANGE THE TIME AXIS SCALE. 

7. CHANGE PLOT SIZE. (DEFAULT IS 8.5 X 6.0) 

8. CHANGE THE LETTERING HEIGHT. 

9. CHANGE POSITION OF THE LEGEND. 

10. EDITING COMPLETE. 

SELECT A NUMBER BETWEEN 1 AND 10. 



THE FOLLOWING OPTIONS ARE AVAILABLE. 

1. EEGIN NEW GRAPH OF OTHER CONTROLS. STATES, OR ESTIMATES. 

2. EEPLCT PREVIOUSLY SAVED GRAPH DATA. 

3. EDIT THE CURRENT GRAPH. 

4. PICT REVISED GRAPH ON THE TEK618. 

5. QUIT AND/OR MAKE METAFILE OF THE CURVES. 

PREVIOUSLY SAVED. 

SELECT A NUMBER BETWEEN 1 AND 5. 
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TEE F01L0WING OPTIONS ARE AVAILABLE 



1 . 

2 . 

3. 

4. 

5. 



BEGIN NEE GRAPH OF OTHER CONTROLS. STATES. 
EEFLCT PREVIOUSLY SAVED GRAPH DATA. 

EDIT THE CURRENT GRAPH. 

PICT REVISED GRAPH ON THE TEK618. 

QUIT AND/OR HAKE METAFILE OF THE CURVES. 
PREVIOUSLY SAVED. 



OR ESTIMATES. 



SELECT A NUMBER BETWEEN 1 AND 5. 

-p 

3 



THE GRAPH EDIT MENU 

1. CHANGE VARIABLES OR ADD A CURVE ON THE CURRENT PLOT. 

2. DELETE CURVE FECM CURRENT PLOT. 

3. EDIT CURVE TITLE (S) . 

4. EDIT PAGE HEADING (5). 

5. CHANGE THE Y-AXIS SCALE. 

6. CHANGE THE TIME AXIS SCALE, 

7. CHANGE PLOT SIZE. (DEFAULT IS 8.5 X 6.0) 

8. CHANGE THE LETTERING HEIGHT. 

9. CHANGE POSITION OF THE LEGEND. 

10. EDITING COMPLETE. 

SELECT A NUMBER EETWEEN 1 AND 10. 



CN WHICH CURVE DC YOU WANT TO CHANGE THE Y-SCALE? 
ENTER CURVE NUMBER- 1, 2, 3, OR 4 



WHAT IS THE NEW Y -MI N VALUE AT THE ORIGIN? 

■? 

-.06 

WHAT IS THE NEW Y-MAX VALUE? 

-p 

6. 1 5 

THE GRAPH EDIT MENU 

1. CHANGE VARIABLES OR ADD A CURVE ON THE CURRENT PLOT. 

2. DELETE CURVE FECM CURRENT PLOT. 

3. EDIT CURVE TITLE (S). 

4. EDIT PAGE HEADING (S). 

5. CHANGE THE Y-AXIS SCALE. 

6. CHANGE THE TIME AXIS SCALE. 

7. CHANGE PLOT SIZE. (DEFAULT IS 8.5 X 6.0) 

3. CHANGE THE LETTERING HEIGHT. 

9. CHANGE POSITION OF THE LEGEND. 

10. EDITING COMPLETE. 

SELECT A NUMBER EETWEEN 1 AND 10. 



THE FOLLOWING OPTIONS ARE AVAILABLE. 

1. BEGIN NEW GRAPH OF OTHER CONTROLS, STATES, OR ESTIMATES. 

2. REFLCT PREVIOUSLY SAVED GRAPH DATA. 

3. EDIT THE CURRENT GRAPH. 

4. PICT REVISED GRAPH ON THE TEK618. 

5. QUIT AND/OR MAKE METAFILE OF THE CURVES. 
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PREVIOUSLY SAVED 



SELECT A NUMBER EETWEEN 1 AND 5. 

4 

THE FOLLOWING OPTIONS ARE AVAILABLE. 

1. EFGIN NEW GRAPH OF OTHER CONTROLS, STATES, OR ESTIMATES. 

2. F.EFLCT PREVIOUSLY SAVED GRAPH DATA. 

3. EDIT THE CURRENT GRAPH. 

4. PICT REVISED GRAPH ON THE TEK618. 

5. QUIT AND/OR MAKE METAFILE OF THE CURVES. 

PREVIOUSLY SAVED. 

SELECT A NUMBER BETWEEN 1 AND 5. 

5 

DO YOU WANT TO SAVE THE CURRENT GRAPH DATA TO 
EE USED LATER TO GENERATE A METAFILE? 

Y OR N 

NOTE: A METAFILE IS REQUIRED FOR SMOOTH VERSATEC PLOTS. 
THERE WILL EE AN OPPORTUNITY TO GENERATE A METAFILE 
JUST BEFORE EXITING THIS PROGRAM, 
filt eron 

YOUR ANSWER MUST BE "Y" OR "N". 

DO YOU WANT TO SAVE THE CURRENT GRAPH DATA TO 
EE USEE LATER TO GENERATE A METAFILE? 

Y OR N 

NOTE: A METAFILE IS REQUIRED FOR SMOOTH VERSATEC PLOTS. 
THERE WILL BE AN OPPORTUNITY TO GENERATE A METAFILE 
JUST BEFORE EXITING THIS PROGRAM. 

y 



WHAT FILE NAME 



filt eron 



EC YOU WANT THE CURVE 
(8 CHARACTERS MAX) 



DATA STORED UNDER? 



THE CURVE EATA IS BEING FILED UNDER FILTERCN DATA 
END OF DISSFLA 9.0 -- 43644 VECTORS GENERATED IN 3 PLOT FRAMES 
PROPRIETARY SOFTWARE PRODUCT OF ISSCO, SAN DIEGO, CA. 

8493 VIRTUAL STORAGE REFERENCES; 9 READS; 0 WRITES. 



THE FOLLOWING OPTIONS ARE AVAILABLE: 

1. MAKE METAFILE OF PREVIOUSLY SAVED CURVE. 

2. QUIT. 

ENTER 1 OR 2 

-> 

i 



WHAT FILE NAME IS THE DATA STORED UNDER? 

filteron 

THE CURVE DATA IS BEING LOADED FROM FILE FILTERON DATA 
» USING A PEE- ALLOCATED DATASET FOR UNIT FT18F001. 
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THE FOLLOWING OPTIONS ARE AVAILABLE: 



1. HAKE METAFILE OF PREVIOUSLY SAVED CURVE. 

2. QUIT. 

ENTER 1 OR 2 



2 

END OF DISSPLA 9.0 -- 14919 
PROPRIETARY SOFTWARE PR 
2874 VIRTUAL STORAGE REFERE 
DASD 121 DETACHED 
DASD 122 DETACHED 
DASD 120 DETACHED 



VECTORS GENERATED IN 1 PLOT FRAMES 
ODUCT OF ISSCO, SAN DIEGO, CA. 

NCES; 9 READS; 0 WRITES. 



********************************************* **#* I,:**#*******#* 

DO YOU WANT A VRSTEC PLOTTER SMOOTH COPY OF THE 
TEE DISSPLA METAFILE THAT YOU JUST CREATED? 

(Y OR N) 

************************************************* ************* 
I (120) E/C 

DASD 001 LINKED R/O ; E/W 3Y MVS 
Z (001) E/C - OS 
DASD 001 DETACHED 
CREATING NEW FILE: 
rPFATTNG NPU HTTP* 

PUN FILS 6910 TO* MVS COPY 001 NOHOLD 

DASD 120 DETACHED 



*** * * ********* ** ** ********* ****************** **** ************* 
YOUR GRAPH (S) CAN BE PICKED UP AT THE COMPUTER CENTER. 

THE GRAPH (S) WILL BE ADDRESSED TO "POP (USER ID)''. 
************************************************************** 



*********** *** ** ** ****** *** **************** ** ********** ** *** * * 
DO YOU WANT TO 

1. RUN OPTSYSX AGAIN 

2. RUN THE PLOT PROGRAM USING THE SAME MATRICES? 

(TO PLOT ANOTHER TYPE OF SYSTEM (OPEN /CLOSED) ) 

3. QUIT 

ENTER 1, 2 OR 3 

************************************************************** 

3 

*********** *** * * ** * * * * ***** ****** ************ *** ******* ** ** * * * 



HAVE A GOOD DAY! ! 



*****************************************£*** ***************** 
P; 7=13.37/21.24 21:23:09 
record off 

END EECCFDING OF TERMINAL SESSION 
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output generated by this 



example follows 



The graphical 
as figure 3.3. 



8 1 



FILTER ONLY CLOSED LOOP 
EXAMPLE 4.1 

LINEAR OPTIMAL CONTROL SYSTEMS 




CJ 

u 

U1 



w 



09i*o 92ro ooro 



9 £0 ' 0 090*0 9 2 0*0 000*0 920*0- 090*0*9£0*0- 

T* - NOIlISOd ilYlflOMV 



I 1 — 

091*0 92T0 



001*0 9£0 ' D 090*0 920*0 000*0 920*0- 09 0* 0 *9 £0 *0 - 

T3$ - 31VNI1S3 NOIlISOd HVIOONV 
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Figure 3.3 Filter Closed-loop Tine Response 



D 



FILTER AND BEGDLATOB CLOSED LOOP TIME RESPONSE 



The following filter and regulator example was taken 
from [Eef. 7 pp. 382 - 384]. 

The full terminal session is recorded below, with user 
input at the left margin in lower case le tters or n u mber s 
below each "?". 



record on 

BEGIN RECORDING OF TERMINAL SESSION 

R; T = 0. 0 1/0 . 02 09: 08 :30 

optsys 



THE OPTSYS EXEC CONTR 

1. OPTS YS X FORTRAN 

2. OPTCALC FORTRAN 

3. OPTPLOT FORTRAN 

EACH PROGRAM PASSES INF 
THROUGH A DATA FILE WRI 
CASE, THESE FILES ARE " 
THE SIZE OF THESE FILES 
CAN USE ABOUT 207 OF TH 
ENSURE THAT SUFFICIENT 

- TYPE "E" TO EXIT, 



OLS A TRIO OF PROGRAMS : 

(SYSTEM ANALYSIS) 

CALCULATE TIME RESPONSE) 
(DISSPLA PLOTTING ROUTINE) 

ORMATION TO THE NEXT PROGRAM 
TTEN TO THE USERS DISK. IN THIS 
OPTMAT DATA” AND ''OPTPLOT DATA". 

VARY WITH THE SYSTEM ORDER, AND 
E USERS DISK SPACE. THEREFORE 
DISK SPACE IS AVAILABLE. 

ANY OTHER ENTRY TO CONTINUE - 



***** ************************************************** ******* 



YOU HAVE A DATA FILE NAMED 'OPTMAT DATA' ON YOUR A 
DISK THAT WAS PREVIOUSLY GENERATED BY THE OPTSYS 
PROGRAM AND CONTAINS THE F, G, H, GAMMA, A AND B 
MATRICES FROM THAT RUN. 

IF YOU WOULD LIKE TO USE THESE SAME MATRICES FOR 
THIS RUM, THE CPTSYS PROGRAM WILL READ IN THE 
DESIRED DATA AT THE APPROPRIATE TIME, 

IF YOU TYPE (Y) ES . 

ANY OTHER INPUT WILL RESULT IN THAT FILE BEING ERASED! 

************************************************************** 

y 



*********** ******* ******************************************** 

DO YOU WANT THE NUMERICAL OUTPUT FROM OPTSYSX TO GO 
TO YOUR TERMINAL S ( CREEN) OR TO A D (ISK) FILE? 

(S OR D) 

**************** ********************************************** 
S 



************************************************************** 
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OUTPUT Will COME TO YOUR TERMINAL SCREEN. 

************ ********************************* ********** ******* 

LOADING OPT SYS 

EXECUTION EEGINS . . . 

OPTSYSX IS A COMPLETELY INTERACTIVE OPTIMAL SYSTEMS CONTROL 
FROGRAM. IT WILL SOLVE NUMEROUS CONTROL PROBLEMS ON THE 
FOLLOWING TYPES OF SYSTEMS CONTROL EQUATIONS: 

XDOT = {F}*X + {G} *U + £G AM} * ( W+K 0) 

MEASUREMENT EQUATION — 

Z = £H} *X + £D} *U + V 

REGULATOR PERFORMANCE INDEX — 

J = 1/2 * INTEGRAL (Y * {A} *Y ♦ U * £B| *U) DT 

STATE FEEDBACK GAIN DEFINITION — 

U = - {C } *x 

DO YOU WISH TO CONTINUE? TYPE "YES" OR "NO". 

Y —DATA ENTRY — 

ALTHOUGH OPTSYSX IS SPECIFICALLY DESIGNED TO READ 
ALL MATRIX DATA INTERACTIVELY, SEVERAL ALTERNATE 
MEIHODS ARE AVAILABLE TO USERS: 

METHOD 1 — T HE "F","G", AND "GAMMA" MATRICES 
MAY BE READ FROM SEPARATE DATA FILES. 

METHOD 2 — THE "F","G", AND "GAMMA" MATRICES MAY BE 
EXPLICITLY DEFINED WITHIN SUBROUTINE "SETUP". 

{NOTE: IN EITHER CASE, THE USER SHOULD OBTAIN A COPY 

OF THE PROGRAM LISTING AND EXAMINE 
THE EXAMPLES CONTAINED IN S/R "SETUP".} 

DO YOU WISH TO CONTINUE? TYPE "YES" OR "NO". 

y 



DO YOU WISH TO INPUT THE "F", "G", AND "GAMMA" 
MATRICES FRCM SUBROUTINE "SETUP" IAW THE 
METHOD DESCRIBED ON THE PREVIOUS SCREEN? 

TYPE "YES" OR "NO". 

GENERAL OPTSYSX OPTIONS: 

OETICN 1 — SYSTEM ANALYSIS WITHOUT 

OPEN-LOOP EIGENS YSTEM CALCULATIONS. 



OPTION 2 — SYSTEM ANALYSIS WITH OPEN- LOOP 
EIGENSYS TEM CALCULATIONS. 

OPTION 3 — OPEN- LOOP EIGENS YS TEM FOUND 
AND PROGRAM TERMINATES. 

£»F"-MATRIX ENTRY FOLLOWS IMMEDIATELY.} 

OPTION 4 — MCDAL DISTRIBUTION MATRICES COMPUTED 
WITHOUT FILTER OR REGULATOR SYNTHESIS 
OR STEADY-STATE ANALYSIS. 
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-Cr *0 



7 



SELECT AN OPTION: 1,2,3, OR 4 



DO YOU DESIRE RMS VALUES OF STATE AND CONTROL? 
TYPE "YES" OR ’’NO”. 



CPTSYSX LQR/CLA SSI CAL OPTIONS: 



OETION 


1 


— OPTIMAL FILTER 
SYNTHESIS WITH 
MATRIX INPUT. 


AND/OR REGULATOR 
NO EXTERNAL "C" OR "K" 


OETION 


2 


— OPTIMAL FILTER 
SYNTHESIS KITH 
MATRIX INPUT. 


AND/OR REGULATOR 
EXTERNAL "C" 


OETION 


3 


— OPTIMAL FILTER 
SYNTHESIS WITH 
MATRIX INPUT. 


AND/OR REGULATOR 
EXTERNAL »K" 


OPTION 


4 


— OPTIMAL FILTER 
SYNTHESIS WITH 
MATRIX INPUT. 


AND/OR REGULATOR 
EXTERNAL "C" AND "K" 


SELECT 


AN 


OPTION: 1, 2, 3, 


OR 4. 



DO YOU NISH TO DETERMINE THE STEADY-STATE RESPONSE 
FOR A CONSTANT DISTURBANCE? 

TYPE '’YES" CR "NO". 

DO YOU WISH TO DETERMINE THE MODAL DISTRIBUTION 
AND GAIN MATRICES? 

TYPE "YES" CR "NO". 



OPEN-LCCP TRANSFER FUNCTION OPTIONS: 

OETION 1 — NO OPEN-LOOP TRANSFER FUNCTIONS COMPUTED. 

CETION 2 — POLES, RESIDUES, AND ZEROS COMPUTED. 
OPTION 3 — ONLY POLES AND ZEROS COMPUTED. 

CETION 4 — ONLY POLES AND RESIDUES COMPUTED. 
SELECT AN OETION: 1, 2, 3, OR 4. 



NOISE TRANSFER FUNCTION OPTIONS: 

OETION 1 -- NO NOISE TRANSFER FUNCTIONS COMPUTED. 
OETION 2 — PCXES, RESIDUES, AND ZEROS COMPUTED. 
OETION 3 — ONLY POLES AND ZEROS COMPUTED. 

OETION 4 — ONIY POLES AND RESIDUES COMPUTED. 
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SELECT AN OPTICN: 1, 2, 3 , OR 4 



COMPENSATOR TRANSFER FONCTION OPTIONS: 



-> 

i 



OPTION 1 — 
OPTION 2 — 
OPTION 3 — 
OPTION 4 — 
{NOTE: 



NO COMP. TRANSFER FUNCTIONS CONFUTED. 

PCIES, RESIDUES, AND ZEROS COMPUTED. 

ONLY POLES AND ZEROS COMPUTED. 

ONLY POLES AND RESIDUES COMPUTED. 

A COMPENSATOR TRANSFER FUNCTION CAN EE 
COMPUTED ONLY IF BOTH A REGULATOR 
AND FILTER ARE SYNTHESIZED 
ANE/OR INPUT.} 



SELECT AN OPTION: 1, 2, 3, OR 4. 



WILL A FEED-FORWARD DISTRIBUTION MATRIX 
{"D" - MATRIX} EE INPUT ? 

TYPE "YES" OR "NO". 



THIS OPTICN DETERMINES THE CRITERIA FOR DECIDING WHEN A 
MARKOV PARAMETER IS ZERO-THE MARKOV PARAMETER INDICATES 
THE ORDER OF THE NUMERATOR POLYNOMIAL OF EACH TRANSFER 
FUNCTION. 

ALL »N" ZEROS OF THIS POLYNOMIAL ARE PRINTED OUT AND 
THIS TEST TELLS HOW MANY EXTRA ROOTS EXIST AT Z = 0. 
LESS THAN 10.0** £-11} IS CONSIDERED ZERO. 



THE DEFAULT VALUE OF THIS PARAMETER {IE} IS 6. 

IN CTHER WORDS, IE = 1.0E-6. 

IF YOU DESIRE A DIFFERENT MARKOV CRITERIA, 

TYPE THE INTEGER VALUE. 

IF YOU DESIRE THE DEFAULT VALUE, TYPE "0" {ZERO} 



DO YOU DESIRE TO SYNTHESIZE A STA3LE FILTER {OR REGULATOR} BY 
DESTABILIZING THE ORIGINAL SYSTEM? 

{NCTEtWORKS FCE FILTER OR REGULATOR BUT NOT FOR BOTH 
IN THE SAME RUN.} 

TYPE "YES" OR "NO". 

DC YOU DESIRE TO PRINT THE EULER-LA GRANGE El GEN SYSTEM 
PRIOR TO DECOMPOSITION {FOR CHECKING THE PROGRAM}? 

TYPE "YES" CR "NO". 



POWER SPECTRAL DENSITY {PSD} OPTION 1 : 

OPTION 1 — COMPUTE THE PSD OF THE OUTPUTS AND/OR THE 
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CONTROLS OF THE CONTROLLED SYSTEM WHEN FORCED BY 
PROCESS AND MEASUREMENT NOISE. f NOTE : DOTH A 

REGULATOR AND A FILTER MUST BE RESIDENT IN THE 
PROGRAM TC USE THIS OPTION.} 

OPTION 2 — SAME AS OPTION 1 A30VE BUT ONLY PRINT THE 
RESIDUES CF EACH TRANSFER FUNCTION 
USED IN THE PSD COMPUTATION. 

CPTICN 3 — NOT DESIRED. 



SELECT AN OPTION: 1, 2, OR 3. 

*> 

3 



DO YOU DESIRE REGULATOR SYNTHESIS ONLY? 
TYPE "YES" CR "NO". 



n 



THE "F", "G", "H" , "GAM", "A" AND " E" MATRICES 

FROM YOUR PREVIOUS OPTSYS RUN WERE SAVED. 

THE FOLLOWING OPTIONS ARE AVAILABLE: 

1. USE ALL OF THE SAME MATRICES AGAIN. 

2. USE SELECTED MATRICES AGAIN. 

3. INPUT ALL NEW MATRICES. 

ENTER 1, 2 , OR 3. 

NOTE; EACH SAVED MATRIX WILL BE REDISPLAYED AT 
THE PROPER INPUT SEQUENCE INTERVAL 
AND YOU Will HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 



IOL 

0 

ISET 

0 



F1AG/PA RAMET ER SETTINGS FOR THIS RUN ARE AS FOLLOWS: 
IQ IR ISS IM ITF1 ITF2 ITF3 IFDFW IE IDEBUG 



0 3 


0 


0 


0 


0 


0 


0 


0 


0 


IDSTAE 


IP SD 


I YU 


INORM 


I REG 


NS 


NC 


NOB 


NG 


0 


0 


0 


0 


0 


2 


1 


1 


1 



ORDER OF S 
NUMBER CF 
NUMBER CF 
NUMBER CF 



YSTEM = 2 

CONTROLS = 1 

OBSERVATIONS = 1 

PROCESS NOISE SOURCES 



THE SYSTEM MATRIX [" F" -MATRIX} .. . 

0.0 1.00000 

0.0 -4.60C00 



DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 



B7 



n 



TYPE "YES" OR "NO" 



OPEN LOOP DYNAMICS MATRIX F.. 

0. 0 0.100 0D+0 1 

0. 0 -0.46 00E+01 

THE MEASUREMENT SCALING MATRIX { "H"-MAT RIX} . . . 

1.00000 0.0 

DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 

TYPE "YES" CR "NO", 
n 

MEASUREMENT SCALING MATRIX H. . 

0. 10C0D+01 0.0 

THE CONTROL DISTRIBUTION MATRIX ["G "-MATRIX} .. . 

0.0 

0.78700 

DC YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 

TYPE "YES" CR "NO", 
n 

ENTER THE FEEDBACK GAIN CONTROL MATRIX ("C"- MATRIX} . 

DIMENSION = # CONTROLS £NC} X # STATES {NS}. 

THE ELEMENT C ( 1, 1) = 

-254. 1 

THE ELEMENT C ( 1, 2) = 

-19.57 

THE FEEDBACK GAIN CONTROL MATRIX {"C"-MATRIX} 
-254. ICO 0C - 19.57000 

DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 

TYPE "YES" OR "NO", 
n 

THE CONTROL DISTRIBUTION MATRIX G.. 

0.0 

0. 7870D+00 

THE OPTIMAL FEEDBACK GAIN CONTROL MATRIX. . .C=BINV*GT*S.. . 
-2.5410D+02 -1.9570D+01 
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THE CLOSED LOOP DYNAMICS MATRIX 



F-G*C 



0.0 1.000000D+00 

-1. 999767D+02 -2 .000 159D+01 

C-LOCP SOBOPT. REG. E- VALUES. .. DE T ( SI- F+G *C) . . 

-1. 000C8D+01, 9. 99 8 0 4E+00 : 

C-LOCP FIGHT EIGENVECTOR MATRIX M 

-5. 0 00980D-02 -4 . 9 9 9 602D-02 
1.000000D+C0 0.0 

CCNTFCL EIGENVECTCR MATRIX C*M.. 

-6.86251 0D+ 0 0 1.270399D+01 

C-ICCP SUBOPT-REG. LEFT E-VECTOR MATRIX. . M-IN V 

0. 0 1.000 OOOD + OO 

-2 . 0 0 015 9D+ 01 -1 .000 276D+00 

THE PROCESS NOISE DISTRIBUTION MATRIX 
{"GAM M A" -MATRIX] 

0.0 

0.10000 

DO YOU KISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 

TYPE "YES" CR "NO", 
n 

ENTER THE PROCESS NOISE PSD WEIGHTING MATRIX 
{"C”MATRIX} . 

DIMENSION = # PROCESS NOISE SOURCES [NG] X 
# PROCESS NOISE SOURCES {NG} 

THE ELEMENT Q ( 1, 1) = 

io 

THE PROCESS NOISE WEIGHTING MATRIX Q.. 

1 0. 00000 

DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 

TYPE "YES" CR "NO". 

I) 

PROCESS NOISE DISTRIBUTION MATRIX GAMMA.. 

0.0 

0. 10C0D+00 
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POWER SPECTRAL DENSITY - PROCESS NOISE Q. . 

0. 1 0C0D+02 

ENTER THE MEASUREMENT NOISE DISTRIBUTION MATRIX { "R "MAT RIX} . 

DIMENSION = # OBSERVATIONS {NO} X # OBSERVATIONS {NO} 
THE ELEMENT R( 1, 1)= 

6. ooooooi 

THE MEASUREMENT NOISE DISTRIBUTION MATRIX R... 

0 . 00000 

DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 

TYPE "YES" OR "NO", 
n 

POWER SPECTRAL DENSITY- MEASUREMENT NOISE.. R. . 

0. 1 0C0D-06 

ENTER THE FEEDBACK GAIN ESTIMATOR MATRIX {"X "-MATRIX} . 
DIMENSION = # STATES {NS} X # OBSERVATIONS {NO}. 



-5 


THE 


ELEMENT 


K( 


1, 


1) = 


95. 4 


~> 

456 1 


THE 


ELEMENT 


K( 


2, 


1) = 



THE FEEDBACK GAIK ESTIMATOR MATRIX {"K"-MATR IX} 

95.40000 

4561.00000 

DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 

TYPE "YES" CR "NO", 
n 

FILTER STEADY STATE GAINS K... 

9.5400C0D+01 
4.56 10 C0D+ 03 

THE CLOSED LOOP FILTER DYNAMICS MATRIX IS.... 

-9. 540000D+C1 1 .0 00000D+00 

-4.561 0 00D+03 -4 . 6 00 OOOD+OO 

EIGENSYSTEM OF OPTIMAL ESTIMATOR 

C-IOCP SU30PT. EST. E- VALUES.. .DET(SI-F+K*H). . 
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-5. 00 0C0D+0 1, 4. 99984D+01: 



C-ICOP EIGHT EIGENVECTOR MATRIX M 

9. 9 5 3957D- 03 - 1 . 0 962 16D-02 
1.00000CD+00 0.0 

MEASUREMENT EIGENVECTOR MATRIX H (EAR) *M. . 

9.953957D-03 - 1 . 096216D-02 

C-ICCP SUBOPT. FILTER LEFT E-VECTOR MATRIX. .» -IN V . . 

0.0 1.0000C0D+00 

-9. 122292D+C1 9.080291D-01 



THE COVARIANCE OF THE ESTIMATION ERROR P.. 

7. 1 5 05 03D-06 2.271 000D-04 

2. 27 1000D- C4 1.18 115 ID-02 



RMS VALUES 
2. 6 7 404 2D- 03 



OF THE ESTIMATION ERROR 
1 . 0 86 8C8D-0 1 



DO YOU WISH TO OBTAIN A TIME RESPONSE 
OF THE SYSTEM YOU ARE EVALUATING? 

(Y OR N) 

NOTE: YOU MUST BE LOGGED ON AT A DUAL SCREEN 
(TEK 6 18) TERMINAL TO UTILIZE THIS MODE. 

THE F (SYSTEM) , G (CONTROL). H (OBS ERV AE LES ) , 3AM (NOISE) , 
A (OUTPUT COST) AND B (CONTROL COST) MATRICES KILL BE 
SAVED FOR REENTRY TO THE MAIN OPTSYS PROGRAM. 

Y 



*********** *** ** ** * ******** **************** ** ** ** ***** *** **** * 
IF YOU ARE DISSATISFIED WITH THE RESULTS 
THUS FAR AND WCULD LIKE TO EXIT TO CMS / 

-TYPE «Y' TO EXIT- 

(ANY OTHER INPUT TO CONTINUE) 

**************** ** *************************************** *** * * 

LOADING CPT CALC. ... 

EXECUTION BEGINS.. . 



DURING THIS SECTION OF THE PROGRAM YOU KILL: 

- SELECT THE TYPE OF SYSTEM RESPONSE TO PLOT 

( CPE N LCOP, CLOSED LOOP, OR FILTER/REGULATOR) 

- PROVIDE START AND STOP TIME FOR PLOTTING CALCULATIONS 

- SELECT THE TYPE CF DRIVING FUNCTION (S) (STEP OR RAMP) 

- PROVIDE START AND STOP TIMES FOR THE DRIVING FUNCTION (S) 

- PROVIDE DRIVING FUNCTION MAGNITUDE (S) . 
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CLEAR THE SCREEN TO CONTINUE 



THE E MATRIX 



0.0 
0 . 0 



1 .00000 

-4.60000 



THE G MATRIX 



0 . 0 

0.78700 



THE C MATRI X 



-254.10000 -19.57000 



THE H MATRI X 



1.00000 



0.0 



THE K MATRIX 



95.40000 
4561. 00000 



TEE FOLLOWING PLOTTING OPTIONS ARE AVAILABLE IF THE 
REQUIRED MATRICES WERE CALCULATED IN GPTSYSX: 

1. OPEN LOOP TIME RESPONSE 

XDOT = {F} *X + {G} * UC 

2. CICSED LOOP TIME RESPONSE 

XDOT = {F-G*C} *X + {G} *UC, U = {C} *X 

3. OPTIMIZED FILTER CLOSED LOOP SYSTEM RESPONSE. 



4. FILTER + REGULATOR CLOSED LOOP SYSTEM RESPONSE. 

XDOT = f F+G*C} *X + {G} *UC. Z = {H}*X 

XHDCT = {F} * XH + {G} U + {K}*{Z - H*XH} , U = {C)*XH 

SELECT 1 , 2, 3 OK 4. 




4 



THE (G*C) MATRIX 



0 . 0 

-199.97670 



0.0 

- 15.40159 



THE (K*H) MATRIX 

95.40000 0.0 

4561. 00000 0.0 



THE COMBINED SYSTEM F MATRIX (2*NS X 2* NS) 



0.0 
0 . 0 



1.00000 0.0 

-4.60000 -199.97670 



0.0 

-15.40159 
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0*0 



S 5. 40000 0.0 - 95.40000 1.00000 

4561.00000 0.0 - 4760.97670 - 20.00159 



THE AUGMENTED G MATRIX (2*NS X NC) 

0 . 0 

0.78700 

C. 0 

0 . 78700 



AT WHAT TIME DO YOU WANT TO START 
THE TIME RESPONSE CAICUIATIONS? 

INPUT START TIME IN SECONDS. (NORMALLY 0.0) 



AT WHAT TIME DO YOU WANT TO STOP 
TEE TIME RESPONSE CALCULATIONS? 

INPUT STOP TIME IN SECONDS. 

0.6 



THIS PROGRAM DIVIDES THE TIME INTERVAL YOU HAVE 
JUST SPECIFIED INTO UP TO 500 SMALL INTERVALS FOR 
THE INTEGRATION AND PLOTTING ROUTINES. IN ORDER 
10 SAVE COMPUTER TIME. THE NUMBER OF POINTS CAN BE 
CAN BE REDUCED WITH SOME LOSS IN CURVE FIDELITY. 

HCW MANY POINTS DO YOU WANT TO CALCULATE? 

500 

DOES THE SYSTEM UTILIZE A DRIVING FUNCTION (CONTROL INPUT)? 

( Y ) ES OR (N) 0 
n 



DOES THE SYSTEM START WITH ALL INITIAL CONDITIONS = 0.0 ? 



(Y) ES OR (N) 0? 
n 





WHAT 


IS 


THE 


INITIAL 


CONDITION 


FOR 


x ( 


D 


7 




6 *i 






















7 

n 


WHAT 


IS 


THE 


INITI AL 


CONDITION 


FOR 


XHAT ( 


D 


7 


u 


WHAT 


IS 


THE 


INITIAL 


CONDITION 


FOR 


x ( 


2 ) 


7 




6 
























WHAT 


IS 


THE 


INITI AL 


CONDITION 


FOR 


XHAT ( 


2 ) 


7 


6 
























THIS 

MAKE 


IS YOUR 
CHANGES 


LAST OPPORTUNITY 
IN THE FOLLOWING 


TO 

AR EAS . 









1. SELECT ANOTHER TYPE OF SYSTEM TO PLOT 



93 



(OPEN, CLOSED, FILTER OR FILTER/REGULATOR) 

2. STAR 1 AND STOP TIKES 

3. DRIVING FUNCTIONS 

4. INITIAL CONDITIONS 

5. CONTINUE 

SELECT A NUMEER BETWEEN 1 AND 5. 

THE FOLLOWING INFORMATION IS PROVIDED ONLY 
FOR AN INDICATION OF PROPER PROGRAM OPERATION. 

ALL CONTROLS, STATES AND STATE ESTIMATES CAN BE PLOTTED. 

TIME U(1) X(1) X (2) X (3) 



0.0 


0. 0 


0. 1 000000D+00 


0. 01 


0. 0 


0. 9816154D-01 


0.02 


0. 0 


0. 8962622D-01 


0. 04 


0. 0 


0.7 5 42258D - 0 1 


0.05 


0. 0 


0.5 899229D-0 1 


0.06 


0. 0 


0.451 0993D -0 1 


0.07 


0. 0 


0 . 29204 10D-0 1 


0.08 


0. 0 


0. 1767465D-01 


0. 10 


0. 0 


0.8 388221 D -02 


0. 1 1 


0. 0 


0. 1 034654D-02 


0. 12 


0. 0 


-0. 4696796D-02 


0. 13 


0. 0 


-0. 9C63199D-02 


0. 14 


0. 0 


-0. 1 227031 D-0 1 


0. 16 


0. 0 


-0. 1448987D-01 


0. 17 


0. 0 


-0. 1587379D-01 


0. 18 


0. 0 


-0. 16561 21D-01 


0. 1 9 


0. 0 


-0. 1668037D-01 


0. 20 


0. 0 


-0. 1634808D-01 


0.22 


0. 0 


-0. 1 5668 75D -0 1 


0.23 


0. 0 


-0. 1 473380D-0 1 


0.24 


0. 0 


-0. 1 3621 84D-0 1 


0.25 


0. 0 


-0. 1239918D-01 


0.26 


0. 0 


-0. 1 112072D-01 


0. 28 


0. 0 


-0. 9831066D-02 


0.29 


0. 0 


-0. 8565537D-02 


0. 30 


0. 0 


-0.7 351302D-02 


0.31 


0. 0 


-0. 6208421D-02 


0.32 


0. 0 


-0. 5150854D-02 


0.34 


0. 0 


-0.4 187403D-02 


0. 35 


0. 0 


-0 . 3 222593D - 02 


0.36 


0. 0 


-0. 2557462D-02 


0.37 


0. 0 


-0. 1890282D-02 


0.38 


0. 0 


-0. 1317192D-02 


0.40 


0. 0 


-0.8 227506D-03 


0.41 


0. 0 


-0. 4304142D-03 


0.42 


0. 0 


-0. 1 029384D-03 


0.43 


0. 0 


0. 1 572873D-03 


0. 44 


0. 0 


0. 3579670D-03 


0.46 


0 . 0 


0. 50668 18D -03 


0.47 


0 . 0 


0.6 107220D-03 


0.48 


0 . 0 


0. 6769634D-03 


0.49 


0 . 0 


0.7 1 17801D-03 


0.50 


0 . 0 


0. 7209884D-03 


0. 52 


0 . 0 


0.7 098 1 76D-03 


0. 53 


0 . 0 


0.6 8290 17D-03 


0.54 


0 . 0 


0.6 442879D-03 


0.55 


0 . 0 


0. 5S74586D-03 


0. 56 


0 . 0 


0. 5453629D-03 



0.0 0.0 
•0 . 4069 147D + 00 0. 8 100300D-0 1 
0. 9917358D+00 0. 1 0 420 1 4D+0 0 
0. 1323812D+01 0. 93 794 55D-0 1 
0. 1375381D+01 0. 7 1 245 87D-0 1 
0. 1252435D+01 0. 48677 1 5D-0 1 
0. 1060393D+01 0. 30 55684D-0 1 
0. 8635309D+00 0. 17 22332D-0 1 
0. 6887453D+00 0. 75 7 1 9 09D-0 2 
0.5412268D+00 0. 431 079SD-03 
0.4175739D+00-0. 4997717D-02 
0. 3129904 D+ 0 0- 0.91 54016D-0 2 
0 .2238 93 0D+ 00-0. 12261 95D-01 
0.1481 34 6D+00-0. 14 45 494D-0 1 
0. 84 44 92 9D-0 1-0. 15 84 4 67D-01 
0.31 90280 D- 01-0. 1654535D -01 
0. 10 41 88 0D-01-0. 16 674 77D-01 
0.4351 168D-01-0. 1634779D-01 
0. 6844 1 37 D- 01-0. 15 67 017D-01 
0. 8629 163D-0 1-0. 1473518D-01 
0. 9811 538D- 01-0. 1 3622 66D-01 
0. 1049 01 5D + 0 0-0. 1 239950D-01 
0.1075564D+00-0. 1 1 12077D-01 
0. 10 6897 0D+ 0 0-0. 98310 12D-02 
0. 10 36 48 6D+ 00-0. 85 6 54 74D-0 2 
0. 9844669D-01-0.7351261D-02 
0. 9184092D-01-0.6208403D-02 
0. 8429 932D- 01 -0.5 150849 D-0 2 
0. 7621404D-01-0.41874C5D-02 
0. 6790754 D- 01-0. 3322596D -02 
0. 5963914 D- 01-0. 2557464D-02 
0.5161 190D-01-0. 1390283D-02 
0. 4397 938D- 01-0. 1317192D-02 
0. 3685230D-01-0.8327506D-03 
0. 30 3048 4D- 01- 0.43041 41 D-0 3 
0. 2438 053D- 01-0. 1029383D-03 
0. 1909772D-01 0. 1 5 72 873D-03 
0. 1445443D-01 0. 3579670D-03 
0. 1043279D-01 0. 50 66 8 1 8D-0 3 
0.70 0286 2D- 02 0. 6 1 07 220D-0 3 
0.4125 9 33D- 02 0. 676 9634D-03 
0. 1757615D-02 0. 7 1 1 78 01 D-0 3 
0. 1502632D-03 0. 7 20 9884D-03 
0. 1647187D-02 0. 70981 76D-0 3 
0. 2782618D-02 0. 68 29 0 17D-0 3 
0. 36 04 6 93D- 02 0. 64 42879D-03 
0.41 59284D- 02 0. 59 74586D-03 
0 . 44 89 3 04D- 02 0. 54 536 29 D-0 3 
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0.58 0.0 0. 4904551D-03-0.4634245D-02 0. 4904551D-03 

0. 5 S 0. 0 0.4347377D-03-0.4629897D-02 0. 4 34 7377 D-0 3 

0.60 0.0 0. 3 798064D-03-0. 4508230 D- 02 0. 3 7980 64D-0 3 



** * * *********** *** ********* ****************** **** ♦** + ** 4 * ** * * * 



IF YOU ARE DISSATISFIED WITH THE RESULTS 
THUS FAR AND WOULD LIKE TO EXIT TO CMS, 

-TYPE 1 Y * TO EXIT- 

JANY OTHER INPUT TO CONTINUE) 



**************** ********************************************** 



B 

C 

E 



120 

121 

122 



R/O 

R/O 

R/O 



... Your Fortran program is now being loaded ... 
... execution will soon follow ... 
EXECUTION BEGINS... 



THIS PORTION OF THE PROGRAM PLOTS: 

_ the states 

- EXTERNAL C6NTROL INPUTS, 

- FEEIEACK CONTROL INPUTS, 

- STATE ESTIMATES AND 

- RECONSTRUCTION ERRORS 

FROM THE DATA THAT YOU JUST CALCULATED. 

TEE CAPABILITY IS ALSO AVAILABLE TO REVIEW ANY 
GRAPHS THAT YOU HAD PREVIOUSLY SAVED AS DATA 
FILES ON YOUR DISK. 

CLEAR THE SCREEN TO CONTINUE. 

THE FOLLOWING OPTIONS ARE AVAILABLE: 

1. PLOT THE DATA YOU JUST CALCULATED. 

2. PLOT A CURVE THAT YOU PREVIOUSLY SAVED. 

ENTER 1 OR 2 

•? 

i 



YCU MAY PLOT UP TO 4 SYSTEM VARIABLES VS TIME. 
HOW MANY VARIABLES DO YOU WISH TO PLOT? 

3 



WHICH TYPE CF VARIABLE DO YOU WISH TO PLOT AS CURVE NUMEEK 1? 



•5 

i 



1. STATE VARIABLE (IE., Xl r X2, ETC) 

2. FEEDBACK CCNTRCL (IE. , U = -C*X) 

3. CCNTRCL INPUT (IE. . U 1, U2 , ETC.) 

4. STATE ESTIMATE (OBSERVER) (IE., XHAT1, XHAT2 , 

5. STATE RECONSTRUCTION ERROR (IE., Xl-XHATl, 
X2-XHAT2, ETC) 

ENTER 1,2, 3,4 OR 5 



ETC.) 
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i 



WHAT IS THE SUBSCRIPT OF THE STATE VARIABLE THAT 
YOU WANT TC PLOT AS THE NUMBER 1 CURVE VS TIME? 



WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 



NOTE: 



1. 40 CHARACTERS MAX LENGTH 

2. GREEK SYMBOLS WILL BE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 



IE. 



angular position - 



A) => ALPHA 

E) => BETA 

F) => PHI 

Q) => THETA 
X 1 



WHICH TYPE CF V ARI ABIE DO YOU WISH TO PLOT AS CURVE NUMBER 2? 

1. STATE VARIABLE (IE. , XI, X2, ETC) 

2. FEEDBACK CONTROL (IE., U = -C*X) 

3. CCNTRCL INPUT (IE., Ul, U2, ETC.) 

4. STATE ESTIMATE (OBSERVER) (IE., XH AT 1 , XHAI 2 , ETC.) 

5. STATE RECONSTRUCTION ERROR (IE., X1-XHAT1, 

X2-XHAT2, ETC) 

ENTER 1,2, 3,4 OR 5 

4 



i 



WHAT IS THE SUBSCRIPT OF THE STATE ESTIMATE THAT 
YCU WANT TC PLOT AS THE NUMBER 2 CURVE VS TIME? 



WHAT IS THE CURVE LABEL FOR THIS VARIABLE? 

NOTE; 1. 40 CHARACTERS MAX LENGTH 

2. GREEK SYMBOLS WILL EE PRINTED FOR ANY LETTERS 
ENCLOSED IN PARENTHESES. 

IE. (A) => ALPHA 
'B) => BETA 
F) => PHI 
lQ\ => THETA 
angular position estimate - (x)el 



WHICH TYPE CF VARIABLE DO YOU WISH TO PLOT AS CURVE NUMBER 3? 



7 

2 



1 . 
2 . 

3. 

4. 

5. 



STATE VARIABLE (IE., 
FEEDBACK CCNTRCL (IE 



XI, X2, ETC) 
v ^. , U = -C*X) 

CCNTRCL INPUT (IE., Ul, U2, ETC.) 

STATE ESTIMATE (OBSERVER) (IE., XH AT 1 , XHAI 2 , ETC.) 
STATE RECONSTRUCTION ERROR (IE., X1-XHAT1, 

X2-XHAT2, ETC) 

ENTER 1,2, 3,4 OR 5 



WHAT IS THE SUBSCRIPT OF THE FEEDBACK CONTROL THAT 
YOU WANT TO PLOT AS THE NUMBER 3 CURVE VS TIME? 
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1 



WHAT IS THE CURVE LABEL 



NOTE 



1. 

2 . 



40 CHARACTERS 
GREEK SYMBOLS 
ENCLOSED IN P 



IE. 



input voltage - v 



= > AL 
=> EE 
= > PH 
=> TH 



3 



YOU MAY USE UP TO 
HOW MANY HEADINGS 

0 , 1 , 2 



WHAT IS THE DESIRED HEA 



NOTE: 1. 40 CHARACTERS 

2. GREEK SYMBOLS 
ENCLOSED IN P 
IE. (A) => AL 
= > BE 
=> PH 
=> TH 

filter + regulator closed 1 
WHAT IS THE DESIRED HEA 



NOTE: 



1 . 

2 . 



40 CHARACTERS 
GREEK SYMBOLS 
ENCLOSED IN P 



IE. 



example 5. 1 



=> AL 
=> BE 
= > PH 
=> TH 



WHAT IS THE DESIRED HEA 

NOTE: 1. 40 CHARACTERS 

2. GREEK SYMBOLS 
ENCLOSED IN P 
IE. (A) => AL 
B) => BE 
Fj => PH 
,Q) => TH 
control syst 



linear optimal 
» USING A PRE -ALLOCATED "DA 



THE FOLLOWING OPTIONS 

1. BEGIN NEW GRAPH OF CT HER 

2. REPLOT PREVIOUSLY SAVED 

3. EDIT THE CURRENT GRAPH. 

4. PLOT REVISED GRAPH ON TH 

5. QUIT AND/OR MAKE METAFIL 
PREVIOUSLY SAVED. 



SELECT 

3 



A NUMBER EETWEEN 1 A 



FOR THIS VARIABLE? 

MAX LENGTH 

WILL BE PRINTED FOR ANY LETTERS 
ARENTHESES. 

PHA 

TA 

I 

ETA 



3 HEADINGS. 

DO YOU DESIRE ON THIS GRAPH? 
OR 3 



DING NUMBER 1? 

MAX LENGTH 

WILL BE PRINTED FOR ANY LETTERS 
ARENTHESES. 

PHA 

TA 

I 

ETA 

oop system 
DING NUMBER 2? 

MAX LENGTH 

WILL BE PRINTED FOE ANY LETTERS 
ARENTHESES. 

PHA 

TA 

I 

ETA 



DING NUMBER 3? 

MAX LENGTH 

WILL BE PRINTED FOR ANY LETTERS 
ARENTHESES. 

PHA 

TA 

I 

ETA 

ems 

TASET FOR UNIT FT17F001. 

ARE AVAILABLE. 

CONTROLS. STATES, OR ESTIMATES. 
GRAPH DATA. 

E TEK618. 

E OF THE CURVES. 



ND 5. 
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N )*0 



THE GRAPH EDIT MENU 



1. CHANGE VARIABLES OR ADD A CURVE ON THE CURRENT PLOT. 

2. DELETE CURVE FROM CURRENT PLOT. 

3. EDIT CURVE TITLE (S) . 

4. EDIT PAGE HEADING (S) . 

5. CHANGE THE Y-AXIS SCALE. 

6. CHANGE THE TIME AXIS SCALE. 

7. CHANGE PLOT SIZE. (DEFAULT IS 3.5 X 6.0) 

8. CHANGE THE LETTERING HEIGHT. 

9. CHANGE POSITION OF THE LEGEND. 

10. EDITING COMPLETE. 

SELECT A NUMBER EETWEEN 1 AND 10. 



•? 

6 



HOW MANY 
(LEFT OR 
MOVE THE 

NOTE: 1. 

2 . 

3. 



INCHES IN THE X DIRECTION 

RIGHT), DO YOU WANT TO MOVE 

LEGEND BOX FROM ITS PRESENT POSITION 

DEFAULT PLOT SIZE IS 6.5 X 6.0 
LEFT IS NEGATIVE 
RIGHT IS POSITIVE 



HOW MANY INCHES IN THE Y DIRECTION 
(UP OR DOWN), DO YOU WANT TO MOVE 
MOVE THE LEGEND BOX FROM ITS PRESENT POSITION 

NOTE: 1. DEFAULT PAGE SIZE IS 8.5 X 6.0 

2. DOWN IS NEGATIVE 

3. UP IS POSITIVE 



THE GRAPH EDIT MEND 

1. CHANGE VARIABLES OR ADD A CURVE ON THE CURRENT PLOT. 

2. DELETE CURVE FROM CURRENT PLOT. 

3. EDIT CURVE TITLE (S). 

4. EDIT PAGE HEAEING(S) . 

5. CHANGE THE Y-AXIS SCALE. 

6. CHANGE THE TIME AXIS SCALE. 

7. CHANGE PLOT SIZE. (DEFAULT IS 8.5 X 6.0), 

8. CHANGE THE LETTERING HEIGHT. 

9. CHANGE POSITION OF THE LEGEND. 

10. EDITING COMPLETE. 

SELECT A NUMBER BETWEEN 1 AND 10. 



ON WHICH CURVE DO YOU WANT TO CHANGE THE Y- SCALE? 
ENTER CURVE NUMBER- 1, 2, 3, OR 4 



WHAT IS THE NEW Y-MIN VALUE AT THE ORIGIN? 

•> 

-.025 

WHAT IS THE NEW Y-MAX VALUE? 

6. 125 

THE GRAPH EDIT MENU 
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1. CHANGE VARIABLES OR ADD A CURVE ON THE CURRENT PLOT. 

2. DELETE CURVE FROM CURRENT PLOT. 

3. EEIT CURVE TITLE (S). 

4. EDIT PAGE HEALING (5) . 

5. CHANGE THE Y-AXIS SCALE. 

6. CHANGE THE TIME AXIS SCALE. 

7. CHANGE PLOT SIZE. (DEFAULT IS 8.5 X 6.0), 

8. CHANGE THE LETTERING HEIGHT. 

9. CHANGE POSITION OF THE LEGEND. 

10. EDITING COMPLETE. 

SELECT A NUMBER EETWEEN 1 AND 10. 

•7 

io 

TEE FOLLOWING OPTIONS ARE AVAILABLE. 

1. BEGIN NEW GRAPH OF OTHER CONTROLS. STATES, OR ESTIMATES. 

2. REPLCT PREVIOUSLY SAVED GRAPH DATA. 

3. EDIT THE CURRENT GRAPH. 

4. PLCT REVISED GRAPH ON THE TEK618. 

5. QUIT AND/OR MAKE METAFILE OF THE CURVES. 

PREVIOUSLY SAVED. 

SELECT A NUMBER EETWEEN 1 AND 5. 

*7 

4 



THE FOLLOWING OPTIONS ARE AVAILABLE. 

1. BEGIN NEW GRAPH OF OTHER CONTROLS, STATES, OR ESTIMATES. 

2. REPLCT PREVIOUSLY SAVED GRAPH DATA. 

3. EDIT THE CURRENT GRAPH. 

4. PLOT REVISED GRAPH ON THE TEK618. 

5. QUIT AND/OR MAKE METAFILE OF THE CURVES. 

PREVIOUSLY SAVED. 

SELECT A NUMBER EETWEEN 1 AND 5. 



DC YOU WANT TO SAVE THE CURRENT GRAPH DATA TO 
EE USED LATER TO GENERATE A METAFILE? 

Y OR N 

NOTE: A METAFILE IS REQUIRED FOR SMOOTH VERSATEC PLOTS. 
THERE WILL BE AN OPPORTUNITY TO GENERATE A METAFILE 
JUST BEFORE EXITING THIS PROGRAM. 



WHAT FILE NAME DC YOU WANT THE CURVE DATA STORED UNDER? 

(8 CHARACTERS MAX) 

filt ereg 

THE CURVE IATA IS BEING FILED UNDER FIL TEE EG DATA 



END OF DISSPLA 9.0 -- 26332 VECTORS GENERATED IN 2 PLOT FRAMES 
PROPRIETARY SOFTWARE PRODUCT OF ISSCO, SAN DIEGO, CA. 

5020 VIRTUAL STORAGE REFERENCES; 9 READS; 0 WRITES. 



THE FOLLOWING OPTIONS ARE AVAILABLE: 

1. MAKE METAFILE OF PREVIOUSLY SAVED CURVE. 

2. QUIT. 
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ENTER 1 OR 2 



WHAT FILE NAME IS THE DATA STORED UNDER? 

filt ereg 



THE CURVE DATA IS BEING LOADED FROM FILE FILT EREG DATA 
» USING A ERE- ALLOCATED DATASET FOR UNIT FT18F031. 



THE FOLLOWING OPTIONS ARE AVAILABLE 



1. MAKE METAFILE OF PREVIOUSLY SAVED CURVE 

2. QUIT. 



ENTER 1 OR 2 



2 

END CE DISSELA 9.0 -- 13201 VECTORS GENERATED IN 1 PLOT FRAMES 
PROPRIETARY SOFTWARE PRODUCT OF ISSCO, SAN DIEGO, CA. 

2772 VIRTUAL STORAGE REFERENCES; 9 READS; 0 WRITES. 

DASD121 DETACHED 
DAS D 122 DETACHED 
DASD 120 DETACHED 



**** ** ***** *** ******************************* ** ** ************* 

DO YOU WANT A VRSTEC PLOTTER SMOOTH COPY OF THE 
T F.E DISSPLA METAFILE THAT YOU JUST CREATED? 

(Y OR N) 

**** ******* *** ********* *** * ********************** ******** ** * * * 
I (120) E/O 

DASD 001 LINKED R/O ; E/W BY 'MVS; R/O BY 0700P 
Z (00 1) E/C - OS 
DASD 001 DETACHED 
CREATING NEW FILE: 

CREATING NEW FILE: 

PUN FILE 8317 TO MVS COPY 001 NOHOLD 

DASD 120 DETACHED 



**** ************ ** ********* ********************** ******** ** * * * 
YOUE GEAPH(S) CAN BE PICKED UP AT THE COMPUTER CENTER. 



*********** ********* ****************************************** 
**** *************** ******** ****************** ***************** 



TEE GRAPH (S) WILL BE ADDRESSED TO "POP (USER ID) " 



DO YOU WANT TO 




2 . 



RUN OPTS YSX AGAIN 
RUN THE PLOT PROGRAM US IN 
(TO PLOT ANOTHER TYPE OF 
QUIT 



G THE SAME MATRICES? 
SYSTEM (OPEN/CLOSED)) 



ENTER 1 , 2 OR 3 
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**************** ** ********* ********************** ************* 



3 

**************** *** ******** ****************** **** ******** ** * * * 

HAVE A GOOD DAY! ! 

*************************** *********************************** 
R; T=19. 00/31. 53 09: 37:38 
record off 

END RECORDING OF TERMINAL SESSION 

The graphical output generated by this example follows 
as figure 3.4. 
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FILTER + REGULATOR CLOSED LOOP SYSTEM 

EXAMPLE 5.1 

LINEAR OPTIMAL CONTROL SYSTEMS 



© 
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Figure 3.4 Filter plus Regulator Closed-loop Time Response 






IV. CONCLUSIONS AND RECOMMENDATIONS 
A. CCNCIO SIONS 

As an ultimate evaluation of the computational abilities 
of OFTSYSX, the program was tested using an 82 X 82 matrix 
of aircraft longitudinal motion equations for the X-29A 
experimental forward-swept wing Fighter aircraft prototype, 
provided by NASA-Edwards. 

For a system of equations of this magnitude, all program 
arrays were re-dimensioned, and a 2-megabyte virtual machine 
size was required. 

The graphical time response curves generated from the 
X-29A longitudinal system matrix follow as figures 4.1, 4.2, 
4.3 and 4.4. The accuracy of these time response curves is 
mixed. All of the states shown have the correct waveforms, 
but differ in a scale factor of approximately times 10.0. 
Unfortunately the data supplied by NASA was not explicit 
regarding hew the control input was applied, and whether any 
additional gains were used in their simulation of the 
system. Time constrains did not allow the clarification of 
these problem areas, however the results of the X-29A longi- 
tudinal system are encouraging and should be a topic of 
further research. 
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29A LONGITUDINAL SYSTEM 
82 X 02 MATRIX 
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Figure 4.1 X-29A Longitudinal Time Response 
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Figure 4.2 X-29A Longitudinal Time Response 
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figure 4.3 X-29A longitudinal Time Response 
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figure 4.4 X-29A Longitudinal Time Besponse 



It is hoped that control system instructors will 
encourage their students to use this interactive graphical 
time response program for all applicable class projects; and 
that its enhanced capabilities will stimulate both interest 
in and research on basic systems control problems, as well 
as more advanced designs. 

E. R ECC EMENDATIONS 

Eased on the results of this thesis, three areas emerged 
as possibilities for further research and study: 

1 . P ro gra m Availabili t v 

Ihe use of OPTSYSX and similar design programs 
should be encouraged in all undergraduate and graduate level 
courses involved in the analysis and design of control 
systems. Toward this end, it is recommended that the OPTSYS 
time response programs be placed in the non-IMSL library of 
subroutines, making it easily available to all potential 
users. 

2 . P rogra m Memory R equ i re m ents 

Nhen configured for large matrix operations (98X98) 
the CPTSYSX program requires over 2 megabytes of virtual 
memory. Virtual machines of. this size are not normally 
available to a user. The memory usage for matrix storage is 
a possible area of improvement in the efficiency of the 
CPTSYSX program design. All matrices calculated in OPTSYSX 
(except DUMMY matrices) are still available when the run is 
finished. This simplifies program operation but uses an 
excessive amount of memory. Memory usage should be studied 
and program modifications should be made to reduce the 
excessive memory requirements. 
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3* Further Modifications 

Program sequencing during optimal filtar synthesis 
should be modified. At the present time a regulator must be 
designed or supplied when a filter is designed. Various 
test runs indicate that this problem can be overridden if 
the number cf controls (Nc) is given as zero, but this is 
not a viable solution for systems which use a driving 
function . 



APPENDIX A. OPTSYS EXEC LISTING 



STE ACE OFF 

*********** ******* ********* ****************** ** ** ******** ** * 
* 

* THE OPTSYS EXEC 

* 

* CONTROLS THE OPTSYSX, OPTCALC AND OPTPLOT 

* TC DETERMINE THE TIME RESPONSE OF A SYSTEM. 

* 

* BY H. A. DIEL 

* VERSION 1.0 16 JULY 1984 

* 

*#*#£*****# *** * * ********* ********************** ******** ** 
**** * * * * * * * *** * * ** *************************** ** ******** ** ** 

* 

* CHECK FOR USEE’S VM SIZE = > THAN 1 MEGBYTE 

* 

*** ******** *** ****** ****** ******************* ** * * * ******* *** 
VMSIZE 

SIF SEC GE 1024 SGOTC -TWO 
CLP. SCEN 

SBEGTYPE -ENDTHREE 

*********** *** ********************************************** 



YOU MUST HAVE A 1M OR LARGER VIRTUAL MACHINE 
TO RUN THIS OPTSYS PROGRAM 



TO 


DEFINE A 


1 M 


VIRTUAL 


MACH IN 


E: 


DEFINE 


STORAGE 


1 M 




[PRESS 


ENTER) 


I CMS 








'PRESS 


ENTER 


OPTSYS 








PRESS 


ENTER) 


FOR 


SYSTEMS 


LA 


RGER THAN 32 X 


32 


OETAIN 


A LISTING 


OF THE OPTSYS PROGRAM 



AND FOLLOW INSTRUCTIONS CONTAINED IN THE LISTING. 

*********** *** * * ** ********* ****************** ** £**$£#**$*#*# 

-ENDTHREE 

SEXIT SRC 

-TWO 

CLRSCRN 

SBEGTYPE -ENDZERO 



THE OPTSYS EXEC CONTROLS 



1 . 

2 . 

3. 



OFTSYSX 

OPTCALC 

OPTPLOT 



FORTRAN 

FORTRAN 

FORTRAN 



EACH PROGRAM PASSES INFO 
THROUGH A DATA FILE WRIT 
CASE, THESE FILES ARE "0 
THE SIZE OF THESE FILES 
CAN USE ABOUT 20X OF THE 
ENSURE THAT SUFFICIENT D 



A TRIO OF PROG 

SYSTEM ANALYSIS 
CALCULATE TIME 
DISSPLA PLOTTIN 

RMATION TO THE 
TEN TO THE USER 
ETM AT DATA" AND 
VARY WITH THE S 
USERS DISK SPA 
ISK SPACE IS AV 



RAMS: 



RES PONS E) 

G ROUTINE) 

NEXT PROGRAM 
S DISK. IN THIS 
"OPTPLOT DATA". 
YSTEM ORDER, AND 
CE. THEREFORE 
AILABLE. 
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tt * -fr * X 



- TYPE "E" TO EXIT, ANY OTHER ENTRY TO CONTINUE - 
-END ZERO 
SREAD VARS SANS 
SIF .SANS EQ .E SEXIT SRC 

**************** ** ********* ********** *********************** 

3je * 

* Allow THE USE OF AN OID "OETM AT DATA A1 M * 

* * 

##*##$#**** *** * * #$***#*$$****##****:*#$: 3$C# J$C 

RENAME OPTMAT DATA A1 OPTSY S DATA A1 
SIF SRC KE 0 S GO TO -FIRST 
RENAME OPTSYS DATA A1 OPTMAT DATA A1 
CLRSCRN 

SBEGTYPE -ENDONE 



## # # $$$$$$$ ## 3$C 3$C 3$C ## # $£$$$$$$ $ £ 3$C £ 

YOU HAVE A DATA FIIE NAMED ’OPTMAT DATA’ ON YOUR A 
DISK THAT WAS PREVIOUSLY GENERATED BY THE OPTSYS 
PROGRAM AND CONTAINS THE F, G, H, GAMMA, A AND B 
MATRICES FROM THAT RUN. 

IF YOU WOULD LIKE TO USE THESE SAME MATRICES FOR 
THIS RUN, THE OPTSYS PROGRAM WILL READ IN THE 
DESIRED DATA AT THE APPROPRIATE TIME, 

IF YOU TYP E (Y) ES. 

ANY OTEER INPUT WILL RESULT IN THAT FILE BEING ERASED! 

************************************************************ 

-ENDCNE 

SREAD VARS SANS 

SIF .SANS EC . Y SGOTC -ONE 

-FIRST 

*************************** ********************************* 



$ ♦ 

* ERASE THE OLD "OPTMAT DATA A1" DATA FILE * 

* PLACE "000 0" IN TEE NEW "OPTMAT DATA FILE" * 

* TO ACT AS A FLAG FOR OPTS YSX AND OPTCALC * 

* * 



ERASE OPTMAT DATA A1 
SSTACK 000 0 

FILESTCK OPTMAT DATA A1 F 80 1 

-ONE 

-THIRD 

CLRSCRN 

SBEGTYPE -ENDFOUR 



3Qc 3$C 3$C 3fC3fC3fC3fC3^C3{C3fC3fC3fc^C jjt^C 3jC 3}C3}C3$C3}C3jC3}c3jC3jC35C35c 3$C 3$C 3$C 3$C 5$C 3$C jjc^c ^C^C3jc ^ 3$C 2$C 3$C 3$C 3$C 

DC YOU WANT THE NUMERICAL OUTPUT FROM CPTSYSX TO GO 
TO YOUR TERMINAL SfCREEN) OR TO A D (ISK) FILE? 

(S OR D) 

**** ********** ***** ******* * ****************** **** * ***** ** *** 

-ENDFOUR 

SREAD VARS SANS 

SIF .SANS EQ . S SGOTC -FOURTH 
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SIF .SANS EC . D SGOTC -FIFTH 
CLRSCEN 

SBEGTYPE -ENDFIVE 



********************************************* t **** *********** 
YOU HOST ANSHER S(CREEN) OR D (ISX) . 

-ENDFIVE 
CP S1EEP 3 SEC 
SGOTO -THIRD 
-FOURTH 
CLRSCEN 

SBEGTYPE -ENDSIX 



*** ******** * ** **** * * ** ***** ********************** ****** ** *** 
OUTPUT NIL! COME TO YOUR TERMINAL SCREEN. 
************** ********************************************** 

-ENDSIX 

CP SLEEP 1 SEC 

STYPE LOADING OPTSYS 

FILEDEF 06 TERM (RECFM FA BLKSIZE 133 

FILEDEF 8 DISK OPTPLOI DATA A1 f PER M 

FILEDEF S DISK OPTMAI DATA A1 (PERM 

GLOBAL TXTLI3 FORTMOD2 MOD2 EEH IMSLDP NONIMSL 

LOAD OPTSYSX (START 

SGOTC -FIVE 

-FIFTH 

n p crtM 

SBEGTYPE -ENDSEVEN 



**** ******* *** ****** ******* ********************************* 

OUTPUT WILL GO TC DISK FILE 'OUTPUTX LISTING A1’ 

**** * ****** *** ********* *** *********************** ******** *** 

-ENDSEVEN 
CP SLEEP 1 SEC 

STYPE LOADING OPTSYS 

FILEDEF 06 DISK OUTPUTX LISTING A1 

FILEDEF 8 DISK OPTPLCT DATA A1 (PERM 

FILEDEF 9 DISK OPTMAT DATA A1 (PERM 

GLOBAL TXTLIB FORTMOD2 MOD2 EEH IMSLDP NONIMSL 

LOAD OPTSYSX (START 

-FIVE 
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CLRSCEN 

SBEGTYPE -ENDEIGHT 



*********** ******* ******************* ** ****** **** * ********** 
IF YOO ARE DISSATISFIED WITH THE RESULTS 
THUS FAR AND WOULD LIKE TO EXIT TO CMS, 

-TYPE ’ Y’ TO EXIT- 

(ANY OTHER INPUT TO CONTINUE) 

*** ***************** ******* ****************** ** ************* 

-ENDEIGHT 

SREAD VARS SANS 

SIF .SANS EC .1 SEXIT SRC 

-SIXTH 

*********** *** ** ** ********* **************** ** **** *********** 



* * 

* CHECK FOR DATA IN THE FILE "OPTMAT DATA " * 

* BEFORE LOADING OPTCALC * 

* * 



************ ** ************* ********************** *********** 

FINDSTAK OPTMAT DATA A1 C0 1 0 LIM002 ALL G ECU PI 

SREAD VARS SA1 S A2 

SIF . SA2 EQ .0 SEXIT SRC 

STYPE LOADING OPTCALC 

FILEDEF 8 DISK OPTPLCT DATA A1 (PERM 

FILEDEF 9 DISK OPTMAT DATA A1 (PERM 

GLOBAL TXTLIB FORTMOD2 MOD2EEH IMSLDP NOMIMSL 

LOAD OPTCALC (START 

-SEVENTH 

CLRSCEN 

SBEGTYPE -ENDNINE 



************************************************************ 



IF YOU ARE DISSATISFIED WITH THE RESULTS 
THUS FAR AND WOULD LIKE TO EXIT TO CMS, 

-TYPE ' Y' TO EXIT- 

(ANY OTHER INPUT TO CONTINUE) 

*********** *** ******************************************* ** * 

-ENDNINE 

SREAD VAES SANS 

SIF .SANS EQ .Y SEXIT SRC 

FILEDEF 8 DISK OPTPLCT DATA A1 (PERM 

EXEC DISSP1A OPTPLOT 

-EIGHTH 

**************** ** ***************************************** 

* 

* CHECK FOR FILE "DI SSPLA METAFILE A4" ON 

* THE USER'S DISK BEFORE GOING TO DISSPOP 

* 

*********** *** ******************************* ** ************ 
RENAME DISSPLA METAFILE A4 OPTSYS METAFILE A4 
SIF SRC KE 0 S GO TO -TENTH 

RENAME OPTSYS METAFILE A4 DISSPLA METAFILE A4 
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# # # * * # 



CLRSCEN 

SBEGTYPE -ENDTEN 



************************************************************ 

DC YOU KANT A VRSTEC PLOTTER SMOOTH COPY OF THE 
THE DISSPLA METAFILE THAT YOU JUST CREATED? 

(Y OR N) 

************** ************************** *** ** *************** 



-ENDTEN 

SREAD VARS SANS 

SIF .SANS EQ . Y SGOTC -NINTH 

SIF .SANS EC . N SGO TO -TENTH 

CLRSCEN 

SBEGTYPE -ENDELEVEN 



*********** *** * * ** ***** *** ********************************** 
YOU MUST ANSWER Y(ES) OR N (0) . 

************** ********************************************** 

-END ELEVEN 
CP SLEEP 4 SEC 
SGOTC -EIGHTH 
-NINTH 

EXEC DISSPCP VRSTEC 
CLRSCEN 

S3EGTYEE -ENDTWEIVE 



**** ******* *** * * ** ********* ********************** ******** ** * 

YOUR GRAPH (S) CAN EE PICKED UP AT THE COMPUTER CENTER. 

THE GRAPH (S) Will BE ADDRESSED TO "POP (USER ID)". 

*********************** ************************************* 

-ENDTWELVE 
CP SLEEP 5 SEC 
-TENTH 
CLRSCEN 

SBEGTYPE -ENDTHIRTEEN 



********************************************* ** ** *********** 
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DO YOU WANT TO 



1. RON OPTSYSX AGAIN 

2. * RON THE PLOT PROGRAM USING THE SAME MATRICES? 

(TC PLOT ANOTHER TYPE OF SYSTEM (OPEN/CLOSED)) 

3. QUIT 

ENTER 1, 2 OR 3 

*********** **************** ********************************* 

-ENDTEIRTEEN 

BREAD VARS SANS 

SIF .SANS EQ . 1 SGOTO -THIRD 

SIF .SANS EQ .2 SGOTC -SIXTH 

pT per dm 

SEEGTYPE -ENDGOODBY 



**** *********************** ****************** ***** ********** 



HAVE A GOOD DAY! ! 



-ENDGCODEY 
CP SLEEP 3 SEC 
CLRSCEN 
SEXIT SRC 
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TFESE ARE EXAMPLES CF SEVERAL POSSIBLE METFODS GF ARRAY G EAERAT1GN 
WlTHIA SUE f CUT INE SETUP. TFE "GAM" ARRAY WAS SET TO ZERO STALE AO 
41 AO IS E " WAS PRESENT, A AD THE AON-ZERG ELEAEATS OF THE ,, G ,, ARRAY WERE 
EXPLICITLY CEFINED. THEY COULD ALSO BE REAC FRCM FILES AS ABOVE. 
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140 FORMAT (//,5X,49H ********** INCONSISTENT PSD INPUT FLAGS ******** 
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IF (IANS.EQ.INO) I SAB = 
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8R0LTINE RDMAT — REAOS THE F, G, H, GAP « A AND 6 MATRICES FROM 
TRICES FROM THE DATA FILE OPTMAT ON FILECEF S. 
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